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ABSTRACT 


At  cyclic  loading  frequencies  below  that  of  the  transitional  frequency  level,  the 
elevated  temperature  fatigue  crack  growth  process  in  Alloy  718  is  viewed  to  be  fully 
environment-dependent.  Of  all  the  crack  growth  stages,  this  stage,  while  is  the  most 
critical  in  high  temperature  applications  due  to  its  highly  accelerated  crack  growth  rate, 
is  the  least  studied  or  understood.  The  objective  of  this  study  is  to  focus  on  the 
understanding  of  the  controlling  mechanisms  of  this  environment-dependent  crack  growth 
stage  in  Alloy  718  in  order  to  develop  the  ability  to  predict  the  crack  growth  performance 
under  different  loading  and  environment  conditions.  For  this  purpose,  two  major  studies 
have  been  carried  out  here.  The  first  was  to  provide  evidence  of  the  existence  of  the  fully 
environment-dependent  stage  in  which  the  crack  growth  rate  is  viewed  to  be  equal  to  the 
oxygen  penetration  rate  at  the  crack  tip.  This  has  been  achieved  by  indirectly  estimating 
the  depth  of  oxygen  diffusion  at  the  crack  tip  of  compact  tension  specimens  tested  at 
650°C  under  cyclic  loadings  with  and  without  hold  time  durations  at  minimum  load  level. 
As  a  result,  the  relationship  between  the  intergranular  oxygen  diffusion  rate  and  the  value 
of  the  stress  intensity  factor  has  been  established.  This  relationship,  when  integrated  over 
the  cycle  effective  oxidation  time,  results  in  a  closed  form  solution  describing  the 
environment-dependent  fatigue  crack  growth  rate.  A  good  agreement  was  obtained  when 
a  comparison  was  made  between  the  results  of  this  solution  and  the  corresponding 
experimental  results. 

The  second  study  was  to  establish  a  micromechanical  based  quantitative  model  to 
predict  the  environmentally-dominated  crack  growth  performance.  This  model  wus  bused 
on  the  concept  of  the  intergranular  two-stage  crack-tip  oxidation  mechanism.  According 
to  this  mechanism,  the  oxygen  partial  pressure  controls  the  preferential  formation  of  the 
oxide  layers  at  the  crack  tip,  and  it  suggests  that  intergranular  oxidation  under  air 
environment  occurs  in  two  stages  with  the  formation  of  a  Ni-Fe  rich  oxide  followed  by 
that  of  the  more  protective  film  formed  by  Cr203  oxide.  Thus,  the  reduction  of  the  grain 
boundary  ductility  along  the  grain  boundary  fracture  path  depends  on  the  rate  of 
formation  of  the  chromia  layer  in  relation  to  the  build-up  of  other  oxide  types  at  the  crack 
tip.  The  determination  of  the  amount  of  Cr203  depends  on  the  amount  of  both  oxygen 
diffused  along  the  affected  grain  boundary  and  chromium  transported  via  a  mobile 
dislocation  network.  The  reduction  in  grain  boundary  ductility  due  to  oxidation  is 
balanced  by  considering  the  effective  strain  at  the  crack  tip  resulting  from  the  external 
loading.  This  balance  defines  the  fracture  criterion  of  the  model  and  permits  the 
calculation  of  the  crack  advance  per  cycle,  The  model  was  then  used  to  predict  the  crack 
growth  rate  in  Alloy  718  at  650  °C  for  two  different  loading  conditions,  and  u  good 
agreement  was  achieved  when  the  results  obtained  from  the  model  compared  with  those 
obtained  experimentally. 

The  significant  contribution  of  this  three-year  study  is  that  it  provides  the  jet 
engine  manufactures  with  a  mechanical  approach,  parallel  to  alloying,  by  which  the 
intrinsic  resistance  of  jet  engine  components  to  environment-dependent  crack  growth 
process  can  be  controlled. 
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CHAPTER  1 


INTRODUCTION 


Alloy  718*  is  a  precipitation  strengthened  nickel-base  superalloy  with  a  excellent 
resistance  to  creep  and  stress  rupture  up  to  650°C.  Because  of  its  good  mechanical 
properties,  weldability,  fabricability  and  cost  effectiveness,  Alloy  718  is  widely  used  in 
aircraft  engines,  nuclear  structures,  oil  fields  and  gas  production  industries.  Damage 
tolerance  methodologies  needed  in  some  of  these  applications  require  accurate  predictions 
of  the  crack  growth  behavior  under  typical  service  conditions.  For  this  reason,  and  as 
part  of  material  characterization,  extensive  studies  [1-31,37,47,49,50,69,82]  have  been 
carried  out  in  the  past  three  decades  investigating  the  influence  of  vaiious  operating 
parameters  on  the  crack  growth  mechanisms  under  different  loading  and  environment 
conditions.  In  order  to  recognize  the  role  of  each  of  these  parameters  and  to  establish  the 
objectives  of  this  research  program,  an  attempt  will  be  made  in  this  chapter  to  review 
several  related  studies.  The  parameters  which  will  be  analyzed  in  this  nr  iew  are 
temperature,  frequency,  wave  shape,  hold  time,  microstructure  and  environment. 

Effects  of  Temperature 

A  comprehensive  evaluation  of  the  effect  of  temperature  on  the  crack  growth  in 
Alloy  718  is  that  of  Weerasooriya  et  al  [  1  ].  Constant  AK  tests  at  R=0.1  were  conducted 

1  Alloy  718  is  a  generic  name  of  Inconel  718  which  is  tire  trademark  of  Huntington  Alloys  Inc.. 


at  Kn,M  =  27.8  and  40.0  MPaVm  for  frequencies  ranging  from  0.001  to  1.0  Hz.  Results 
of  the  tests  at  =  40  MPaVm  in  the  form  of  the  fatigue  crack  growth  rate  versus 
reciprocal  absolute  temperature  are  shown  in  Fig.  1.1.  The  frequency  effects  are 
discussed  in  the  next  section.  For  temperatures  between  room  temperature  and  650°C, 
the  growth  rate  increases  with  increase  in  temperature.  The  rate  of  increase  in  the  lower 
temperature  regime  (c)  is  limited  and  is  attributed  to  a  combination  of  the  change  in 
modulus  and  plastic  deformation  characteristics  of  the  material  with  temperature.  The 
fracture  mode  is  transgranular  and  the  behavior  is  characterized  as  cycle  dependent.  At 
some  transition  temperature,  depending  upon  frequency,  the  growth  rate  increases 
significantly  with  further  increase  in  temperature.  In  this  region  (b)  the  fracture  is  a 
combination  of  transgranular  and  intergranular  and  the  behavior  is  characterized  as  mixed 
mode.  As  can  be  observed  in  Fig,  1.1,  the  temperature  at  which  this  transition  occurs  is 
higher  for  higher  frequencies.  At  the  highest  temperatures,  the  behavior  of  this  material 
becomes  purely  intergranular,  the  mode  is  characterized  as  time  dependent,  and  the  rate 
of  change  of  growth  rate  with  temperature  is  maximum.  This  region  is  denoted  by  (a) 
in  Fig.  1.1.  The  transition  temperature  from  mixed  mode  to  time  dependent,  denoted  in 
this  figure  as  T^,  again,  depends  on  test  frequency  and  has  a  higher  value  for  higher 
frequencies.  Similar  observations  of  the  fatigue  crack  growth  rate  behavior  as  a  function 
of  frequency  where  the  mode  of  fracture  and  the  rate  of  crack  growth  can  be  correlated 
are  presented  in  the  next  section. 

The  effects  of  temperature  on  the  fatigue  crack  growth  rate  can  also  be  observed 
when  Alloy  718  is  tested  in  vacuum  environment.  However,  the  fatigue  crack  growth  rate 
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temperature  on  crack  growth  rate  at  =  40  MPaVrn. 


only  increases  slightly  with  Increase  in  temperatures  compared  to  the  changes  observed 
in  air  [2], 


Effects  of  Frequency 

The  effects  of  cyclic  frequency  on  the  fatigue  crack  growth  behavior  of  Alloy  718 
has  been  the  subject  of  many  researches  [3,4, 5-8],  James  [3]  studied  the  frequency  effects 
for  Alloy  718  at  538°C  over  the  range  of  1.39xl0‘3  Hz  to  6.67  Hz,  the  results  are  shown 
in  Fig.  1,2,  It  can  be  seen  from  this  figure  that  the  fatigue  crack  growth  rates  increase 
as  cyclic  frequencies  decrease.  This  trend  was  also  noted  by  Clavel  and  Pineau  [4]  in 
their  work  at  550°C  over  the  range  of  5xlO-3  Hz  to  20  Hz.  Crack  growth  rates  in 
laboratory  air  at  a  temperature  of  650°C  are  also  found  to  depend  strongly  on  frequency. 
Increasing  the  frequency  generally  results  in  a  decrease  in  crack  growth  rate  per  cycle. 
For  example,  Weerasooriya  [6,7]  has  shown  that  at  two  constant  AK  conditions  the 
behavior  ranges  from  purely  time  dependent  at  frequencies  below  0.01  Hz  to  purely  cycle 
dependent  at  frequencies  in  excess  of  10  Hz.  His  work  also  show  that  different 
mechanisms  of  crack  growth  are  activated  in  each  regime;  striation  formation  is  generally 
observed  in  the  cycle-dependent  regime,  intergranular  crack  growth  is  dominant  in  the 
time-dependent  regime,  and  a  mixed  transgranular-intergranular  mode  of  crack  extension 
is  observed  in  the  mixed  regime.  Similar  observations  of  the  frequency  dependence  have 
been  documented  by  Florcen  and  Kane  [5]  and  Pedron  and  Pineau  [9].  Although  the 
magnitude  of  the  change  in  growth  rate  with  decrease  in  frequency  varies  depending  on 
the  alloy  microstructure,  temperature,  and  frequency,  the  pattern  of  behavior  is  widely 
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The  effects  of  cyclic  frequency  on  crack  growth  rate  at  538°C  in  air. 


noted  in  all  Nickel-based  superalloys  and  is  generally  attributed  to  decreasing 
environmental  resistance  with  increasing  exposure  time  as  cyclic  frequency  decreases. 

While  most  of  frequency  studies  were  conducted  in  air  environment,  Floreen  and 
Kane  [5]  studied  the  fatigue  crack  growth  behavior  of  Alloy  718  at  650°C  over  the  range 
of  0.01  Hz  to  1  Hz  in  air  and  helium  (99.995%  purity)  as  well.  In  helium  environment, 
as  the  frequency  decreased,  the  fatigue  crack  growth  rate  only  slightly  increased,  while 
in  air  environment,  the  fatigue  crack  growth  rate  increased  substantially  as  frequency 
decreased. 

Similar  to  other  high  strength  superalloys,  the  frequency  effects  for  Alloy  718  are 
only  observed  at  elevated  temperatures  and  not  at  room  temperature  [8'J.  Further,  the 
cycle  dependent  (high  frequency)  growth  rate  at  high  temperature  corresponds  closely 
with  the  rate  observed  at  room  temperature  as  well  as  those  observed  in  vacuum, 
indicating  a  lack  of  environmental  contributions  to  the  growth  rate  process  when  the 
transgranular  fatigue  mode  dominates  [7]. 

Effects  of  Wave  Shape 

Although  the  behavior  of  crack  growth  in  Alloy  718  as  a  function  of  frequency 
is  widely  noted  arid  accepted,  this  behavior  when  more  complex  wave  shapes  are  involved 
is  not  as  well  understood  nor  are  the  results  consistent  among  the  various  investigators. 
There  are  a  number  of  observations  involving  other  wave  shapes  which  illustrate  the 
variable  nature  of  the  elevated  temperature  crock  growth  behavior.  James  [10]  found  that 
the  crack  growth  rate  at  538°C  in  Alloy  718  is  significantly  faster  for  a  square  wave 


form  than  for  a  sawtooth  wave  form  of  the  same  cyclic  duration  when  the  loading 
frequency  was  much  lower  than  the  unloading  frequency.  In  these  experiments,  the  hold 
time  contribution  in  the  square  wave  appears  to  dominate  the  cyclic  contribution,  even  for 
a  low  loading  frequency  in  the  sawtooth  wave.  Somewhut  conflicting  results  were 
obtained  by  Clavel  and  Pineau  [4]  who  examined  the  wave-shape  effects  on  the  fatigue 
crack  growth  rate  using  square  and  triangular  wave  shapes  at  550°C.  Three  cyclic  wave 
shapes  were  chosen  (shown  in  Fig,  1.3)  for  differentiating  the  effect  of  loading  and 
unloading  rate  as  well  as  the  effect  of  hold  time  at  maximum  and  minimum  load. 
Wave-shape  3  cot  responds  to  a  frequency  of  2  Hz,  whereas  wave-shape  2  leads  to  a 
frequency  of  0.05  Hz.  Wave-shape  1  has  the  same  loading  and  unloading  rate  as 
wave-shape  3  but  superposes  a  10  s  hold  time  at  the  maximum  load  level  and  at  the 
minimum  load  level,  and  its  cyclic  duration  is  the  same  as  that  of  wave-shape  2.  Their 
results  showed  that  the  fatigue  crack  growth  rate  of  wave-shape  2  is  higher  thun  those  of 
the  other  two  wave  shapes,  and  the  fatigue  crack  growth  rate  of  wave-shape  1  is  only 
slightly  higher  than  that  of  wave-shape  3.  They  inferred  that  a  hold-time  of  10  s  both  at 
the  maximum  and  at  the  minimum  load  has  no  significant  effect  on  the  fatigue  crack 
growth  rate,  and  that  the  main  parameter  is  not  the  frequency  per  se  but  the  loading  and 
the  unloading  rate.  As  already  pointed  out  by  Bathias  and  Pclloux  [11]  the  slip  character 
strongly  affects  the  aspect  of  die  fracture  surfaces,  In  fact,  their  fractogTaphic  examination 
showed  that,  at  550°C,  the  cracks  propagate  by  a  transgranular  mode  when  the  frequency 
is  higher  than  0,5  Hz,  while  at  lower  frequency  (0.05  Hz)  the  crack  growth  occurs  by  a 
mixture  of  transgranular  and  intergranular  fracture.  The  well  known  fact  for  this  material 
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is  that  homogeneity  in  plastic  deformation  promotes  transgranular  fracture  mode  while 
inhomogeneity  in  plastic  deformation  favors  intergranular  fracture  mode  [12,13-17],  This 
subject  will  be  discussed  in  more  detail  later.  Obviously,  for  wave-shape  1  the  loading 
and  unloading  rates  promoted  the  homogeneity  in  plastic  deformation  which  generally  has 
less  time-dependent  sensitivity,  thus  the  hold  times  after  the  loading  and  unloading 
portions  are  not  expected  to  produce  a  large  influence  on  the  fatigue  crack  growth  rate 
which  is  why  the  fatigue  crack  growth  rate  of  wave-shape  1  is  only  slightly  higher  than 
that  of  wave-shape  3.  For  wave-shape  2,  the  loading  and  unloading  rates  produced  the 
inhomogeneity  in  plastic  deformation  near  crack  tip  which,  in  turn,  favors  the  occurrence 
of  intergranular  fracture,  and  certainly  the  fatigue  crack  growth  rate  is  much  higher  than 
that  of  wave-shape  1. 

Ashbaugh  [18]  conducted  a  series  of  tests  at  649°C  using  triangular  wave  shapes 
on  CT  specimens  and  R*(),l.  In  addition  to  symmetric  wave  shapes  with  frequencies 
covering  the  range  from  0.01  Hz  to  10  Hz,  asymmetric  wave  shapes  were  used  in  which 
the  rise  and  fall  times  had  different  values  corresponding  to  the  baseline  frequencies  of 
the  symmetric  wave  form  tests.  For  various  combinations  of  rising  and  falling 
frequencies  from  1  Hz  to  0.0 1  Hz,  growth  rates  per  cycle  were  essentially  equivalent  to 
those  obtained  using  symmetric  wave  shapes  where  the  frequency  was  that  of  the  rising 
portion  of  the  asymmetric  wave.  This  observation  implies  that  the  fatigue  damage  occurs 
almost  exclusively  during  the  rising  load  portion  of  the  fatigue  cycle  At  the  temperature 
and  .frequencies  at  which  these  experiments  were  conducted,  the  fracture  mode  is  a  mixed 
mode  combination  of  intergranular  and  transgranular  [7], 
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Effects  of  Hold  Time 


It  is  widely  reported  that  the  addition  of  a  hold  time  at  maximum  load  to  a  fatigue 
cycle  tends  to  increase  the  crack  growth  rate  per  cycle  in  Alloy  718  as  well  as  in  a 
number  of  other  nickel-base  superalloys  [4,19-28].  Fig.  1.4  shows,  for  example,  results 
from  the  work  of  Pedron  and  Pineau  [24]  involving  Alloy  718  at  650°C  where  hold  times 
are  added  to  baseline  fatigue  cycles.  This  results  in  an  increase  in  growth  rate  per  cycle. 
In  all  three  cases  it  can  be  seen  that  the  behavior  at  long  hold  times  is  essentially  purely 
time  dependent,  that  is,  growth  rate  increases  linearly  with  increase  in  hold  time.  There 
are,  however,  a  number  of  observations  which  do  not  show  this  same  trend.  Shahanian 
and  Sadananda  [20]  show  that  at  760°C,  the  addition  of  hold  times  of  0.1,  1,  or  10  min. 
does  not  increase  the  crack  growth  rate  in  Alloy  718.  The  results  of  this  investigation 
also  showed  that  growth  rates  obtained  under  constant  load  range  (increasing  AK) 
conditions  depend  on  load  amplitude,  thereby  demonstrating  the  inadequacy  of  AK  as  a 
correlating  parameter  at  this  temperature.  At  593°C,  VanStone  et  al.  [28]  showed  that 
Alloy  718  still  exhibits  an  increase  in  growth  rate  with  increase  in  hold  time,  but  the 
effect  is  not  as  pronounced  as  at  650°C. 

There  are  several  instances  where  hold  times  retard  crack  growth  instead  of 
accelerating  it  and  even  contribute  to  total  crack  arrest.  This  has  been  observed  when 
hold  times  occur  at  other  than  maximum  load  [22],  In  these  cases,  either  significant  creep 
occurs  at  the  test  temperature,  thereby  blunting  the  crack  tip,  or  the  stress  intensity  during 
hold  is  generally  less  than  the  creep  crack  growth  threshold. 

Several  observations  have  been  made  in  Alloy  718  at  649°C  of  the  effect  of  a  hold 
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The  effects  of  hold  time  on  crack  growth  rate  at  650°C  in  air. 


time  at  minimum  load  on  the  crack  growth  rate.  In  these  studies,  the  minimum  load  has 
been  below  the  threshold  for  crack  growth  under  sustained  load,  thereby  producing  no 
contribution  in  a  linear  summation  model  prediction.  Nicholas  and  Wecrasooriya  [23] 
found  no  effect  of  hold  times  between  50  and  500  s  under  constant  K  conditions. 
Shahinian  and  Sadananda  [22]  and  Sadananda  and  Shahinian  [26],  on  the  other  hand, 
observed  a  definite  increase  in  growth  rate  due  to  the  addition  of  a  hold  time  at  minimum 
load.  They  attributed  this  acceleration  of  growth  rate  to  the  additional  time  available  for 
oxidation.  Similar  observations  have  been  made  by  Diboine  and  Pincau  [25]  where  hold 
times  at  minimum  load  between  3  and  3000  s  showed  an  increase  in  growth  rate 
compared  to  that  obtained  under  continuous  cycling.  In  comparing  these  apparently 
contradictory  results,  Ghonem  et  al.  [17]  have  noted  that  the  frequency  of  the  cyclic 
loading  was  different  in  the  three  cases  cited.  In  the  work  of  Nicholas  and  Weerasooriya 
[23],  the  cyclic  frequency  was  1  Hz  which  produces  tTansgranular  crack  growth.  In  [22] 
and  [25],  the  frequencies  were  0.1  and  0.05  Hz,  respectively,  both  of  which  produce 
intergranular  crack  growth.  Thus,  the  effect  of  frequency  is  important  in  this  situation 
because  of  the  differences  in  mechanisms  resulting  from  different  rates  of  diffusion  of 
oxygen  into  the  grain  boundaries.  As  discussed  previously,  the  strain  rate  for  the  1  Hz 
cycle  must  promote  homogeneity  in  plastic  deformation  near  the  crack  tip,  and  this  kind 
of  deformation  mode  has  less  time-dependent  sensitivity,  thus  the  hold-time  effect  does 
not  appear.  However,  strain  rates  for  0.17  Hz  and  0.05  Hz  certainly  favor  the 
inhomogeneity  in  plastic  deformation  near  the  crack  tip,  which  will  enhance 
time-dependent  behavior  and  thus  will  show  strong  hold-time  effects. 
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Saturation  effects  of  the  hold  time  at  minimum  load  level  on  the  fatigue  crack 
growth  rate  have  been  noted  by  both  Sadananda  and  Shahinian  [26]  and  Diboine  and 
Pineau  [25],  The  former  observed  the  saturation  effect  appears  for  hold  times  about  1 
minute,  the  latter  reported  that  a  saturation  effect  seems  to  occur  for  hold  times  longer 
than  about  1000s. 

If  the  test  temperature  is  not  sufficiently  high,  the  hold-time  effect  may  disappear. 
For  example,  Sadananda  and  Shahinian  [19]  showed  that  at  425°C  the  1-tnin.  hold  time 
period  has  a  negligible  effect  on  the  crack  growth  process.  Furthermore,  they  found  that, 
at  425°C,  hold-times  up  to  10  minutes  still  produced  purely  cycle-dependent  behavior. 
Clearly,  the  effects  of  temperature  and  hold  time  arc  synergistic. 

Very  few  studies  [2,5,9,21,27-30]  have  dealt  with  the  effects  of  hold  time  on 
fatigue  crack  growth  rate  in  an  environment  other  than  in  air.  Ghonem  and  Zheng  [30] 
have  examined  the  effects  of  hold  time  in  vacuum  by  conducting  tests  of  continuous 
cycling  at  a  frequency  of  0.05  Hz  and  cycling  with  a  hold  time  of  300  s  at  maximum  load 
level.  They  found  that  the  fatigue  crack  growth  rate  for  ?  hold  time  test  is  higher  than 
that  for  a  continuous  cycling  test.  Similar  results  have  been  reported  by  Smith  and 
Michel  [21],  Sadananda  and  Shahinian  [29]  and  Vai.Stone  et  al.  [28]. 

Effects  of  Microstructure 

Large  amount  of  observations  have  revealed  that  the  fatigue  crack  growth 
properties  of  Alloy  718  at  elevated  temperatures  are  strongly  dependent  on  their 


microstructure.  In  air  environment,  microstructural  parameters,  particularly  grain  size, 
grain  boundary  morphology  and  precipitate  particle  size,  can  be  altered  to  produce  order 
of  magnitude  changes  in  high  temperature  crack  growth  behavior  [9,10,13,33-43]. 
Thamburaj  et  al  [36]  observed  that  high  temperature  futigue  crack  growth  rate  was  much 
faster  in  fine-grained  structure  than  in  coarse-grained  structure  for  Alloy  718  at  650°C, 
thus  they  deduced  that  more  severe  grain  boundary  oxidation  can  occur  in  relatively  fine 
grained  Alloy  718  due  to  the  greater  grain  boundary  area  available  along  which  rapid, 
stress  assisted  diffusion  of  oxygen  can  occur.  The  work  done  by  James  [10],  Smith  and 
Michel  [37],  Pedron  and  Pineau  [9]  also  showed  that  the  coarse-grained  structure 
improves  crack  growth  resistance,  and  the  necklace  microstructure  renders  die  highest 
fatigue  crack  growth  resistance.  Again,  they  attributed  this  grain  size  effect  on  crack 
growth  rate  to  the  larger  grain  boundary  area  exposed  to  oxygen  penetration  in  the  fined- 
grained  structure  comparing  with  that  in  the  coarse-grained  structure. 

Grain  boundary  morphology  is  another  important  microstructural  factor  that 
ci  ntrols  the  high  temperature  crack  growth  behavior.  Thamburaj  et  al  [36]  have  pointed 
out  that  an  inegular  grain  boundary  morphology  developed  by  grain  boundary  5  phase 
precipitates  will  enhance  creep  ductility  and  resistance  of  high  temperature  futigue  crack 
growth  rate.  Andrieu  [35]  observed  that  the  presence  of  5  precipitates  along  the  grain 
boundaries  leads  to  an  improvement  in  creep  crack  growth  resistance.  Apart  from  the 
factor  of  irregular  grain  boundary  morphology  created  by  8  phase  particles,  he  suggested 
that  this  beneficial  influence  of  8  phase  might  be  related  to  grain  boundary  oxidation 
behavior  either  because  of  the  intrinsic  oxidation  resistance  of  Ni3Nb  phase  or  because 
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of  the  existence  of  oxygen  traps  formed  along  the  5  phase  interfaces. 

While  the  effects  of  grain  size  and  grain  boundary  morphology  are  important,  there 
have  been  a  number  of  investigations  [10,39,41,42]  in  which  indirect  evidence  of  an 
inverse  relationship  between  matrix  y'  and  y"  particle  size  and  the  sensitivity  to 
environmental  embrittlement  has  been  demonstrated.  When  the  matrix  precipitate  is  very 
fine,  dislocations  cut  the  precipitate  particles  instead  of  bypassing  them  and  are  confined 
to  move  in  coplanar  arrays,  causing  planar  pile-ups  and  consequently,  intense  stress 
concentrations  at  grain  ,1 boundaries  and  other  obstacles  to  their  motion.  If  at  the  same 
time,  the  grain  boundary  cohesive  strengths  have  been  weakened  as  a  result  of  oxidation 
effect,  intergranular  cracking  would  be  promoted.  Besides,  it  is  also  thought  that  in 
materials  exhibiting  predominantly  planar  slip  character,  oxygen  diffusion  [43,44]  into  the 
material  can  occur  quite  rapidly  along  planes  of  intense  deformation  at  elevated 
temperatures  and  be  swept  to  the  grain  boundaries  as  a  result  of  a  solute-dislocation 
interaction.  In  the  presence  of  more  homogeneous  slip,  the  number  of  dislocations 
arriving  at  a  grain  boundary  would  be  reduced  and  consequently,  the  stress  concentrations 
and  oxygen  transport  to  the  grain  boundaries  occurring  due  to  a  dislocation  "sweep-in" 
mechanism  will  be  reduced  as  well.  Examined  thin  foils  prepared  from  regions 
immediately  below  fracture  surfaces  by  means  of  transmission  electronic  microscope, 
Thamburaj  et  al  [36]  observed  that  planar  deformation  was  predominant  in  the  case  of 
existing  the  finest  y"  particle  size  distributions,  and  the  poor  resistance  of  fatigue  crack 
growth  in  this  case  is  consistence  with  the  mechanism  discussed  above. 

It  has  been  reported  that  the  overaging  bea.  treatment  or  the  modified  heat 


treatment  (higher  solutioning  temperature  and  longer  aging  time)  for  Alloy  7 1 8  appears 
to  improve  crack  growth  resistance  significantly  in  air  environment  [5,12,45].  Sadananda 
and  Shahinian  [45]  observed  that  the  prolonged  aging  increases  the  tensile  ductility  of  the 
alloy  by  nearly  fifty  percent,  they  believed  that  this  increased  ductility  causes  the  crack 
tip  to  blunt  readily  and  indirectly  decreases  the  environmental  sensitivity  to  the  crack 
growth  [46],  Floreen  and  Kane  [5]  have  noticed  that  the  fracture  surface  of  overaging 
heat  treatment  specimen  seems  to  reduce  the  amount  of  intergranular  crack  growth 
compared  to  that  of  conventional  heat  treatment  specimen.  They  then  speculated  that 
homogenization  of  slip  and  consequently  minimization  of  localized  stress  concentrations 
at  the  grain  boundaries  through  changes  in  precipitate  morphology  at  the  grain  boundaries 
is  a  way  in  which  overaging  heat  treatment  could  increase  the  resistance  to  crack  growth. 
Smith  and  Michel  [12]  suggested  that  the  improvement  in  crack  growth  resistance 
produced  by  the  modified  heat  treatment  is  likely  to  be  a  result  of  the  alteration  of 
dislocation  interactions  with  y"  particles  and  grain  boundaries.  Comparison  of  the  slip 
line  appearance  for  the  conventional  and  modified  heat  treatments  at  427°C  in  vacuum 
revealed  that  the  modified  heat  treatment  produced  a  microstructurai  condition  conducive 
to  homogeneous  defotmation  which  has  been  found  by  Fournier  and  Pineau  [15]  to  result 
in  significant  decrease  in  crack  growth  rate. 

In  sum,  one  can  conclude  that  the  types  of  microstructure  that  would  be  most 
susceptible  to  environmental  embrittlement  would  be  those  are  fine  grained,  have  very 
few  grain  boundary  strengthening  agents,  and/or  exhibit  pronounced  planar  deformation 
characteristics  at  elevated  temperatures.  Incidently,  it  has  been  reported  by  several 
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investigators  [5,34,36]  that  crack  growth  rate  in  vacuum  environment  was  almost 
independent  of  microstructure. 

Effects  of  Environment 

As  reviewed  previously,  the  high  temperature  fatigue  crack  growth  in  Alloy  718 
is  influenced  by  thermally-activated  and  time-dependent  processes.  These  processes  are 
generally  ascribed  to  the  phenomena  including  creep  and/or  oxidation  degradation 
processes.  The  relative  importance  of  these  two  processes  depends  on  a  number  of 
factors  such  as  the  aggressiveness  of  the  environment,  the  strengthening  characteristics 
of  the  material,  the  operating  temperature,  the  duration  of  exposure  time  and  the  applied 
load,  level.  Nickel  base  superalloys,  in  particular  Alloy  718,  are  a  class  of  high  strength 
alloys  that  have  been  designed  as  highly  creep  resistant  materials.  Brooks  and  Bridges 
[47],  for  example,  have  shown  that  Alloy  718  is  microstructurally  stable  after  10,000 
hours  exposure  at  temperatures  up  to  600°C.  Nicholas  ct  al  [23]  and  Pineau  [48]  have 
estimated  the  Ricdel-Ricc  characteristic  time  for  transition  from  small  scale  yielding  to 
extensive  creep  in  Alloy  718  at  650°C  to  be  in  the  range  of  20  hours  to  3  years.  This 
time  period  is  much  larger  than  the  possible  cyclic  periods  and  hold  time  durations  in 
practical  applications  of  this  material  in  gas  turbine  engines,  for  example,  but  the  material 
may  be  insufficient  for  applications  in  the  power  industry  where  lifetimes  arc  measured 
in  tens  of  years. 

When  examining  the  relationship  between  load  and  load  line  displacement  during 
creep-fatigue  tests  in  Alloy  718  at  tiSO^C,  Pedron  and  Pineau  [33]  found  that  the  specimen 
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behavior  is  predominantly  elastic,  thus  they  suggested  that  the  creep-fatigue  crack  growth 
behavior  of  Alloy  718  at  650°C  can  be  described  in  terms  of  linear  elastic  fracture 
mechanics,  in  other  words,  the  stress  intensity  factor  K  is  the  controlling  parameter  in  this 
case.  Sadananda  and  Shahinian  [45]  also  concluded  similarly  that  when  there  is  a 
significant  contribution  from  environmental  interactions,  creep  effects  are  minimized 
because  of  the  rapid  growth  of  cracks  and  under  these  conditions  the  stress  intensity 
factor  K  can  adequately  characterize  the  crack  growth  rates.  In  addition,  fractographic 
observations  of  Alloy  718  [5,6,25]  showed  that  there  was  no  evidence  of  any  creep 
cavitation  along  grain  boundaries  in  cycling  tests.  All  these  observations  demonstrated 
that  the  time-dependent  effects  in  Alloy  718  are  principally  due  to  the  environmental 
degradation  process. 

Besides  other  environments,  air  environment  has  received  more  attentions  and  is 
commonly  considered  as  an  aggressive  environment  for  the  resistance  of  high  temperature 
crack  growth  [5,21,26,30,33,34,38,49,50].  In  general,  the  air  environmental  effects  on 
crack  growth  behavior  are  studied  by  comparing  the  crack  growth  rates  in  both  air  and 
vacuum,  and  under  the  certain  loading  conditions  where  the  environmental  effects  can  be 
shown  up,  it  has  been  found  that  crack  growth  rates  in  air  are.  several  orders  higher  than 
their  counterpan  in  vacuum  [12,26,49]. 

One  of  the  major  issue  of  studying  the  influences  of  air  environment  on  fatigue 
life  or  fatigue  crack  growth  behavior  is  to  identify  which  constituents  in  air  actually  cause 
the  reduction  in  the  fatigue  life  or  the  resistance  of  the  fatigue  crack  growth  rate 
[26,51,52],  After  studying  the  fatigue  properties  of  pure  lead  in  air  at  various  reduced 
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pressures  and  in  separate  gas  atmospheres,  Snowden  [51,53,54]  found  that  oxygen  is  the 
element  in  air  environment  responsible  for  the  fatigue  life  reduction,  while  environments 
of  nitrogen,  hydrogen  or  water  vapor  had  much  less  effect.  Same  conclusion  has  been 
reached  by  Sadananda  and  Shahinian  [26],  Moreover,  in  static  load  tests  carried  out  on 
Alloy  718  at  650°C,  in  which  the  environment  was  alternated  from  air  to  vacuum  to  air 
[26],  the  history  of  environment  was  found  to  be  of  no  importance  and  the  environmental 
effect  was  found  to  be  limited  to  the  crack  tip  region.  This  localized  oxidation  around 
crack-tip  region  was  also  observed  by  Andrieu  et  al  [35,55]. 

It  has  been  realized  that  mechanical  loading  conditions  can  considerably  influence 
the  environmental  degradation  process  [56-68],  Coffin  [56]  is  probably  the  earliest  to 
recognize  the  effect  of  mechanical  loading  on  oxidation.  His  work  documented  the 
observation  of  highly  localized  regions  of  surface  oxidation  in  certain  high-temperature 
alloys  when  subjected  to  cyclic  plastic  strain  at  elevated  temperatures.  The  explanation 
of  the  observed  effects  was  based  on  localized  and  reversed  grain  boundary  deformation, 
leading  to  repeated  rupture  of  the  protective  oxide  film,  and  accelerated  oxidation  in  the 
region  of  deformation.  In  order  to  determine  if  oxygen  had  been  transported  into  the 
material  during  fatigue  crack  propagation,  Swanson  and  Marcus  [63]  have  analyzed  the 
fracture  surfaces  of  specimens  by  comparing  fatigue  crack  growth  rate  in  vacuum  with 
that  in  isotope  oxygen- 18  environment.  Their  results  indicated  that  fatigue  process 
enhances  the  oxygen  penetration.  Rcmy  and  associates  [66,67]  studied  fatigue  oxidation 
interaction  in  a  cobalt-base  superalloy  MAR-M509  and  a  nickel-base  superalloy  IN  100 
by  quantitative  metallography,  they  found  that  low  cycle  fatigue  at  high  temperature  can 
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strongly  enhance  matrix  oxidation  kinetics  in  these  superalloys,  and  the  enhancement  of 
matrix  oxidation  is  by  cyclic  stress  for  MAR-M509  but  is  by  cyclic  straining  for  IN10G. 
The  work  done  by  Ward  et  al  [62]  in  a  low  carbon  steel  and  by  Bueklow  and  Skelton 
[58]  in  a  Cr-Mo-V  steel  has  also  shown  that  cyclic  stress  has  a  strong  influence  on 
oxidation  kinetics.  Balsone  and  Nicholas  [68]  examined  the  effect  of  stress  on  the 
magnitude  of  oxidation  in  Rene '80,  their  results  demonstrated  that  oxidation  is 
significantly  greater  under  an  applied  stress  than  in  the  unstressed  condition,  and  the  depth 
of  oxidation  penetration  increases  as  the  stress  level  increases. 

So  far,  several  parameters  which  influences  the  fatigue  crack  growth  rate  in  Alloy 
718  have  been  reviewed,  and  the  factors  that  can  alter  the  fatigue  crack  growth  rate  in  this 
material  arc  shown  to  be  complex,  many  of  these  variables  may  interact  syncrgistically 
to  produce  a  variety  of  fatigue  crack  growth  behaviors.  The  above  review,  however,  can 
lead  to  the  following  conclusions: 

1)  The  high  temperature  fatigue  crack  growth  behavior  is  related  to  the  cyclic 
loading  frequency.  In  this  regard,  three  crack  growth  stages  can  be 
identified:  cyclic-dependent  stage  (at  high  frequency),  cyclic-  plus  time- 
dependent  stage  (around  transitional  frequency)  and  time-dependent  stage 
(at  low  frequency). 

2)  The  time-dependent  crack  growth  behavior  is  primarily  due  to  the 
environmental  influence. 

3)  Oxygen  is  the  element  responsible  for  the  oxidation  degradation  process. 


4)  Mechanical  loading  conditions  can  strongly  influence  the  oxidation 
process. 

It  should  be  noticed  that  the  time-dependent  crack  growth  stage  is  the  most  critical 
in  high  temperature  applications  due  to  its  highly  accelerated  crack  growth  rate,  it  is, 
however,  the  least  studied  or  understood  in  all  of  the  high  temperature  crack  growth 
stages.  This  time-dependent  crack  growth  stage  is  the  focus  of  the  present  study  with  the 
objective  of  understanding  the  crack  tip  controlling  mechanism  in  order  to  develop  the 
ability  to  predict  its  fatigue  crack  growth  performance  under  different  loading  conditions. 
Thus,  following  this  introductory  chapter,  the  subject  of  Chapter  2  is  to  analyze  the  nature 
of  the  time-dependent  crack  growth  process  by  studying  the  depth  of  intergranular  oxygen 
diffusion  in  the  crack  tip  region.  Chapter  3  focuses  on  the  intergranular  crack-tip 
oxidation  mechanism  and  outline  experimental  studies  carried  out  to  verify  the  proposed 
mechanism.  A  quantitative  crack  growth  model  based  on  the  knowledge  of  grain, 
boundary  oxidation  mechanism  and  localized  stress  and  strain  analyses  will  be  detailed 
in  Chapter  4.  Chapter  5  describes  the  numerical  simulation  of  this  model  and  examines 
its  predictive  capability  by  comparing  the  results  generated  from  the  model  with 
corresponding  experimental  data.  Chapter  6  summarizes  the  results  and  presents  the 
conclusions  of  this  research  program. 


CHAPTER  2 


ENVIRONMENT-DEPENDENT  FATIGUE  CRACK  GROWTH  STAGE 


In  this  chapter  attentions  will  be  made  to  establish  the  fact  that  the  environment- 
dependent  crack  growth  stage  in  Alloy  7 1 8  is  fully  controlled  by  the  depth  of  oxygen 
penetration  at  the  crack  tip  along  an  intergranular  fracture  path.  This  notion  will  be 
derived  mathematically  .and  supported  by  experimental  results  as  well. 

2.1  High-temperature  fatigue  Crack  Growth  Stages  in  Alloy  718 

As  discussed  in  the  previous  chapter,  attempts  to  predict  elevated  temperature  fatigue 
crack  growth  behavior  in  Alloy  718  have  been  considered  as  the  combined  effects  of  the 
cyclic  and  the  oxidation  damage  in  the  crack  tip  region.  In  this  regard,  the  damage 
process  ranges  from  purely  cyclic  dependent  at  high  frequency  levels  to  fully 
environment-dependent  at  very  low  frequency  levels.  The  relationship  between  loading 
frequency  and  environmental  effects  on  the  acceleration  of  the  fatigue  crack  growth  rate 
in  Alloy  718  can  be  explained  in  terms  of  the  intergranular  oxygen  diffusion  process  in 
the  crack  tip  region.  One  of  the  governing  factors  of  this  process  is  the  grain  boundary 
diffusivity  of  oxygen  —  a  subject  that  has  not  been  extensively  studied  particularly  in 
the  case  of  Alloy  718.  It  is  recognized,  as  discussed  in  last  chapter,  that  intergranular 
oxygen  diffusion  depends  on  the  stress  and  strain  states  along  affected  grain  boundaries, 
thus  the  diffusivity  of  a  stressed  grain  boundary,  Dg,  is  proposed  to  be: 


D  =  Dexp 


-Sl 

RT 


(2.1) 


where  D  is  a  diffusi vity  constant,  R  is  the  gas  constant,  T  is  the  temperature  in  Kelvin, 
and  Qg’  is  the  effective  activation  energy  of  the  grain  boundary  diffusion  which  could 
be  expressed  as: 


Qi  =  Q,-f(Wp)  <i2> 

where  Qg  is  the  activation  energy  of  the  grain  boundary  diffusion  in  the  stress-frees  state 
and  Wp  is  the  inelastic  strain  energy  density.  Through  this  definition  the  influence  of 
loading  frequency  and  the  associated  deformation  mode  on  the  magnitude  of  Dg  and  f(Wp) 
also  consequently  on  the  crack  growth  response  can  be  interpreted  qualitatively  as 
follows.  High  frequency  loading,  which  is  generally  characterized  by  high  slip  density 
and  a  homogeneous  form  of  deformation,  would  result  in  both  strain  accommodation  us 
well  as  stress  relief  along  affected  grain  boundaries  in  the  crack  tip  region.  Hence 
variations  in  f(Wp),  and  consequently  D0,  tend  to  be  minimal.  In  addition,  the  increase 
in  slip  density,  generally,  leads  to  an  increase  in  the  lateral  matrix  diffusion  across  the 
affected  grain  boundaries.  These  two  combined  effects  would  result  in  limited  or  no 
acceleration  of  intergranular  oxygen  diffusion  rate.  In  this  situation,  the  influence  of 
crack  tip  oxidation  is  minimal  and  crack  tip  damage  becomes  generally  dominated  by 
cycle-dependent  effects,  giving  rise  to  transgranular  fracture  mode.  On  the  other  hand, 
low  frequency  loading  accompanied  by  low  slip  density  would  promote  grain  boundary 
stress  concentration  resulting  in  an  increase  in  the  magnitude  of  both  f(Wp)  and  Dg.  This 
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is  magnified,  particularly,  if  stress  relieve  by  grain  boundary  sliding  is  not  permitted,  as 
in  the  case  of  the  highly  creep  resistant  Alloy  718.  Furthermore,  the  decrease  in  slip 
density  would  limit  the  grain  boundary  lateral  matrix  diffusion  process.  Here,  the 
expected  increase  in  the  grain  boundary  diffusivity  and  associated  increase  in  depth  of 
grain  boundary  oxidation  results  in  an  increase  in  the  crack  tip  damage  due  to 
environment  effects  giving  rise  to  intergranular  fracture  process. 

The  relationship  between  slip  line  homogeneity  and  fracture  mode  characteristics, 
as  described  above,  has  been  investigated  in  the  work  of  Ghonem  et  al  [7]  using  compact 
tension  specimens  made  of  Alloy  718  with  large  grain  size  (50-120  pm)  tested  at  650°C 
in  air.  In  this  work,  fatigue  fracture  mode  as  well  as  slip  line  density  were  compared  for 
loading  frequency  of  30  Hz  versus  that  of  0.05  Hz  which  includes  300  seconds  hold  time 
at  maximum  load  level.  Slip  line  traces  at  and  below  the  fracture  surface  of  the  test 
specimens  were  obtained  using  the  decoration  technique  described  in  [69],  These  traces 
for  different  AK  levels  in  both  test  conditions  are  shown  in  Figs.  2.1  and  2.2.  In  Fig,  2.1, 
which  corresponds  to  a  fully  transgranular  fracture  mode,  the  degree  of  homogeneous 
deformation  is  evident  by  the  high  slip  line  density  and  the  confinement  of  the  reversed 
plastic  zone  to  a  narrow  band  near  the  fracture  surface.  This  could  be  compared  with 
Fig.  2.2(b)  which  corresponds  to  a  fully  intergranular  fracture  and  displays  a  lower  slip 
line  density  as  well  as  larger  plastic  zone  size.  One  should  also  observe  that  in  the  three 
intergranular  locations  in  Fig.  2,2,  the  degree  of  slip  homogeneity,  measured  as  the 
inverse  of  the  slip  line  Interspacing,  increases  as  AK  increases. 

On  the  basis  of  above  argument  the  crack  growth  response  of  Alloy  718  with  respect 
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to  loading  frequency,  f,  has  been  divided,  following  the  work  of  Pineau  [48]  and  Nicholas 
and  co-workers  [6,7],  into  three  distinctive  types  as  shown  in  Figs.  2.3  and  2.4,  The  first 
type  is  associated  with  high  frequency  loading  in  which  the  deformation  mode  is 
governed  by  a  high  degree  of  slip  homogeneity.  Cracking  proceeds  primarily  in  the 
matrix  material  (in  contrast  to  the  grain  boundary),  resulting  in  a  predominantly 
transgranular  fracture  mode.  The  value  of  the  frequency,  f,  required  to  produce  this  type 
of  environment-independent  behavior,  decreases  as  the  magnitude  of  AK  increases.  This 
type  of  crack  growth  behavior  is  generally,  predicted  by  the  use  of  a  Paris-type  equation. 

As  loading  rate  decreases,  the  degree  of  slip  line  homogeneity  in  the  crack  tip  zone 
is  lowered.,  resulting  in  a  relative  increase  of  the  intergranular  oxygen  diffusion.  Under 
this  condition,  the  cruck  lip  damage  becomes  a  combination  of  oxidation  and  cycle- 
dependent  components.  In  this  type  of  response,  the  total  crack  tip  damage  has  been 
described  by  many  authors  [20,66,70-73]  using  various  models  which  are  generally 
reduced  to  a  form  of  the  linear  summation  rule  and  can  be  expressed  as: 


r  da  ^ 
dN 


/total 
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/time- dependent 


where,  in  the  absence  of  sustained  loading  effects  (hold  times),  the  time-dependent  term 
represents  the  contribution  of  oxidation  to  crack  tip  damage  and  could  be  written  as  a 
time  integral  in  the  form: 


^  da  ^ 


dt  (2.4) 

where  t„  represents  the  time  period  of  the  cycle  during  which  the  oxidation  process  is  an 
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active  damaging  component.  As  shown  in  Figs,  2.3  and  2.4,  this  oxidation  enhanced 
process  is  characterized  by  a  mixed  trans granular/intergranular  fracture  mode.  The  degree 
of  contribution  of  each  of  the  cycle-  and  time-dependent  terms  in  the  above  equation 
depends  on  both  the  frequency  and  AK  values.  For  the  same  frequency,  as  AK  increases 
the  conunbution  of  the  cycle-dependent  damage  also  increases,  since  increasing  AK  leads 
to  nn  increase  in  the  degree  of  slip  line  homogeneity,  see  Fig.  2.1.  The  increase  in  the 
cycle-dependent  damage  is  measured  by  the  increase  in  percentage  of  the  transgranular 
features  along  the  fracture  surface.  For  the  same  AK  value,  however,  the  influence  of  the 
time-dependent  damage  increases  as  frequency  decreases. 

The  third  type  of  response  mentioned  above  occurs  for  loading  frequencies  below 
a  transitional  level,  f0,  where  the  crack  tip  damage  becomes  mainly  an  environment- 
dependent  process  in  which  crack  growth  is  largely  intergranular.  The  value  of  the 
transitional  frequency,  for  a  particular  AK,  in  Alloy  718  was  found  to  depend  on  both 
temperature  and  microstructure  [7,17,48].  As  shown  in  Fig.  2.4,  for  a  grain  size  around 
100  (im  at  650°C  and  AK=40  MPaVm,  the  value  of  f0  is  on  the  order  of  0.1  Hz.  The 
influence  of  air  environment  on  fatigue  crack  growth  behavior,  where  the  test  frequency 
is  in  the  environment-dependent  regime  (f  <  fc),  is  depicted  in  Fig.  2.5  by  comparing  the 
fatigue  crack  growth  rates  tested  in  air  and  vacuum  (10'8  torr)  conditions.  It  clearly  shows 
that  the  fatigue  crack  growth  rate  in  vacuum  is  about  5  to  10  times  slower  than  that  in 
air.  Furthermore,  the  fracture  surfaces  in  both  cases  are  shown  in  Figs.  2.6(a)  and  2.6(b) 
where  they  exhibit  fully  intergranular  fracture  in  the  case  of  air  testing  and  fully 
transgranular  fracture  in  the  case  of  vacuum  testing. 


28 


TEMPERATURE,  R 
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Fig.  2.3  Schematic  of  high  temperature  fatigue  crack  growth  stages. 
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Ig.  2.4  Effect  of  frequency  on  fatigue  crack  growth  rate  of  Alloy  718  at  R  =  0.1  for  different  temperatures,  grain  sizes 
(S=20_50pm,  L=  150pm)  and  AK  levels. 
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Fig.  2,5  Comparison  of  fatigue  crack  growth  rate  in  air  and  vacuum  (frequency 
0.05  Hz,  temperature  650°  C). 
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(a) 


(b) 


Fig,  2.6  Scanning  electron  micrographs  of  fracture  surfaces  of  Alloy  718  tested 
with  0.05  Hz  at  650°C  in  (a)  air  and  (b)  vacuum. 
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2.2  Environment-Dependent  Fatigue  Crack  Growth  Stage 
2.2.1  Intergranular  Depth  of  Oxygen  Diffusion 

If  one  makes  the  assumption  that  in  the  environment-dependent  stage  the  crack 
growth  increment  per  cycle  is  equal  to  the  intergranular  depth  of  oxygen  diffusion,  X, 
occurring  during  the  cycle  effective  oxidation  time,  tM,  then  the  fatigue  crack  growth  rate 
can  be  expressed  as: 


=  _[X(AK,t)dt  (2.5) 

where  X  is  the  oxygen  intergranular  diffusion  rate,  generally  it  varies  with  time.  The 
recognition  that  oxygen  diffusion  is  an  energy  activated  process,  it  could  then  be  treated 
as  a  function  of  the  stress  intensity  factor  range,  AK,  acting  on  the  crack  tip  during  the 
cycle  effective  oxidation  time.  Testing  the  validity  of  this  assumption  (i.c.  the  ability  of 
the  above  equation  to  describe  the  environment-dependent  stage),  requires  the  knowledge 
of  both  tox  as  well  as  the  relationship  between  X  and  AK.  These  two  requirements  are  not 
readily  available  in  literature  and  therefore  an  attempt  will  be  made  here  to  determine 
these  requirements  experimentally, 

2.2.2  Cycle  Effective  Oxidation  Time 

The  cycle  effective  time,  t„t  is  defined  here  as  being  the  period  of  the  cycle  during 
which  the  oxidation  effects  take  place.  Several  authors  have  assigned  different  measures 
to  tox.  For  example,  Achter  et  al  [74]  in  their  study  of  the  effect  of  oxygen  partial 


pressure  on  crack  growth  rate  in  Type  316  stainless  steel  at  elevated  temperature, 
proposed  a  calculation  method  which  was  based  on  the  assumption  that  the  time  required 
for  adsorbing  a  gas  atom  monolayer  at  the  crack  rip  is  equal  to  half  of  the  tensile  part  of 
the  cycle  or  a  quarter  of  the  whole  cycle  period.  Wei  and  associates  [75,76]  in  their 
attempt  to  predict  environment  assisted  crack  growth  behavior  in  AISI  4340  steel, 
assumed  that  the  value  of  the  maximum  load  is  the  controlling  factor  for  crack  surface 
reaction  rate.  The  cycle  oxidation  time  in  their  work  was  assumed  to  be  equal  to  half  of 
the  loading  period  plus  half  of  the  unloading  period,  Nicholas  ct  al  [23,71],  in  their  work 
on  Alloy  718,  developed  a  model  to  predict  purely  time-dependent  fatigue  crack  growth 
behavior  by  integrating  the  sustained  load  growth  rate.  They  argued  that  the  loading  pan 
of  the  cycle  is  the  part  responsible  for  the  environment-assisted  effects.  Similar 
conclusions  were  made  by  Florcen  [52]  in  his  work  on  grain  boundary  diffusion  in  nickel- 
base  superalloys.  Other  investigators,  specially  Liu  an  McGowan  [77],  Rechet  et  al  [66], 
Antolovich  [78],  Romanoski  [79]  and  Saxena  [73]  employed,  for  different  materials, 
different  cycle  period  to  represent  t0„,  In  the  face  of  these  different  definitions  of  t,H>  an 
experimental  attempt  was  made  here  to  determine  the  effective  oxidation  time  for  Alloy 
718  at  temperature  level  of  650°C  and  for  loading  frequencies  less  than  fc.The  fatigue 
crack  growth  tests  were  conducted  under  constant  load  range,  AP,  and  R=0. 1  where  R  is 
the  load  ratio.  These  tests  involved  triangular  wave  forms  with  different  cycle  durations, 
all  satisfying  the  condition  that  f  <  fe.  These  durations  include  27.5  seconds  (25s-2,5-;), 
50  seconds  (25s-25s),  100  seconds  (25s-75),  anti  100  seconds  (90s- 10s);  labelled  here  as 
cases  A,  B,  C  and  D,  -wSpccrively.  Fracture  surfaces  corresponding  to  all  tests  exhibited 
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fully  intergranular  fracture  features  (confirming  that  f  <  f0),  as  seen  in  Fig.  2.7.  Results 
in  the  form  of  fatigue  crack  growth  rate,  da/dN  versus  AK  are  plotted  in  Fig.  2.8.  They 
show  the  fatigue  crack  growth  rate  for  case  C  to  be  identical  to  that  of  case  D  indicating 
that  for  the  same  frequency,  varying  the  ratio  of  loading  and  unloading  portions  of  the 
cycle  will  not  influence  the  fatigue  crack  growth  behavior.  This  result  is  contrary  to 
results  obtained,  for  example,  by  Coffin  [80]  in  his  work  on  fatigue  crack  initiation  in 
OFHC  copper  at  400°C  and  on  304  stainless  steel  at  650°C.  His  results  showed  that  in 
an  asymmetric  loading  cycle,  slow-fast  loading  is  more  damaging  than  fast-slow  loading. 
Similar  results  have  been  obtained  by  Ghonem  et  al  for  crack  growth  behavior  in  various 
ductile  materials  [81],  However,  their  results  have  been  interpreted  as  time-dependent 
crack  growth  caused  by  grain  boundary  cavitation.  In  this,  cracking  mechanism,  cavity 
growth  and  consequently  crack  growth  is  aided  by  slow  rate  loading  while  unloading 
results  in  cavity  healing  with  no  contribution  to  the  crack  growth  process.  This  cavity 
growth  related  mechanism  is  not  operative  in  the  highly  creep  resistant  Alloy  718.  In  this 
alloy,  for  time-dependent  environmental  effects  which  are  governed  by  temperature  and 
load  levels,  both  segments  of  the  loading  cycle  should  exert  the  same  damage  effects  at 
the  crack  tip.  Furthermore,  cases  of  A,  B  and  C,  which  are  of  different  frequencies  but 
identical,  loading  times,  did  not  result  in  similar  crack  growth  behavior  so  that  an  increase 
in  the  total  cycle  duration  yields,  as  expected,  an  increase  in  the  crack  growth  rate.  The 
conclusion,  based  on  these  observations,  is  that  for  loading  frequencies  lower  than  the 
transitional  frequency,  the  cycle  effective  oxidation  time  is  equal  to  the  total  time  of  the 
loading  cycle. 


Fig.  2,7a  Scanning  electron  micrographs  of  fracture  surfaces  of  Alloy  718  tested  at 
0.05  Hz  and  650°C  for  (a)  25s-2.5s  (b)  25s-25s. 
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Fig.  2.7b  Scanning  electron  micrographs  of  fracture  surface  of  Alloy  718  tested  at 
0.05  Hz  and  650°C  for  90s-  10s. 
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Fig,  2.8  da/dN  versus  AK  for  different  frequencies  with  different  loading  and 
unloading  rates. 
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In  comparing  these  results  with  those  of  Ashbaugh  [18],  it  should  be  noted  that  the 
variation  in  loading  frequency  between  case  D  and  the  other  three  is  less  than  a  factor  of 
4  whereas  Ashbaugh  investigated  the  effect  using  frequencies  which  differed  by  either  one 
or  two  orders  of  magnitude.  Further,  the  crack  growth  rates  differ  by  much  less  than  an 
c  da  of  magnitude  at  low  AK  values  in  the  results  of  the  present  work,  whereas 
Ashbaugh  reported  that  the  growth  rates  were  "comparable."  Although  the  loading 
frequency  in  the  asymmetric  wave  produced  growth  rates  in  Ashbaugh’s  study  which  were 
nearly  identical  to  those  from  symmetric  wave  forms  of  the  same  frequency  at  high  values 
of  AK  in  excess  of  30  MPaVm,  differences  appeared  at  lower  AK.  In  particular,  the 
slow-fast  (50s-5s)  wave  had  a  lower  fatigue  crack  growth  rate  than  the  symmetric  0.01 
Hz  (50s-50s)  while  the  fast-slow  (5s-50s)  had  a  considerably  higher  fatigue  crack  growth 
rate  than  the  (5s-5s)  symmetric  wave.  This  demonstrates  that  total  cycle  time  as  well  as 
loading  frequency  can  influence  growth  rate,  depending  on  the  particular  conditions. 

These  results  indicate  that  as  long  as  the  inhomogeneity  in  plastic  deformation  near 
the  crack  tip  has  been  produced  by  slow  strain  rates,  the  fatigue  crack  growth  rate  is 
dependent  on  frequency.  In  general,  however,  crack  growth  rates  in  nickel-base 
superalloys  due  to  combinations  of  hold  times  and  cyclic  loading  involve  a  complex 
interaction  which  can  be  attributed  to  a  variety  of  mechanisms  including  fatigue,  creep 
or  stress  relaxation,  environmental  degradation,  and  crack  tip  blunting  or  crack  branching. 
Further,  the  AK  value  can  influence  the  results  because  of  the  differing  contributions  of 
cycle  dependent  and  time  dependent  mechanisms. 
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2.2.3.  Relationship  between  X  and  Stress  Intensity  Factor  Range  AK 
2.2.3.1  Concept 

Work  of  Nichoi...  and  Weerasooriya  [23],  has  demonstrated  that  subjecting  Alloy 
718  to  elevated  temperature  fatigue  testing  with  a  loading  cycle  having  a  frequency  f  > 
f0,  and  an  imposed  hold  time  at  the  minimum  load  level  for  periods  up  ;  ■  000  seconds 
did  not  result  in  any  measurable  acceleration  in  the  crack  growth  rate  when  compared  to 
the  crack  growth  rate  due  to  the  base  cycle  without  the  hold,  time  periods.  By  iowering 
the  loading  frequency  to  levels  below  the  transitional  frequency  of  the  alloy,  an 
accelerated  intergranular  crack  growth  was  detected  when  addition  of  hold  times  at 
minimum  load  level  were  imposed  for  periods  as  small  as  30  seconds,  see  Diboine  et  al 
[25].  This  obsetved  increase  in  the  crack  growth  rate  is  interpreted  in  the  work  of 
Ghonem  et  al  [17],  as  being  a  result  of  further  intergranular  oxidation  taking  place  at  the 
crack  tip  during  the  hold  time  period,  provided  that  the  minimum  load  level  does  not 
contribute  to  the  mechanical  driving  force  of  the  crack  tip.  The  crack  growth  rate  i”"  ler 
this  type  of  loading  calculated  using  a  damage  summation  form  similar  to  that  of  eq.  (2.3) 
can  be: 


where  the  first  term  of  the  right  hand  side  of  the  above  equation  represents  the 
contribution  to  the  crack  growth  rate  due  to  the  reversed  part  of  the  cycle  and  the  second 
term  is  the  contribution  due  to  the  hold  time  period.  The  above  equation  can  be 
expressed  in  terms  of  (da/dt)  as: 
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'da  ' 

^  ycyo+li 


:|  A.^dt  J  A^dt 


dt 


(2.7) 


The  term  (da/dt)  is  viewed  as  being  equal  to  the  intergranular  oxygen  diffusion  rate  X 
which,  as  mentioned  before,  is  assumed  to  be  a  function  of  the  stress  intensity  factor 
range,  AK.  Therefore,  eq.  (2.7)  could  be  rewritten  as: 


t  A,  \ 


dt 


fii) 

UnJ.,. 

The  intergranular  depth  of  oxygen  diffusion,  X,  could  then  be  written  as: 


da 

V  _/cye+h 


(2.8) 
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The  differential  form  of  this  eauation  is  expressed  as: 


(2.9) 
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Where  At,  in  this  case,  is  the  hold  time  duration  V  Therefore 
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(2.11) 


The  determination  of  X  could  thus  be  achieved  through  the  knowledge  of  the  hold  rime 
period  th  and  the  crack  growth  rates,  (da/dN)syc+b  and  (da/dN)cyo  which  correspond  to  the 
total  loading  cycle  and  the  reversed  pan  of  the  cycle;  respectively.  In  order  to  solve  this 
equation  an  experimental  program  was  carried  out  to  provide  these  data  for  different  of 
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K  and  ^  values.  This  experimental  program  as  well  as  its  results  and  analysis  will  be 
described  in  the  following  section. 

2. 2.3.2  Experiments,  Results  and  Analysis 

In  this  program,  a  set  of  crack  growth  experiments  were  performed  on  compact 
tension  specimens  made  of  the  Ailoy  718  previously  described  at  temperature  level  of 
650°C  and  a  stress  ratio  of  0.1.  This  set  of  experiments  consisted  of  constant  AK  tests 
with  and  without  hold  time  periods  imposed  at  the  minimum  load  level.  Three  different 
values  of  AK  were  selected;  30,  40,  and  50  MPaVm.  For  each  level  of  AK  four  different 
hold  time  periods  were  investigated,  these  are  50,  100,  680  and  3600  seconds.  The  load 
level  during  the  hold  time  period  was  determined  such  that  the  hold  time  period  would 
produce,  at  the  longest  expected  crack  length,  a  stress  intensity  factor  with  a  value  lower 
than  that  of  the  fatigue  threshold  stress  intensity  factor  for  this  material  which  was 
estimated  to  be  12  MPaVm.  Under  this  condition,  an  observed  increase  in  the  crack 
length  would  be  purely  correlated  with  the  oxygen  influence  at  the  crack  tip.  Results  of 
these  tests  in  the  form  of  the  ratio  (da/dN)e},c+it  /(dn/dN)eyc  versus  th  for  different  values  of 
AK  are  shown  in  Fig.  2.9.  They  indicate  that  the  influence  of  the  hold  time  at  minimum 
load  level  <5n  the  crack  growth  rate  is  measurable  for  hold  time  periods  as  small  as  50 
seconds.  This  observation  confirms  similar  findings  reported  by  Diboine  and  Pineau  [25]. 
The  influence  of  the  hold  time  reaches  a  saturated  level  at  ^  values  that  increases  as  AK 
decreases.  These  results  in  the  form  of  (da/dN)cyo+h  /(da/dN)cyo  when  substituted  in 
eq.(2.11)  yield  values  of  X  calculated  as  a  function  of  th  and  AK.  This  type  of 
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load  level  on  (da/dN^^/tda/dN)—  for  different  Kh  levels. 


relationship  h.  ilbmted  in  Fig.  2.10.  It  shows  that,  for  all  values  of  AK,  X  decreases 
as  tb  increase  /.  I •;*  relationship  between  Xandth  could  be  interpreted  on  the  basis  of 
the  two-stage  oxidation  mechanism,  see  [82].  At  the  onset  of  the  oxidation  process,  the 
oxygen  diffusion  rate  reaches  its  peak  since  no  barrier  to  diffusion  exits.  As  the  oxidation 
time  increases  thus  permitting  the  formation  of  the  dense  chromia  (Cr203)  layer,  the 
oxygen  penetration  rate  decreases.  When  the  formation  of  Cr203  is  completed  after  a 
certain  transition  oxidation  time  which  is  microstructure  dependent  [50],  X  would 
approach  zero.  The  relationship  between  X  and  both  th  and  AK,  as  illustrated  in  Fig.  2.10, 
is  fitted  into  the  following  mathematical  form 


X  =  Gj(t)  AKm 


(2.12) 


where  m  is  a  constant  with  the  value  of  3.144,  and  G,(t)  assumes  the  following 
polynomial  form: 


Gj(t)  =  At'** e  T 


(2.13) 


A,  at  and  a2  are  coefficients  having  the  values  of  1.787x10'*,  0.632  and  28.67, 
respectively.  One  should  observe  that  in  the  above  equation  when  time  t  approaches 
infinity,  both  G^t),  and  consequently  X,  approaches  zero. 

The  validity  of  eq.  (2.12)  can  not  be  tested  through  direct  measurements  of  X,  it 
could,  however,  be  verified  through  the  use  of  X  to  derive  a  fatigue  crack  growth  rate 
expression,  which  can  be  tested  by  comparing  its  predictive  results  to  those 
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Fig.  2.10  Oxygen  intergranular  diffusion  rate,  X,  vs  hold  time,  th,  imposed  at 
minimum  load  for  different  AK  levels. 
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experimentally  generated.  This  is  achieved  as  follows.  One  can  write  the  crack  growth 
rate  equation  as: 


(2.14) 


where  t,x  is  the  effective  oxidation  time  during  u  loading  cycle.  Substituting  eq.  (2.12) 
into  eq.  (2.14)  and  considering  the  oxidation  time  of  a  loading  cycle  to  be,  as  previously 
discussed,  equal  to  the  total  cycle  time,  the  fatigue  crack  growth  rate  can  then  be  derived 


as: 


where 


X  =  ii  nG3(f)AKffl 
dN  2 


(2.15) 


G2(f)  »Af‘le'‘,f 

f  is  the  loading  frequency  while  A,  at  and  %  are  constants  having  the  same  numerical 
values  as  those  in  eq.  (2.12).  Eq.  (2.15)  could  now  be  tested  by  comparing  its  predictive 
results  with  those  experimentally  obtained  for  the  same  loading  frequency.  The  degree 
of  matching  between  these  two  sets  of  results  can  then  determine  the  validity  of  eq.(2.12). 
Eq.  (2.15)  was  applied  for  three  different  frequencies,  i.e.  0.01,  0.02  and  0.05  Hz  all  of 
which  are  below  f„  for  the  alloy  718  used  in  this  study.  Results  of  this  application  are 
shown  in  Fig.  2.11.  It  is  observed  that  a  reasonable  agreement  exists  in  the  full  range  of 
AK  between  the  values  and  trends  of  the  three  experimentally  obtained  data  sets  and  those 
theoretically  predicted  through  the  use  of  eq.  (2.15).  This  result  is  taken  in  support  of 
eq.(2.12)  as  being  a  valid  expression  for  estimating  the  oxygen  diffusion  rate  in  Alloy  718 
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Fig.  2.11  cla/dN  vs  AK  for  both  experimental  and  theoretically  predicted  results 


♦ 


* 


« 


o 


♦ 


♦ 


and  also  as  an  indirect  support  to  the  notion  that  the  growth  process  governed  by 
eq.  (2.15)  is  a  fully  environment-dependent  process. 

An  important  feature  of  eq.  (2.12)  is  in  its  use  to  directly  calculate  the  oxygen 
diffusivity  of  grain  boundaries,  Dg,  a  parameter  required  for  the  physical  understanding 
and  the  quantitative  modelling  of  the  crack  tip  oxidation  mechanism  in  Alloy  718,  see 
[83],  This  could  be  achieved  as  follows:  a  parabolic  diffusion  law  when  applied  to 
estimate  the  intergranular  depth  of  oxygen  diffusion,  X,  during  a  time  interval  t  could  be 
written  as: 

•  x.a/57  <2-16> 

where  a  is  a  geometrical  constant,  Dg  is  the  oxygen  diffusivity  of  grain  boundaries.  The 
intergranular  depth  of  oxygen  diffusion,  X,  during  the  time  interval  corresponding  to  one 
loading  cycle,  could  thus  be  calculated  as: 


X  a  a 

N 


(2.17) 


where  f  is  th<'  loading  frequency.  Another  identification  of  X  is  that  described  in 
eq.  (2.15).  By  equating  these  two  X-expressions,  i.c..  eqs.  (2.15)  and  (2  17),  an  explicit 
relationship  linking  the  oxygen  diffusivity  of  grain  boundaries,  Dg,  with  both  AK  and 
frequency  f  cun  be  derived.  This  relationship  is  obtained  here  as 


Dg  =  G3(f)  AK2tn  ;  AK>0  (2.18) 

where 
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This  relationship  between  D#,  AK  and  f  as  illustrated  in  Fig.  2.12  shows  that,  for  the  same 
AK  range,  the  diffusivity  Dg  is  inversely  proportional  to  the  loading  frequency.  This  cun 
be  explained,  as  qualitatively  discussed  in  eq.  (2.1),  by  considering  the  influence  of 
frequency  on  the  local  stress  field  near  the  grain  boundary  regions.  For  example,  the 
decrease  in  the  loading  frequency  is  expected  to  result  in  an  inhomogeneous  form  of 
plastic  deformation  which  consequently  leads  to  the  increase  of  stress  concentration  across 
the  affected  grain  boundary  paths  near  the  crack  tip.  This  would  then  result  in  the 
decrease  of  the  effective  dlffusional  activation  energy,  Qg',  of  these  affected  boundaries. 
The  decrease  in  Qg  is  expressed  in  cq.  (2.1)  through  the  introduction  of  an  inelastic  strain 
energy  function  f(Wp)  which  could  now  be  shown,  through  the  comparison  between 
eqs.  (2.1)  and  (2.18),  to  be  a  AK-  and  f-dependent  function.  Furthermore,  Fig.  2,12  shows 
that,  for  the  same  loading  frequency,  the  diffusivity  Dg  increases  as  AK  increases.  This, 
again,  can  be  explained  in  terms  of  the  direct  influence  of  the  crack  tip  stress  field  on  the 
diffusional  characteristic  of  the  grain  boundary  path.  Now,  if  one  make  the  assumptions 
that  the  ratio  of  matrix  diffusivity,  Dm,  to  grain  boundary  diffusivity,  Dg,  is  maintained 
at  lO4  [84]  and,  furthermore,  Dm  is  a  frequency-independent  parameter,  one  could  thus 
establish  profiles  of  Dm  for  different  AK  values.  These  profiles  are  shown  as  dashed  lines 
in  Fig.  2.1.2,  The  point  of  intersection  between  two  corresponding  Dm  and  Dg  lines  would 
then  identify  the  transitional  frequency,  f0,  for  a  particular  AK  level.  Below  f0,  Dg 
dominates  giving  rise  to  intergranular  crack  growth  mode,  while  above  f0,  Dm  dominates 
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resulting  in  transgranular  crack  growth  path.  In  this  case,  while  f0  decreases  as  AK 
increases,  the  variation  of  fg  is,  as  described  by  V/eerasooriya  et  al  [6,7],  to  the  narrow 
range  between  0.1  to  0.2  Hz. 


CHAPTER  3 


INTERGRANULAR  CRACK-TIP  OXIDATION  MECHANISM 


In  the  previous  chapter  it  is  established  that  the  environment-dependent  crack 
growth  stage  in  Alloy  718  is  fully  linked  to  the  crack  tip  oxygen  diffusion  process.  In 
this  chapter,  the  crack  tip  oxidation  process  will  be  detailed  and  will  be  then  used  in 
chapter  4  in  a  quantitative  approach  to  predict  the  crack  growth  behavior  in  this  alloy. 

3.1  Crack  Tip  Oxidation  Mechanisms 

As  mentioned  before,  extensive  evidence  exists  in  literature  showing  that  the  time 
dependency  of  high  temperature  crack  initiation  and  propagation  is  the  result  of  the 
aggressive  effect  of  environment.  The  early  work  of  Smith  et  al  [85]  has  clearly 
demonstrated  that  environmental  degradation  effect  is  a  result  of  oxygen  penetration  at 
the  crack  dp.  Efforts  have  been  made  in  the  last  two  decade  to  detail  the  crack  tip 
oxidauon  mechanism  and  provide  correlations  between  these  mechanisms  and  both 
material  and  loading  parameters.  When  not  considering  simple  adsorption  of  oxygen  at. 
the  crack  dp,  oxidation  mechanisms  could  be  identified,  in  general,  in  terms  of  oxygen 
short  and  long  range  diffusion  processes.  In  the  short  range  diffusion  process,  oxygen 
forms  an  oxide  layer  at  the  crack  tip  with  a  depth  that  depends  on  many  operating  and 
muterials  parameters.  The  forme, don  of  this  layer,  under  die  restricted  concave  crack  tip 
geometry,  results  in  high  stresses  that  could  easily  be  transmitted  to  the  substrate.  The 


important  aspect  of  this  oxidation  mechanism,  however,  is  the  possible  formation  of 
wedge-shaped  oxide  intrusions  along  the  crack  front.  The  rupture  of  these  wedges  at 
grain  boundary  intersections,  could  result  in  an  accelerated,  intergranular  crack- growth 
rate,  see  Bricknell  and  Woodford  [S6]  McMahon  and  Coffin  [87],  Coffin  [88]  and 
McMahon  [89].  In  the  long  range  diffusion  process,  oxygen  penetrates  the  crack-tip 
material  along  rapid  diffusion  paths,  such  as,  slip  planes  and  grain  boundaries.  The 
internal  oxidation  process  taking  place  along  these  path  could  occur  in  the  form  of 
internal  oxide  sites,  cavity  formations,  and/or  solute  segregation.  As  pointed  out  by 
Woodford  and  Bricknell  [90],  if  is  also  possible  for  oxygen  to  take  part  in  chemical 
reactions  releasing  known  embrittlement  agents  onto  grain  boundaries.  Each,  or  all,  of 
these  may  be  operative  in  any  particular  alloy  under  a  given  set  of  conditions.  These 
processes,  in  particular,  along  grain  boundaries  result  in  the  inhibition  of  the  sliding  and 
migration  of  these  boundaries  and  thus  reduce  their  ability  to  relieve  'ocal  stresses  built 
up  during  deformation. 

While  it  is  recognized  that  the  oxidation  mechanisms  associated  with  the  short  and 
long  range  diffusion  are  not  completely  separated,  experimental  observations  indicate  that, 
for  Alloy  718,  at  an  intermediate  temperature  range  and  at  lower  cyclic  loading 
frequencies,  the  short-range  oxygen  diffusion  contributes  primarily  to  the  occurrence  of 
an  intergranular  fracture  mode.  Recognizing  that  oxygen  partial  pressure  plays  an 
important  role  in  the  development  of  this  oxidation  process,  Andrieu  [50]  has  suggested 
a  qualitative  model  which  depends  on  the  type,  and  sequence,  of  oxide  layers  formation 
at  the  crack  tip.  This  model  will  be  termed  as  the  two-stage  crack-tip  oxidation 
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mechanism  in  the  following  chapters  of  this  report. 

3.2  Two-stage  crack-tip  oxidation  mechanism 

The  damage  process  associated  with  short-range  oxjgen  diffusion  at  the  crack  tip 
has  been  studied  by  several  authors  [34,58,9 1J.  important  factors  to  be  considered  in  this 
process  are  those  related  to  temperature,  frequency,  alloy  chemistry  and  oxygen  partial 
pressure,  P02.  The  first  clear  demonstration  of  the  effect  of  P02  was  exhibited  in  the 
work  of  Smith  et  al  [85].  They  observed  that  reducing  the  oxygen  partial  pressure  during 
fatigue  crack  propagation  in  3 16. stainless  steel  at  both  500  and  800°C  led  to  a  lower 
crack-growth  rate.  Transition  to  higher  crack  growth  rates,  as  F02  increased,  was  also 
observed,  among  others,  by  Smith  and  Michel  [92]  in  their  work  on  a  Co-base  alloy  at 
427°C,  by  Stegmann  and  Shahinian  [93]  on  nickel  alloy  at  500°C  and  by  Smith  and 
Shahinian  [94]  on  Silver  at  20,  150  and  350°C.  This  transition  was  interpreted  by  Achter 
[95]  in  terms  of  the  impingement  rate  of  oxygen  molecules  on  successive  rows  of  freshly 
exposed  metal  atoms  at  the  crack  tip.  As  discussed  by  Ericsson  [96],  the  geometry  of  the 
impingement  concepts  has,  in  general,  been  described  in  vague  terms  facing  the  difficulty 
of  providing  an  explanation  for  the  severity  of  damage  that  a  monomyer  ei  oxide  could 
cause  at  the  crack  tip.  These  concepts,  furthermore,  cannot  support  the  experimental 
observations  made,  for  example,  by  Smith  and  Shahinain  [94].  They  reported  that,  the 
effect  of  oxygen  partial  pressure  on  fatigue  life  of  Silver  at  350°C  is  small  compared  to 
the  effect  of  that  at  20°C  and  150°C.  An  increase  in  the  initial  rate  of  oxygen  adsorption 
has,  however,  been  observed  with  increased  temperature  at  the  range  of  20°C  to  350°C. 


Another  interpretation  of  the  influence  of  P02  on  crack  growth,  is  related  to  the 
ability  of  oxygen  at  different  pressures  to  control  the  preferential  formation  of  certain 
oxides  that  may  shield  or  contribute  to  the  crack  tip  damage  through  the  passivation  or 
the  enhancement  of  the  oxidation  process,  respectively.  It  is  known  that  oxidation  of 
some  elements,  e.g.  chromium,  can  produce  protective  oxide  layers,  and  chromia  (Cr20-j) 
is  such  oxide.  While  oxides  of  less  noble  elements,  e.g.  the  iron  and  nickel  oxides  and 
their  spinei  type  oxides,  are  less  protective  [87].  If  chromium  diffusion  in  the  alloy  is  too 
low  to  sustain  the  supply  of  chromium  at  the  scale/metal  interface,  less  noble  elements 
of  the  alloy  will  be  oxidized  at  the  scale/metal  interface.  This  will  result  in  an  increased 
reaction  rate  [95] .  In  this  regard,  Andrieu  [50]  has  experimentally  examined  the 
influences  of  P02  on  the  formation  of  selective  oxides.  In  his  work,  small  discs,  with  a 
thickness  of  3  mm  and  diameter  of  20  mm,  were  machined  from  Alloy  718  with  an 
average  grain  size  of  150  |lm.  These  electro-  and  mechanically  polished  discs  were 
subjected  to  different  P02  levels  ranging  from  3  x  10‘‘  to  10'5  torr  for  periods  of  480 
seconds  at  650°C  using  a  high-vacuum  devise  equipped  with  a  mass  spectrometer.  An 
Auger  spectrometry  analysis  was  subsequently  used  to  identify  the  oxide  type  formed  on 
the  free  surfaces  of  these  discs  as  a  function  of  the  operating  P02  level.  Sputter  profiles 
of  the  formed  oxides  are  summarized  in  Fig.  3.1.  It  is  interesting  to  observe  that  oxygen 
partial  pressure  as  high  as  104  torn  promotes  selective  oxidation  of  chromium  and  that, 
under  the  pressure  from  10'1  to  10'2  torr,  the  first  oxides  to  appear  are  nickel  and  iron 
oxides.  The  relationship  between  time  required  to  form  the  chromium  oxide  as  a  function 
of  oxygen  partial  pressure  is  schematically  shown  in  Fig.  3.2;  at  atmospheric  pressure  this 
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Fig.  3.1 


Sputter-depth  profiles  of  the  surface  oxides  (oxidation  time  8  min, 
temperature  650°C). 
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Alloy  718  (free  surface) 


Fig.  3.2  Types  of  oxides  as  a  function  of  transition  time  tp  and  oxygen  partial  pressure. 


passivation  time  is  estimated  to  be  in  the  range  of  5-8  minutes.  Throughout  this  study 
it  was  difficult  to  distinguish  whether  the  surface  being  analysis  was  of  a  matrix  or  of  a 
gain-boundary  interface.  Transmission  Electron  Microscopy  work  was,  therefore,  carried 
out  on  a  foil  made  of  the  same  material  in  which  grain  boundary  was  identified.  This  foil 
was  subjected  to  a  temperature  of  650°C  for  240  seconds  in  room  environment.  Results 
of  this  TEM  work  are  shown  in  Fig,  3.3.  They  indicated  that,  under  atmospheric 
pressure,  microcrystallized  oxides  of  nickel  and  iron  developed  rapidly  along  grain 
boundaries. 

Furthermore,  attempts  have  been  made  to  identify  the  role  of  P02  in  determining 
the  cracking  mode.  This  was  achieved  by  carrying  out  a  series  of  slow  strain-rate,  tensile 
tests  on  oxisymmetrically  precracked  fatigue  round  bars  made  of  the  same  alloy.  These 
tests  were  performed  at  650°C  in  an  ultra  high  vacuum  apparatus  detailed  in  Ref.  [50]. 
Using  a  mapping  technique  employed  on  the  resulted  fracture  surfaces,  examples  of  which 
are  shown  in  Fig.  3.4,  the  percentage  of  transgranular  fracture  area,  with  respect  to  the 
total  fracture  surface  area,  was  measured.  These  percentages  in  terms  of  oxygen  partial 
pressure  are  reported  in  Fig.  3.5.  They  indicate  that  the  intergranular  fracture  ratio 
increases  at  the  P02  increases  and  under  atmospheric  pressure,  60%  of  the  fracture  surface 
of  the  test  specimens  was  found  to  be  intergranular.  Also,  a  transition  exists  between  10'2 
and  1  torr,  below  which  the  failure  mode  is  fully  transgranular.  'Phis  range  of  transitional 
pressure  is  consistent  with  observations  made  on  different  alloys  [92-94]. 

These  results,  when  linked  to  those  in  Fig.  3.1  showing  the  influence  of  P02  in 
determining  the  type  and  sequence  of  selective  oxides,  could  provide  a  qualitative  view 


Fig.  3.4a  Fracture  surfaces  corresponding  to  different  oxygen  partial  pressures  (a) 
lO-4  torr  and  (b)  1  torr. 
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Fig.  3.4b 


Fracture  surface  corresponding  to  the  oxygen  partial  pressure  of  4  torn 
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Fig.  3.5  Ratio  of  transgranular  fracture  versus  oxygen  partial  pressure. 


of  the  crack-tip  oxidation  along  a  preferred  grain  boundary  path.  This  view  suggests  that 
oxidation  occurs  in  two  stages,  schematically  depicted  in  Fig.  3,6.  The  first  stage, 
influenced  by  atmospheric  level  of  oxygen  partial  pressure,  would  result  in  the  formation 
of  FeO  and  NiO  and  their  spinels.  For  these  oxides  which  are  permeable  to  oxygen,  the 
gas  diffusion  through  the  pores  plays  an  important  role  in  determining  the  reaction  rate 
as  a  function  of  time.  The  build-up  of  these  oxide  products  would  result  in,  at  the 
oxide/metal  interface,  lowering  the  oxygen  partial  pressure.  This,  consistent  with  the 
Ellingham/Richardson  diagram,  would  provide  the  reaction  kinetics  required  for  the 
formation  of  a  dense  Cr203  oxide' sub-layer.  This  dense  layer  is  considered  protective  in 
the  sense  that  it  limits  oxygen  diffusion  to  the  grain  boundary  material  and  thus  decreases 
the  metal/oxygen  reaction  rate  [25].  The  complete  formation  of  this  layer  along  the  crack 
front  is  regarded  as  un  attainment  of  oxide  depth  saturation.  The  time  required  to  reach 
this  stage  is  identified  as  the  transition  time,  tp,  which  is  a  function  of  temperature,  P02 
and  localized  stress  ad  strain  fields.  When  this  time  is  reached,  an  assumption  is  mude 
that  the  condition  for  oxygen  passivation  is  met  and  no  further  oxygen  penetration  can 
take  place  along  the  effected  grain  boundary  puths.  One  should  observe  that  the 
passivation  effect  is  not  assumed  here  to  be  a  discrete  process,  it  is  a  continuous  gradual 
decrease  in  the  oxide  depth  formation  rate  as  the  ratio  of  the  Cr203  to  the  total  oxide 
weight  increases. 

3.3  Experimental  Verification 

This  two  stage  crack-tip  oxidation  mechanism  was  verified  by  Andrieu  et  al  [82], 
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Fig.  3.6  Suggested  mechanism  of  grain  boundary  oxidation. 


64 


In  their  work,  two  sets  of  high  temperature  fatigue  crack  growth  tests  were  carried  out. 
These  experiments  utilized  compact  tension  type  specimens  made  of  Alloy  718  with  grain 
size  of  about  20-50(im.  The  objective  of  the  first  set  of  experiments  was  to  investigate 
the  influence  of  crack-tip  oxidation  during  a  time  duration  equal  to  transition  time,  tp,  on 
the  subsequent  fatigue  crack  growth  rate.  This  was  accomplished  by  comparing  the  crack 
growth  rate  of  continuous  cycling  (labeled  case  A)  with  that  of  cycling  imposed  a  hold 
time  at  minimum  load  level  (labeled  case  B).  The  duration  of  this  hold  time  was  selected 
to  be  600  seconds  which  is  in  the  order  of  the  transition  time,  tp,  of  this  alloy  [25,49], 
In  order  to  insure  intergranular  fracture  mode,  the  frequency  used  in  the  study  was  0.05 
Hz  which  is  below  the  transitional  frequency  of  this  material.  Results  of  these  tests  are 
shown  in  Fig.  3.7.  They  show  a  distinct  increase  in  the  crack  growth  rate  as  a  result  of 
Imposing  a  hold  time  at  minimum  load  during  which  crack-tip  oxidation  is  assumed  to 
be  the  only  active  process.  These  results  could  be  explained  on  the  basis  that  the  hold 
time  has  permitted  unintemruped  metal/oxygen  reaction  leading  to  the  formation  of  the 
protective  chromia  layer.  Upon  applying  a  loading  cycle  with  a  frequency  that  induces 
intergranular  cracking,  the  depth  of  ruptured  oxides  along  the  intergranular  fracture  path 
would  represent  an  additional  jump  of  the  crack  tip.  Furthermore,  the  type  of  oxides 
formed  on  the  fracture  surface  of  the  air  tested  specimen  was  identified  by  analytical 
transmission  electron  microscopy  on  a  thin  foil  of  a  fracture  surface  edge.  Fig,  3.8 
indicates  the  existence  of  two  oxide  layers  -  the  inner  layer  (A)  at  the  metal/oxide 
interface  is  identified  as  a  Cr203  while  the  outer  layer  (B)  has  been  identified  as  a  spinel 
oxide  Ni(Fe  Cr)204  in  the  same  order  as  suggested  in  the  proposed  model. 
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Fig.  3.8  TEM  observation  and  microanaiysis  on:  oxides  types  developed  on 
fracture  surface. 
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Ghonem  et  al  pointed  out  that  the  depth  of  oxide  layer  under  the  hold  time  tests 
condition  would  correspond  to  a  saturated  depth  since  the  hold  time  is  of  the  order  of  the 
transition  time.  However,  if  a  high  frequency  minor  cycle  is  imposed  on  the  hold  time 
period  in  case  B,  then  the  metal/oxygen  reaction  processes  will  be  disturbed  by 
continuously  rupturing  the  rapid  formation  of  nickel  and  iron  oxides.  This  would 
continuously  alter  the  oxygen  partial  pressure  on  the  metal/oxides  interface,  consequently, 
the  oxygen/metal  reaction  rate  will  keep  high  which,  in  turn,  prolongs  the  duration  of  the 
transition  time,  tp,  and  delays  the  formation  of  protective  chromia  layer.  As  a  result,  the 
oxide  depth  along  the  intergranular  path  will  be  deeper.  This  idea  was  examined 
experimentally  by  repeating  case  B  with  a  high  frequency  (1  Hz)  minor  cycle  imposed 
on  the  hold  time  duration  (illustrated  in  Fig.  3.9),  labeled  case  C,  and  the  corresponding 
results  in  terms  of  da/dN  versus  AK  are  shown  also  in  Fig.  3.7.  These  results,  which 
clearly  indicate  an  increase  in  tiie  crack  growth  rate  when  compared  with  those  results  of 
tests  A  and  B,  are  seen  here  as  an  indirect  support  of  the  two  stage  crack- tip  oxidation 
concept. 

The  two  stage  crack-tip  oxidation  concept  was  then  applied  in  the  present  study 
to  develop  a  quantitative  model  to  predict  the  high  temperature  fatigue  crack  growth  rate 
for  Alloy  718,  this  will  be  the  subject  of  Chapter  4, 
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Fig.  3.9  Load  spectrum  including  a 
hold  time  at  minimum  loac 
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CHAPTER  4 


MATHEMATICAL  MODELLING  OF  INTERGRANULAR  FATIGUE 
CRACK  GROWTH  BEHAVIOR  AT  650°C 


In  this  chapter,  a  mathematical  model  employing  the  basic  concepts  of  the  two  - 
stage  crack-tip  oxidation  mechanism  will  be  developed.  As  a  First  step  towards  this 
objective,  an  attempt  will  be  made  in  the  following  section  to  review  the  existing  fatigue- 
oxidation  crack  growth  models.  ' 

4.1  Review  of  Existing  Fatigue-Oxidation  Crack  Growth  Models 

As  mentioned  in  the  previous  chapter,  one  of  the  objective  of  understanding  the 
crack  tip  oxidation  mechanism,  particular  in  the  high  temperature  alloys,  is  to  provide 
basis  for  the  development  of  quantitative  models  that  predict  crack  growth  for  these  alloys 
under  different  temperatures  and  loading  conditions.  As  a  result,  a  wide  spectrum  of 
prediction  models  have  been  developed  for  different  materials  based  on  different  views 
of  the  fatigue-oxidation  interaction  process.  One  of  the  earlier  models  to  quantitatively 
describe  the  high  temperature  fatigue  crack  growth  behavior  is  that  proposed  by  Haigh 
et  al  [98].  Their  model  can  be  written  as: 

ii=C(f)AKn  (4.1) 

dN 

where  n  is  a  constant  and  C(f)  is  a  frequency  dependent  term  which  reflects  the  time 


dependency  or  the  oxidation  effects.  This  Paris  type  model  was  applied  to  Cr-Mo-V  steel 
in  both  air  and  vacuum  at  550°C  but  was  not  extended  to  other  materials  or  different 
ranges  of  temperature. 

Liu  and  McGowan  [77]  made  the  assumption  that  within  the  environmental  effect 
region,  the  fatigue  crack  growth  rate  is  controlled  by  the  rate  of  the  transport  of 
detrimental  chemical  species.  They  correlated  the  fatigue  crack  growth  rate  with  the 
Arrhenius  relation  in  the  following  form: 
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where  f0  and  T0  arc  frequency  and  temperature  at  a  reference  state,  respectively,  B,  tn  and 
n  are  constants,  Q  is  the  activation  energy  and  R  is  the  gas  constant.  The  data  of  IN  100 
and  304  stainless  steel  supported  the  analysis  based  on  this  empirical  model,  the  data  of 
Waspalloy  was  not,  however,  sufficient  to  draw  any  conclusion. 

Several  other  models  [58,60,61,66]  have  built  their  argument  on  the  assumption 
that  the  high  temperature  fatigue  crack  growth  per  cycle  in  air  environment  can  be 
expressed  as  a  linear  summation  of  the  following  damage  components: 
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where  the  first  term  represents  the  pure  time-independent  mechanical  damage,  while  the 
second  term  represents  the  oxidation  influences  on  the  crack  growth  behavior.  In  Skelton 
and  Bucklow’s  work  on  Cr-Mo-V  steel  at  550°C  [58]  the  oxidation  term  was  written  as: 


(4.4) 
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where  a  is  a  constant,  As  is  the  total  strain  range  and  K  is  a  magnification  factor  which 
depends  on  the  strain  range,  Pilling-Bedworth  Ratio  (PBR)  and  current  crack  length.  If 
oxygen  penetration  becomes  more  significant,  Marshall  [61]  suggested  that  the  above 
equation  can  be  replaced  by 


'  da 


(4.5) 


where  D  is  the  fatigue-enhanced  diffusion  coefficient  and  t  is  the  cycle  period. 

To  account  for  oxidation  damage,  Strangeman  [60]  introduced  the  term  "metal 
recession  rate  by  oxidation  per  cycle"  MRR  which  only  be  calculated  from  experimental 
data.  Again,  based  on  the  linear  summation  damage  concept,  his  equation  is  written  as: 


ii,A(AK)B+MRR  (4.6) 

dN 

where  A  and  B  are  material  constants  and  AKe  is  the  strain  intensity  factor  range.  The 
dependency  of  the  model  on  the  term  AKt  limits  its  application  to  the  condition  of  very 
low-strain-intensity  factor  values  at  which  the  crack  tip  strains  are  insufficient  to  fracture 
the  oxide  scale,  a  condition  which  is  incompatible  with  the  strain  state  in  gas  turbine 
components,  for  example. 

Remy  [66,99]  also  applied  the  linear  summation  rule  to  model  the  fatigue- 
oxidation  interaction  process,  but  more  sophisticated  metallurgical  method  was  used  in 
his  work  to  evaluated  the  oxidation  damage.  He  represented  the  "pure"  fatigue  damage 
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based  on  the  work  of  Tomkin  [100],  while  the  oxidation  effects  was  derived  from  his 
metallurgical  observations  on  partly  fatigue  specimens,  i.e. 


f  da 
^dN 
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(4.7) 


where  fe  is  the  effective  volume  of  carbides  on  the  crack  path,  Bsn  is  a  constant,  At  is  the 
cycle  period,  gCAe,,)  is  a  function  of  plastic  strain  range,  Aep,  ct^  and  ae  are  the  oxidation 
constants  of  the  matrix  and  the  carbides  under  zero  stress,  respectively.  This  equation 
was  applied  successfully  to  model  the  behavior  of  MAR  M509  under  various  temperature 
and  loading  conditions. 

Based  on  the  concept  of  decohesive  strength  of  grain  boundary,  Ramanoski  et  al 
[79]  postulated  a  fatigue  crack  growth  model  which  depends  on  the  effective  stress  range 
Aaoff: 


ii  =  AAcS,  (4.8) 

dN 

where  A  and  m  arc  constants.  If  expressing  the  effective  stress  range  in  a  form  which 
incorporates  a  threshold  stress  for  microcrack  extension  aT,  the  above  equation  can  then 
be  written  as: 


—  =A(o  -<JT)  (4.9) 

dN  m“  T 

This  model  represents  a  balance  between  the  mechanical  driving  force  and  the  resistance 
of  the  material  to  this  force.  The  relationship  between  the  cohesive  strength  of  the  grain 
boundary  and  environmental  factors  was  not,  however,  considered.  In  addition,  while  the 
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oxidation  process  is  time  dependent,  neither  time  nor  test  frequency  is  apparent  in 
eq.(4.9). 

Attempts  also  have  been  made  to  model  hold  time  effects,  during  which  the 
influence  of  oxidation  and/or  creep  might  be  activated  [101-104],  Utah  [103]  developed 
the  MSE  model  which  is  based  on  the  following  equation: 


A  =  exp(B) 
dN 
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(4.10) 


where  AK  is  the  stress  intensity  factor  range,  AK,  and  AKe  represent  the  threshold  and  the 
maximum  values  of  AK,  respectively,  B,  P  and  D  are  constants  for  a  given  set  of 
experimental  conditions.  Annis  et  al  [104]  have  also  proposed  a  SINH  model  as  follows 


log 


^  da  ^ 
vdN, 


=  C.Sinh(C,(logAK  +  C3))  +  C, 


(4.11) 


where  Q,  Q,,  C,  and  C4  are  constants  which  can  be  expressed  as  functions  of  frequency, 
temperature,  hold  time  and  stress  ratio.  Predictions  of  both  MSE  model  and  SINH  model 
were  made  of  the  fatigue  crack  growth  rate  in  alloy  AF115  over  a  temperature  range  from 
538°C  to  760°C,  stress  ratios  from  0.1  to  0.9,  frequencies  from  0.0025  Hz  to  2.5  Hz,  and 
hold  times  up  to  5  minutes.  In  general,  these  two  models  could  represent  the  data  set, 
from  which  the  constants  of  the  model  are  derived,  reasonably  well.  Here  one  should 
realized  that  while  both  MSE  model  and  SINH  model  are,  in  principle,  regression  models 
and  are  not  based  on  any  physical  mechanisms,  they  are  useful  for  fitting  to  a  particular 
data  set. 

Based  on  the  linear  summation  rule,  Nicholas  et  al  [69]  have  developed  an 
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empirical  relationship  for  complex  wave  shapes  as  follows 

^  =  f(K,R)|ci+C2(T)v-J(T)  +  C3(T)|J  4 dt  +  p  J  4dt]j  (4.12) 

where  C,  is  a  constant,  vL  and  vu  are  frequencies  of  loading  portion  and  unloading 
portion,  respectively,  \  =  1/(2Vl)  and  tu  =  l/(2vu)+th,  where  xh  is  the  hold  time.  The 
first  term  of  the  model  represents  the  purely  cycle  dependent  behavior  and  is  obtainable 
from  low  temperature  and/or  high  frequency  test.  The  second  term  represents  the  pure 
cycle  dependency.  The  last  two  terms  represent  the  time  dependent  behavior,  the  integral 
over  the  loading  portion  of  the  cycle  xL  is  always  considered,  but  existence  of  the  integral 
over  the  unloading  and/or  sustained  load  portion  of  the  cycle  tu  is  dependent  on  p  which 
varies  from  zero  to  unity.  There  is  no  method  for  determining  the  functional  variation 
of  p  at  this  stage. 

From  the  above  review,  it  is  clear  that  these  models  could  generally  be  classified 
into  two  main  groups.  The  first  one  is  a  phenomenological  approach  obtained  by  a 
regression  analysis  procedure  and  carried  out  on  a  set  (or  sets)  of  experimental  data  or 
by  modifying  the  Paris-type  law  to  include  parameters  such  as  frequency  and  temperature 
to  account  for  the  time  dependency  nature  of  the  high-temperature  crack  growth  process 
[69,98,103,104].  These  models  are  generally  simple  in  formulation  and  contain  a  limited 
number  of  constants  that  are  not  difficult  to  determine.  However,  the  absence  of  physical 
interpretation  of  their  basic  equations  makes  it  difficult  to  extend  their  applications  to 
materials  or  loading  conditions  different  than  those  used  to  generate  the  model,  'rhe 
second  approach  is  based  on  physical  interpretation  of  metallurgical  and  mechanical 
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observations  of  a  specific  crack  growth  mechanism  [58,60,66,77,79,99].  These  models, 
while  complex,  provide  good  predictions  for  the  same  material  under  all  conditions  that 
generate  an  identical  cracking  mechanism.  Different  observations  and  formulations  will, 
however,  be  required  for  a  different  fracture  mechanism. 

A  general  conclusion  that  could  be  reached  after  reviewing  all  these  models  is  that 
not  one  oxidation-assisted  crack  growth  model  can  be  viewed  as  valid  for  all  materials 
under  different  loading  conditions.  The  validity  of  a  particular  model  must  be  assessed 
in  relation  to  a  specific  material,  temperature  and  cracking  mechanism.  According  to  this 
conclusion,  the  work  described  in  this  chapter  was  undertaken  to  quantitatively  model  the 
intergranular  oxidation  enhanced  fatigue  crack  growth  process  in  Alloy  718  at  650°C 
based  on  the  understanding  of  the  mechanism  of  intergranular  oxidation  and  the  related 
aspects  described  in  Chapters  2  and  3. 

4.2  Model  Concept 

The  concept  of  the  proposed  model  can  be  described  as  follows.  At  the  start  of 
the  oxidation  process,  which  is  taking  place  along  a  preferred  grain  boundary  path  in  front 
of  the  crack  tip,  a  layer  of  selective  oxides  will  replace  the  affected  grain  boundary 
material.  The  depth  of  the  oxide  layer  depends  on  the  degree  of  oxygen  diffusion  along 
the  grain  boundary.  The  survival  of  this  layer,  as  a  load-bearing  material,  is  determined 
from  the  balance  established  between  the  mechanical  driving  force  at  the  crack  tip  and 
the  resistance  of  the  layer  to  this  force.  This  notion  is  represented  in  Fig.  4-1  in  which 
the  grain  boundary  material  which  a  ductility,  ed0,  measured  in  a  high  temperature 
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Fig.  4,1  Schematic  of  the  proposed  model  concept. 
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vacuum  environment,  is  replaced  by  material  with  a  variable  degree  of  oxidation  and 
consequent  variable  amplitude  of  ductility,  ed.  The  distribution  of  ed  is  influenced  by  the 
fact  that  the  closer  the  grain  boundary  material  is  to  the  crack  tip,  the  higher  the  degree 
of  oxidation  and  the  higher  the  degree  in  reduction  of  ed0.  The  changes  in  ed  will  reach 
its  saturation  level  when  the  oxidation  time  equals  or  exceeds  the  transition  time,  tp, 
required  to  achieve  passivation,  i.e.,  time  required  to  complete  the  formation  of  Cr203 
layer.  Acting  in  the  same  zone  near  the  crack  tip  is  the  strain  distribution,  ^  produced 
by  the  externally  applied  load.  One  assumption  made  here  is  that  the  material  dealt  with 
here  will  have  clastic-viscoplastic  behavior  and  the  imposed  loading  conditions  will  only 
result  in  small-scale  yielding. 

The  magnitudes  of  these  two  competing  parameters,  Ed  and  will  determine  the 
length  of  crack  tip  advancement  along  the  affected  grain  boundary  path.  The  crack 
growth  rate  could  therefore  be  calculated  by  determining  the  value  of  r  per  cycle  where 
r  is  the  distance  measured  from  the  crack  tip  which  corresponds  to  position  P  in  Fig.  4.1, 
i.e.,  when  e m  S  ed.  The  present  model  thus  depends  on  determining  the  evolution  of  both 
the  and  ed  distribution  on  a  cycle-by-cycle  basis.  Steps  to  determine  these  distribution 
are  the  subject  of  the  following  section. 

4.3  Mathematical  Derivations 

At  high  temperature  and  in  an  air  environment,  the  oxidation  of  an  alloy  occurs, 
in  genera),  through  the  formation  of  multiple  oxides.  For  Alloy  718,  the  possible 
formation  of  major  oxides  and  their  corresponding  spinels  arc  Cr203,  FeO,  NiO,  FeCr204 


and  NiCr204.  At  the  beginning  of  the  oxidation  process  the  yield  of  Cr203,  which  is  the 
only  protective  layer  in  these  various  oxides,  is  small  and  thus  the  oxidation  effect,  in 
terms  of  the  depth  of  oxygen  penetration  and  the  degree  of  grain  boundary  embrittlement, 
is  relatively  large.  As  the  oxidation  time,  to*,  increases,  the  weight  of  Cr203  in  relation 
to  the  total  weight  of  oxides,  increases.  When  to*  equals  tp  the  formation  of  Cr203  is 
completed  and  the  condition  of  passivation  is  reached.  This  indicates  that  no  further 
oxidation  can  take  place  and  a  saturation  in  the  grain  boundary  ductility  level  of  the 
oxidized  zone  has  been  established.  One  could  thus  correlate  the  amount  of  Cr203  the 
grain  boundary  ductility  by  introducing  a  oxidation  parameter,  p,  which  is  defined  as  the 
ratio  between  the  weight  of  Cr203  and  the  total  weight  of  oxides,  i.e. 


Cr,0: 

m 


(4.13) 


For  Alloy  718  eq.  (4.13)  was  found  to  take  the  form  (see  Appendix  III): 


C1(Y)m& 

C2(Y)m0i~C3(Y)mCr 


(4.14) 


where  Clt  Q,,  and  C3  arc  functions  of  the  parameter  y,  m^  and  rriQ  are  the  amount  of 
oxygen  and  chromium  diffusion  along  the  grain  boundary  and  the  dislocation  network, 
respectively.  The  parameter  y  is  defined  as  mF,  /mN1,  where  mr<  is  the  amount  of  iron 
required  to  form  FeO,  in  the  grain  boundary,  and  mNl  is  the  amount  of  nickel  required  to 
form  NiO,  in  the  grain  boundary.  The  ratio  yean  be  experimentally  determined  by  using 
different  experimental  techniques. 
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From  the  result  obtained  in  Chapter  2  it  is  clear  that  the  effective  oxidation  time 


per  cycle  can  be  expressed  as: 


<fc-f> 

f|fc-fi 


(4.15) 


where  th  is  the  hold  time  component  of  tine  cycle,  f  is  the  loading  frequency,  fc  is  the 
tronsgranular/intergranular  transition  frequency  and  <•>  is  the  MaCauley  bracket.  If  the 
loading  frequency  is  above  the  transition  frequency,  the  crack  advances  along  a 
transgranular  path  and  no  oxidation  effect  will  be  considered,  in  this  case  <*>  =  0.  The 
term  fc  is  temperature  and  grain  size  dependent  [19,22]. 

A  link  could  now  be  established  between  the  oxidation  parameter  P  and  the  grain 
boundary  ductility.  In  the  early  stages  of  the  oxidation  process  where  P  is  small,  a  faster 
degree  of  reduction  in  edo  is  expected;  as  P  increases  the  reduction  in  ed  should  decrease. 
Upon  reaching  tp,  p  reaches  p*  and  assumes  a  saturated  value.  The  relationship  between 
P  and  ed,  can  be  written  as: 

cd=a1arctg(-a2p2  +  aJ)  +  at  (4.16) 

where  a,  to  \  are  material,  temperature  and  loading  condition  dependent. 

Turning  ones  attention  to  the  oxygen  diffusion  along  the  grain  boundary,  one  could 
view  it  as  a  short-circuit  diffusion  process.  This  process,  generally,  embraces  diffusion 
of  metal  or  oxygen  through  extended  defects  in  the  oxide  scale  rather  than  through  the 
bulk  crystal  lattice  [106].  Most  mathematical  treatments  represent  the  short-circuit  path 
as  an  isotropic  slab  of  material  with  a  width  of  26  in  which  diffusion  occurs  according 
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to  Fick’s  law.  In  this  treatment  the  concentration  variation  across  the  slab  is  neglected 
and  the  diffusivity  is  assumed  to  be  much  longer  than  the  diffusion  in  the  surrounding 
matrix  ['07-109].  In  following  this  approach  the  concentration  of  diffusion  species  along 
a  grain  boundary  path  as  a  function  of  position  (x,y)  and  time  t,  can  be  obtained  using 
Whipple’s  solution  [110]. 

Peterson  and  Martin  [108,111,112]  have  examined  Whipple’s  solution  and  pointed 
out  that  whilst  it  is  relatively  more  accurate  than  other  solutions,  such  as  that  of  Fisher 
[106],  it  suffers  from  certain  limitations.  For  example,  it  is  not  appropriate  for  "wide" 
grain  boundaries  since  the  concentration  across  the  width  way  may  not  be  uniform, 
Furthennore  it  can  only  be  applied  to  diffusion  problems  within  the  Harrison’s  type  B 
regime  [1 12],  In  this  paper,  the  concentration  across  a  grain  boundary  path  in  Alloy  718 
is  considered  constant  and  the  matrix  diffusion  distance  is  much  smaller  than  that  of  the 
grain  boundary;  Whipple’s  solution  will  then  be  assumed  appropriate  for  our  applications. 
The  amount  of  oxygen  diffused  along  the  affected  grain  boundary  within  one  loading 
cycle  can  thus  be  obtained  in  the  form 
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Dg  and  Dra  arc  diffusivities  of  oxygen  in  the  grain  boundary  and  the  matrix,  respectively 


whilst  r  is  the  distance  measured  from  the  crack  tip  and  mo2  is  the  oxygen  concentration 
in  an  air  environment  at  650°C. 


During  the  oxidation  process  an  enrichment  of  chromium  was  observed  along  the 
oxidized  grain  boundaries  [25,102],  This  chromium  depletion  from  the  interior  of  the 
grain  structure  can  be  evaluated  by  considering  that  chromium  is  transported  to  the  grain 
boundaries  along  edge  dislocation  lines.  By  assuming  an  effective  dislocation  pipe  with 
a  radius  d  and  a  effective  diffusion  distance  of  chromium,  d^,  similar  to  eq.  (4.17),  the 
amount  of  chromium  transported  through  the  effective  dislocation  pipe  during  one  loading 
cycle  can  be  calculated  as 
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where 


v©7-  D- 

Here  Da  and  Dm  are  the  diffusivities  of  chromium  in  the  effective  dislocation  pipe  and 
matrix,  respectively  while  mCr  is  the  concentration  of  chromium  in  the  matrix. 

An  important  step,  now,  is  to  determine  the  influence  of  the  stress  and  strain  states 
on  the  grain  boundary  diffusivity.  It  is  recognized  that  intergranular  oxygen  diffusion 
depends  on  the  stress  and  strain  states  along  affected  grain  boundaries  [22,56-68].  In 
Chapter  2,  for  a  stressed  grain  boundary  the  oxygen  diffusivity  of  grain  boundaries,  Dg, 
has  been  expressed  as  a  function  of  the  inelastic  strain  energy  density,  f(wp),  of  this 
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boundary.  This  is  written  as  an  Arrhenius  relationship  in  the  following  form: 


Dg«  Dexp 


=  Dexp 


/  , 

RT 


Qg-f(wp) 
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(4.19) 


where  Qg  is  the  activation  energy  of  grain  boundary  diffusion  in  the  stress-f^ee  state,  Qt’ 
is  the  effective  activation  energy  of  grain  boundary  diffusion  for  stressed  material,  D  is 
a  diffusion  constant,  R  is  the  gas  constant  and  T  is  the  temperature  in  Kelvin.  The 
function  f(Wp)  is  expressed  here  as  follows: 


f(Wp)  =  b  +  b2cxp 


(bjW.  +  bj  )2 


(4.20) 


where  bt  to  b3  are  coefficients  which  depend  on  temperature  and  material.  An  equation 
similar  to  eq.  (4.20)  can  be  written  for  the  dislocation  diffutivity  in  the  stressed  state  and 
D,j  can  be  completely  determined. 

The  parameter  Wpl  in  the  most  general  case,  for  loading  conditions  which  include 
both  cyclic  and  hold-time  components,  is  expressed  as: 


Wp  =  Wp  +  Wpc  (4.21) 

where  Wpf  and  Wpc  are  the  components  of  Wp  corresponding  to  the  contributions  of  the 
cyclic  process  and  the  hold-time  duration,  respectively.  The  cyclic  component  is 
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(4.22) 


w;=ix(o*m(odt 

where  am  and  are  the  equivalent  stress  and  strain  rate  near  the  crack  tip  region,  which 
are  caused  by  external  cyclic  loading.  The  hold-time  component,  Wpe,  on  the  other  hand, 
could  be  written  as  [)  13] 

wa_l-|i-2^Ka  (4.23) 

P“~E  T 

where  p.  is  Poisson's  ratio,  E  is  Young’s  modulus,  Km  is  the  stress-intensity  factor  during 
the  hold  time  duration. 

It  is  obvious  from  the  above  two  equations  that  knowledge  of  stress  and 
viscoplastic  strain  fields  in  the  crack-tip  region  is  required.  A  number  of  authors  [114- 
118]  have  worked  in  this  field  and  several  finite  element  codes  are  available  for 
evaluation  stress  and  strain  states  under  a  visco-plasdc  condition  [115-118].  Each  of  these 
codes  has  its  merits  and  limitations.  All  in  general,  however,  require  a  complex 
procedure  and  a  large  set  of  input  parameters.  In  this  paper,  we  will  attempt  to  use  the 
closed-form  solution  derived  by  Riedel  and  Rice  [119]  with  the  assumption  that  we  are 
dealing  with  a  small-scale  yielding  problem.  The  stress,  olj(  and  strain,  e,j,  near  the  crack 
tip  (0  =  0)  could  be  written  as: 

alj  =  AL'1!TTaiJ(e)|e.0  (4.24) 

and 
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(4.25) 
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Here  B  and  n  are  the  coefficients  and  power  of  Norton’s  law,  respectively,  while  Kj  is  the 
stress  intensity  factor,  Ig  is  a  numerical  parameter  which  depends  on  n,  0^(0),  andee(0) 
are  the  functions  of  0.  The  term  $(,(0)  is  the  deviator  of  $,,(0)  Substituting  eq.  (4.24) 


for  the  equivalent  stress  am,  i.  e. 
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and  noting  that  the  deviatoric  stress  tensor  c^’  can  be  written  as: 

°U  =  afij“rrCTkk^ij 


(4.26) 
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eq.  (4.26)  thus  becomes 
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In  the  near-tip  region  the  total  strain  is  composed  of  elastic  and  inelastic 
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components 


Cm  = 


em  +  em 


(4.29) 


and  since  eJ  »  £ m*  in  the  near-tip  region,  £„  -  ej.  The  inelastic  strain,  zj,  which 
actually  is  viscoplastic  strain,  could  consist,  in  the  most  general  case,  of  the  hold-time  and 
cyclic  components,  i.  e. 


1  f  c 
C»  +  Em 

The  material  behavior  is  described  here  by  Norton’s  law 


(4.30) 


=  (4.31) 

During  the  hold  time,  ej  ~  thus  by  combining  eqs.  (4.28)  and  (4.31)  then  integrating 
during  the  hold  time  duration,  one  obtains 
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and  K,  and  th  are  the  stress-intense  :y  factor  and  the  hold  time,  respectively.  The  cyclic 


strain  component,  ej,  could  also  be  calculated  as  follows 
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■bTTv 


eL  -  — f  (KI"7rt  "7r  +  2nt  Kt)dt  (4.33) 

n+1  Jt 

which  when  summed  with  eq.  (4.32)  gives  the  total  inelastic  strain  of  the  form 
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(4.34) 

Noting  that  e*  »  ej,  one  can  substitutes  eq.  (4.34)  into  Norton’s  law,  then  the  equivalent 
stress  am  is  obtained 
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Eqs.  (4.34)  and  (4.35)  represent  the  equivalent  strain  and  stress  near  the  crack  tip. 


The  purpose  of  model  is  to  generate  da/dN  versus  AK,  this  achieved  as  follows: 

(i)  for  a  desired  value  of  AP  and  an  arbitrary  distance  r  measured  from  the  crack  tip, 
AK  can  be  calculated; 

(ii)  the  corresponding  £„,  and  ara  are  calculated  using  eqs.  (4.21),  (4.34)  and  (4.35); 


(iii)  from  knowledge  of  Wp,  the  parameters  Dg  and  E)d  are  determined  using  eqs.  (4.19) 
to  eq.  (4.21); 

(iv)  values  of  m^  and  me,  can  be  obtained  using  eqs.  (4.17)  and  (4.18); 

(v)  using  an  experimentally  determined  value  of  y.  the  oxidation  parameter  [3  can  be 
obtained  from  eq.  (4,14)  and  the  corrcsp."  ling  value  of  ed,  at  the  distance  r,  is 
fully  determined. 

Using  an  iterative  procedure  would  yield  a  value  of  r  which  satisfies  the  failure 
criterion  £„  S  ed.  This  distance  r  is  taken  as  the  crack  length  increment  occurring  at  this 
specific  crack  length  and  this  particular  cycle.  From  this  procedure  the  relationship  of 
da/dN  versus  AK  can  be  generated. 

In  the  following  chapter,  a  numerical  simulation  routine  will  be  conducted  for 
predicting  the  fatigue  crack  growth  rates  under  the  two  different  loading  conditions.  These 
results  then  will  be  compared  with  the  data  generated  from  experiments  for  the  same 
loading  conditions  in  order  to  examine  the  validity  of  the  proposed  model. 


CHAPTER  5 


NUMERICAL  SIMULATION  OF  THE  PROPOSED  MODEL 


In  this  chapter,  the  proposed  model  is  numerically  applied  to  predict  fatigue  crack 
growth  rate  for  two  different  loading  conditions.  The  first  one  is  a  cyclic  load  profile 
including  a  hold  time  duration  for  300  seconds  imposed  at  the  minimum  load  level,  while 
the  second  is  a  continuous  cyclic  load  profile  without  hold  time  durations.  The  frequency 
of  the  reversed  part  of  the  cycle  and  the  stress  ratio  of  these  two  tests  are  0.05  Hz  and 
0.1,  respectively.  A  flow  chan  detailing  the  application  steps  is  shown  in  Fig.  5.1.  The 
computer  code  required  to  execute  this  chart  is  listed  in  Appendix  V.  The  predicted 
results  using  the  proposed  model  are  compared  with  those  generated  from  fatigue  cruck 
growth  tests  under  the  same  loading  conditions.  Details  of  the  data  generation  and 
comparison  are  given  in  the  following  sections. 

5.1  Cyclic  Loading  with  Hold  Time  at  Minimum  Load 

In  this  application,  a  hold  time  duration  of  300  seconds  is  imposed  at  the 
minimum  load  level.  The  cyclic  loading  frequency,  f,  is  selected  as  0.05  Hz  in  order  to 
satisfy  requirements  for  obtaining  an  intergranular  fracture  mode,  see  ref.  [24].  The  stress 
ratio,  R,  is  0.1.  For  this  type  of  loading  cycle,  the  stress-intensity  factor  can  be  expressed 


Fig.  5.1  Flow  chart  showing  the  calculation  procedure  involved  in  performing  the 
proposed  model. 
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where  is  the  stress-intensity  factor  during  the  hold  time  period;  is  equal  to  K^. 
In  this  load  profile,  the  load  level,  Pmlfl)  during  the  hold  time  period  is  determined  by  the 
condition  that,  at  the  longest  expected  crack  length  (a/W  =  0.7),  the  corresponding  stress- 
intensity  factor  would  not  exceed  the  threshold  value  of  the  test  material  which  is 
estimated  to  be  12  MPaVm.  In  this  way  it  is  assumed  that  no  "mechanical"  crack  growth 
could  be  produced  during  the  hold  time  period  and  any  crack  advance  during  this  period 
could  thus  be  entirely  attributed  to  the  effect  of  the  environment. 


STEP  1: 


The  first  step  in  the  numerical  procedure  required  to  determine  the  crack  growth 
rate  is  the  arbitrary  selection  of  a  distance  r  measured  from  the  crack  tip  to  the  point  at 
which  the  localized  equivalent  stress,  om,  and  the  localized  equivalent  strain,  em,  would 

be  calculated.  The  expression  for  am  and  e„  are: 

il  t  f  2n  „  t  n-1  Y.  , 

g=  ,  .  ..  V171  K1"^t~'n^  +  2nt^K,7CTfc.  f  0<t<4  (5.2) 

m  [(n+DBj  ^  v  f 
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-n  „ 1  2n 


t  n-1 


em=  ^r“^j[7(KIWt^+2nt^KI^KI)dt+F(^)r^thirn'  (53) 


The  parameters  F,  B  and  n  are  determined  as  8.34  X  10‘9,  4.5  X  10‘33  and  n  =  10, 
respectively  [25],  When  considering  that  =  RAK/(1-R),  the  above  two  equations 
could  be  simplified  into: 

crm  =  1995.52r_0'M1(  180.0  K/‘8l8t'0,909  +  K[°'8Ut0,091Kr)0,1  O^t^i  (S4) 


em  =  (3.310  +  3. 889t°'091) 


It  should  be  mentioned  here  that  in  the  present  analysis,  the  residual  strain,  due  to  the 
crack  growth  history,  is  neglected.  This  approximation  is  justified  on  the  basis  that  in  the 
small  scale  yielding  condition,  the  viscoplastic  zone  size  in  front  of  the  crack  tip  is  of  the 
same  order  of  magnitude  as  the  distance  of  the  crack  advance  per  cycle  [119]. 

STEP  2: 

Once  arm  and  em  are  determined,  the  inelastic  strain  energy  density,  Wp,  of  the 
grain  boundary  at  the  distance,  r,  can  be  evaluated  from  the  expression 


(5.6) 


wp-wpf+w; 


*  fT  a  m(t)  8  (t)  dt  4.  ™ 

Jo  “w  mW  r 


where  p,  is  Poisson’s  ratio  and  E  is  Young’s  modulus.  The  strain  rate  can  be 
calculated  using  Norton’s  power  law: 


.1  ^  n 

1(1  ••• 


Eq.  (5.6)  can  be  modified  to 


wo  =  _ lIciB-  ^ + B  fT  cC1  dt 

p  tcE  r  Jo 


Now,  taking  E  =  lt>5  GPa  [25]  and  |i  =  0.283  [120],  and  setting  Km  «  K^,  Wp  can  be 
reduced  to 

W  *  2.976 xlO‘6^l  +  2.882x1 0-4^!L 


STEP  3: 

Knowledge  of  Wp  permits  the  determination  of  the  diffusivities  of  both  the  grain 
boundary  path  and  the  effective  dislocation  pipe,  Dg  and  Dd,  respectively.  The  term  Dg 
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accounts  for  oxygen  diffusion  along  the  grain  boundary  while  Dd  accounts  for  the 
chromium  depletion  in  the  matrix  towards  the  grain  boundary.  These  parameters  are 
expressed  as 


Dg  =  D  exp 


(5.10) 


and 


where 


Dd  =  D,exp 


(5.11) 


f(Wp)  =  27.6313  -  25.8863  exp 


0.5202 


(5.6114W  +2.1711)2 


and 


g(Wp)  =  32.2362  -  32.0087  exp 


0.3712 

(35.7822  Wp  -  37.5720)2 


\ 


J 


Other  values  are  D  =  1.477  x  10'6  and  Dj  =  1.450  x  10‘7.  The  activation  energy  of  the 
grain  boundary,  Qg,  and  that  of  the  effective  dislocation  pipe,  Qd,  are  115000  J/mol  [121] 
and  128143  J/mol  [108],  respectively.  The  gas  constant,  R,  is  8.314  J/mol*K. 

STEP  4: 


During  one  loading  cycle,  the  amoum  of  the  oxygen  diffused  along  the  affected 
grain  boundary,  m,2,  and  those  of  chromium  transported  through  the  effective  dislocation 


pipe,  n^,,  are  determined  by  Whipple’s  solution  as  follows 


A,  - 1 
Aj  -u 


(5.12) 

where 

T|1  =  — — -  ,  . — Q  Kt  =  (1  -  At)4t  .  A,  = 

and 


fflft  (^Cr '  ^  )  =  ^Cr 


where 


K2  =  (l-A2)52 


(5.13) 


Tire  grain  boundary  width,  25,  and  the  diameter  of  the  effective  dislocation  pipe,  d,  are 
taken  to  be  comparable  in  size  and  equal  approximately  to  5  A  [107,108].  Furthermore, 
the  diffusivities  of  both  oxygen  and  chromium  in  the  matrix  are  assumed  to  be  the  same 
[108],  thus  here  Dm  and  D0  are  all  taken  as  1.2  x  10'14  cmV  [102,107].  The  term  m,2  and 
are  concentrations  of  oxygen  and  chromium  in  air  environment  and  in  the  matrix  at 
650°C,  respectively.  By  applying  the  ideal  gas  law,  m^  can  be  estimated  as  88.73  g/m\ 
As  for  the  value  of  m^,  since  the  chromium  content  in  Alloy  718  is  about  18.3%(wt),  and 
noting  that  the  density  of  this  alloy  is  8.19  g/m5  [120],  i%  could  thus  be  estimated  about 


1.4499  x  106  g/m3.  The  term  t‘  is  the  effective  oxidation  time,  which  is  defined  in  the 
previous  chapter  as  follows: 


t*  =  th  + 


<f0-f> 


(5.14) 


where  fe  is  the  transgranular/intergranular  transition  frequency,  which  is  estimated  to  be 
in  the  order  of  10’1  Hz,  here  f0  is  selected  as  0.2  Hz.  The  value  of  t*  in  the  present 
application  is  320  seconds.  Within  a  first  approximation,  the  effective  diffusion  distance 
of  chromium,  dcr,  is  estimated,  by  using  the  parabolic  rate  law,  to  be  about  10  pm. 
STEP  5: 

The  oxidation  parameter  (3,  detailed  in  previous  section,  is  expressed  as  a  function 
of  m^  and  mCr.  It  can  be  written  as  follows 


P" 


^•imCr 

C2m0j-C3mcr 


(5.15) 


where 


Ct  =  1.462(0.951  +y) 
C2  =  4.442 +  4.49  8  y 
C3  =  0.659 +0.6107 


and  y  is  the  rate  of  the  amount  of  iron,  mFe,  required  to  form  FeO  and  the  amount  of 
nickel,  mN1,  required  to  form  NiO.  In  the  present  applications  y  is  taken  as  1.125  [102]. 

STEP  6: 

The  grain  boundary  ductility,  Ed,  is  related  to  the  degree  of  oxidation  which  is 


expressed  by  the  oxidation  parameter  p.  The  relationship  between  ed  and  p,  as  discussed 
previously,  follows  a  sigmoidal  curve  which  is  described  by 


ed  =  a,  arctg(~a2p2  +  a3)  +  a4  (5,16) 

where  a„  a2,  a2,  and  a4  are  parameters  which  could  vary  with  material  characteristics, 
service  temperature,  loading  condition  and  other  parameters  that  may  contribute  to  the 
grain  boundary  embrittlement. 

Taking  ed0  =  5.37%  [102]  and  using  high  temperature  creep  crack  test  data  [22], 
the  relationships  between  ed  and  p  at  different  AK  levels  are  obtained  by  applying  a  trial 
and  error  curve  fitting  technique.  The  normalized  ed  /ed0  versus  P  for  AK  values  in  the 
range  of  18  to  40  MPaVm  is  presented  in  Fig.  5.2.  The  determination  of  Ed  for  any  other 
AK  level  within  this  range  can  then  be  performed  by  interpolation. 

STEP  7: 

Steps  1-6  are  carried  out  for  different  r  values.  The  value  of  r  which  results  in 
satisfying  the  fracture  criterion,  namely,  e*  =  Ed,  is  taken  here  to  be  equal  to  the  crack 
advance  per  cycle,  i.e., 

da 

'  dN  <s-17> 

The  pattern  of  the  distributions  of  both  the  mechanical  strain,  £„,,  and  the  gTain 
boundary  ductility,  ed,  calculated  as  functions  of  the  distance  r  are  shown  in  Fig.  5.3.  It 
can  be  seen  that  as  r  increase,  ^  decreases  and  Ed  increases  approaching  the  value  Ed0. 
The  intersection  of  the  E^  curve  and  the  Ed  curve  represents  the  position  of  the  crack  tip 
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Normalized  grain  boundary  ductility  versus  parameter  p  for  the  case  of  cyclic  loading  with  300  seconds 
hold  time  at  minimum  load  level. 
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at  the  end  of  the  corresponding  cycle. 

The  calculated  results  for  Wp,  mo2  /m^,  rtic,  /!%,  [3,  ^  and  da/dN  are  listed  in 
Table  5.1,  and  illustrated  in  Fig.  5.4  which  shows  a  comparison  between  the  experimental 
data  and  the  calculated  results  in  terms  of  da/dN  versus  AK.  The  predicted  results  agree 
very  well  with  the  experimental  data. 

5.2  Cyclic  Loading  without  Hold  Time 

In  this  application,  all  the  test  material  parameters  are  identical  to  those  used  in 
the  previous  application  except  for  the  hold  time,  th,  the  effective  oxidation  time,  t‘,  and 
the  effective  diffusion  distance  of  chromium,  d^,  hence,  the  calculation  steps  arc  the  same 
as  in  the  previous  application.  Here,  the  hold  time,  tj,,  and  consequently  Kn,,  is  taken  as 
zero  since  no  hold  time  duration  is  included  in  the  loading  profile.  The  effective 
oxidadon  time  t*  during  one  cycle  is  20  seconds  and,  the  effective  diffusion  distance  of 
chromium,  dc,,  estimated  by  applying  the  parabolic  rate  law,  is  2.6  pm.  The  relationships 
between  Ej  and  (3  at  different  AK  levels  are  obtained,  by,  again,  udlizing  a  trial  and  error 
curve  fitting  technique.  Fig.  5.5  describes  the  normalized  ejEio  versus  (3  for  AK  levels 
in  the  range  of  20  to  40  MPaVm.  The  pattern  of  the  distributions  of  both  and  Ed  in 
front  of  the  crack  tip,  for  this  case,  is  shown  in  Fig.  5,6,  while  the  calculated  results  for 
Wp,  m0l  /mo2,  m^,  [3,  e^,  and  da/dN  are  listed  in  Table  5.2.  The  comparison  between 

the  experimental  data  and  the  calculated  results  in  terms  of  da/dN  versus  AK  is  also 
presented  in  Fig.  5.4.  Again  in  this  case,  the  predicted  results  agree  well  with  the 
experimental  data. 
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Fig.  5.4  Comparison  between  experimental  and  numerical  simulated  results  in 
terms  of  da/dN  versus  AK. 
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versus  parameter  p  ior  the  case  of  cyclic  loading  without  hold  time. 


AK 

(MPaVm) 

W» 

(J/mol) 

nWmoj 

mcMct 

xlO° 

CQ. 

em 

<%) 

da/dN 

(mm/cyc) 

tested 

da/dN 

(mm/cyc) 

predicted 

18.52 

2.7375 

0.9883 

7.548 

0.5049 

1.687 

6,0xt  04 

5.84X104 

20.00 

2.904 

0.9821 

7.466 

0.5020 

1.342 

9,0xl04 

8.71x10“* 

23.33 

2.592 

0.9670 

7.187 

0.4983 

1.029 

1.64x10° 

1.61x10° 

33.09 

1.863 

0.9221 

6.953 

0.4970 

0.896 

3.84x10° 

3.76x10° 

35.00 

1.844 

0,9120 

6.764 

0.4871 

0.888 

4.34x10° 

4.20x10° 

36.11 

1.739 

0.9002 

6.487 

0.4704 

0.842 

4.29x10° 

40.00 

1.646 

0.8704 

6.148 

0.4592 

0.801 
itao  asamasx 

j§gfg§ 

6.67x10° 

Table  5.1  Results  of  the  model  simulation  for  load  condition  1 
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AK  _ 
(MPaVm) 

Wp 

(J/mol) 

mOi/m01 

rnc,/mCf 

xl0‘J 

P 

£m 

(%) 

da/dN 

(mm/cyc) 

tested 

da/dN 

(mm/cyc) 

predicted 

20,00 

2.180 

0,9561 

3,908 

0.2447 

0,710 

5.46x10"* 

5.43x10"* 

23,33 

1.574 

0.9178 

3.080 

0.2290 

6.26x10"* 

25.00 

1.500 

0.8977 

2.761 

0.2099 

0.506 

1.24xl0'3 

l.lOxlO'3 

32.67 

1.431 

0,8171 

2,334 

0.1912 

2.05x1  O'3 

35,00 

1.424 

0.7896 

2.283 

0.1894 

0,482 

2.56X10'3 

40.00 

1.421 

0,7270 

2,265 

0,1818 

0.481 

3.35xl0'J 

Table  5.2  Results  of  the  model  simulation  for  load  condition  2 
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5.4  Simplified  Approach  and  Discussion 


Applications  of  proposed  model  to  these  two  loading  profiles  showed  a  good 
agreement  with  experimental  results.  This  could  be  taken  as  a  support  for  the  validity  of 
the  proposed  model.  However,  one  should  realized  that,  due  to  the  implicit  functional 
form  of  the  model,  that  an  extensive  and  complex  iterative  procedure  is  required  for  its 
execution.  An  attempt  is  thus  made  to  simplify  the  model  by  deriving  a  closed  form 
equation  of  da/dN  versus  AK  which  could  prove  to  be  easier  to  integrate,  thus  making  it 
possible  to  include  the  model  into  residual  life  prediction  procedures.  This  was  achieved 
by  obtaining  seml-empirical  equation  for  both  e m  and  ed  using  both  eq.  (5.3)  and  Fig.  5.2, 
respectively.  This  approach,  yielded  the  following  expression 


and 


(  Yl 
AK3 


n+l 


(5.18) 


ed  =  edoabrq  (5,19) 

where  Gt(f)  is  a  function  of  frequency,  n  has  the  same  meaning  as  before,  and  parameters 
a,  b  and  q  are  related  to  the  geometrical  shape  of  the  grain  boundary  ductility  curve.  It 
was  found  that  parameters  a  and  q  are  constants  with  values  of  0.9434  and  0.0494  while 
parameter  b  depends  on  both  the  frequency  and  AK  levels,  and  it  could  be  correlated  with 
frequency  and  AK  as  follows: 


(5.20) 


b  =  G2(f)AK_p 

where  G2(f)  is  a  function  of  frequency,  p  is  a  constant  with  a  value  of  1.249.  By 
invoking  the  fracture  criterion  previously  described  and  using  the  above  two  expressions, 
one  obtains 
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which  can  be  rewritten  in  the  form 


(5.21) 


-Ti=C(f)AKm  (5.22) 

dN 

where  m  is  a  constant  having  a  value  of  3.199,  and  C  is  a  function  of  frequency,  it 
should  be  emphasized  the  fact  that  the  micromechanistic  based  model  and  the  approach 
directly  based  upon  the  measurement  of  intergranular  depth  of  oxygen  diffusion  lead  to 
the  idendcal  results,  i.e.  the  Paris-type  law  of  crack  growth  expression  with  a  frequency- 
dependent  coefficient.  This,  again,  implies  that  it  is  the  environmental  influences  that 
govern  the  fadguc  crack  propagation  performance  in  the  environment-dependent  crack 
growth  stage. 

It.  must  be  emphasized  here,  however,  that  the  model  is  also  limited  in  use  to 
frequencies  smaller  than  the  transitional  frequency  of  Alloy  718  at  650°C.  Other 
constraints  of  this  model  lie  in  the  fact  that  the  model  is  valid  only  under  small-scale 


yielding  conditions  and  the  assumption  that  previous  damage  in  the  crack-tip  region,  at 
the  start  of  each  new  load  cycle,  is  neglected.  These  limitations  should  be  recognized 
when  evaluating  the  degree  of  validity  of  the  model  to  conditions  different  than  those 
described  here. 


CHAPTER  6 
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SUMMARY  AND  CONCLUSIONS 
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At  cyclic  loading  frequencies  below  that  of  the  transitional  frequency  level,  the 
elevated  temperature  fatigue  crack  growth  process  in  Alloy  718  is  viewed  to  be  fully 
environment-dependent.  Of  all  the  crack  growth  stages,  this  stage,  while  is  the  most 
critical  in  high  temperature  applications  due  to  its  highly  accelerated  crack  growth  rate, 
is  the  least  studied  or  understood.  The  objective  of  this  study  was  to  focus  on  the 
understanding  of  the  controlling  mechanisms  of  this  environment-dependent  crack  growth 
stage  in  order  to  develop  the  ability  to  predict  the  crack  growth  performance  in  Alloy  718 
under  different  loading  and  environment  conditions.  For  this  purpose,  two  major  studies 
have  been  carried  out  here. 

The  first  was  to  provide  evidence  of  the  existence  of  the  fully  environment- 
dependent  stage  in  which  the  crack  growth  rate  would  be  equal  to  the  oxygen  penetration 
rate  at  the  crack  tip.  This  has  been  achieved  by  measuring  the  depth  of  oxygen  diffusion 
at  650°C  under  the  cyclic  loading  conditions  with  and  without  hold  time  durations  at 
minimum  load  level.  As  a  result,  the  relationship  between  the  intergranular  oxygen 
diffusion  rate,  X,  and  the  value  of  the  stress  intensity  factor,  AK,  has  been  established. 
This  relationship,  when  integrated  over  the  cycle  effective  oxidation  time,  results  in  a 
closed  form  solution  describing  the  environment-dependent  fatigue  crack  growth  rate.  A 
good  agreement  was  obtained  when  a  comparison  was  made  between  the  results  of  this 
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solution  and  the  corresponding  experimental  data.  These  results  were  taken  in  support 
of  the  notion  that  environment-dependent  crack  growth  in  fully  correlated  with  the 
intergranular  depth  of  oxygen  penetration  at  the  crack  tip.  As  such,  the  environment 
dependent  crack  growth  process  is  thus  viewed  as  a  crack  tip  oxidation  phenomenon. 

The  second  study  was  to  establish  a  micromechnicai  based  quantitative  model  to 
predict  the  environmentally-dominated  crack  growth  performance.  This  model  was  based 
on  the  concept  of  the  intergranular  two-stage  crack-tip  oxidation  mechanism.  According 
to  this  mechanism,  the  oxygen  partial  pressure  controls  the  preferential  formation  of  the 
oxide  layers  at  the  cruck  tip,  and  it  suggests  that  intergranular  oxidation  under  air 
environment  occurs  in  two  stages  with  the  formation  of  a  Ni-Fe  rich  oxide  followed  by 
that  of  the  more  protective  film  formed  by  Cr203  oxide.  Thus,  the  reduction  of  the  grain 
boundary  ductility  along  the  grain  boundary  fracture  path  depends  on  the  rate  of 
formation  of  the  chromia  layer  in  relation  to  the  build-up  of  other  oxide  types  at  the  cruck 
tip.  The  determination  of  the  amount  of  Cr203  depends  on  the  amount  of  both  oxygen 
diffused  along  the  affected  grain  boundary  and  chromium  transported  via  a  mobile 
dislocation  network.  The  reduction  in  grain  boundary  ductility  due  to  oxidation  is 
balanced  by  considering  the  effective  strain  at  the  crack  tip  resulting  from  the  external 
loading.  This  balance  defines  the  fracture  criterion  of  the  model  and  permits  the 
calculation  of  the  crack  advance  per  cycle.  The  model  was  then  used  to  predict  the  crack 
growth  rate  in  Alloy  718  at  650  °C  for  two  different  loading  conditions,  and  a  good 
agreement  was  observed  when  the  data  obtained  from  the  model  in  the  form  of  da/dN 
versus  AK  were  compared  with  those  obtained  experimentally. 
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From  the  present  study,  two  major  conclusions  can  be  made: 

The  existence  of  an  environment-dependent  fatigue  crack  growth  stage  in 
Alloy  718  at  elevated  temperature  conditions  is  supported  on  the  basis  of 
comparing  the  results  of  fatigue  crack  growth  tests  in  both  air  and  vacuum 
conditions  and  also  from  comparison  between  the  experimental  fatigue  crack 
growth  results  and  those  obtained  from  a  closed  form  solution  which  is  directly 
derived  from  the  knowledge  of  the  intergranular  diffusion  rate,  X,  in  the  crack  tip 
region  over  the  cycle  effective  oxidation  time,  tox.  Two  interesting  points  should 
be  emphasized  here,  i.e. 

(a)  Within  the  environment-dependent  fatigue  crack  growth  stage,  the  cycle 
effective  oxidation  time,  should  be  token  as  the  whole  cycle  duration 
including  both  the  loading  and  unloading  part  of  the  cycle. 

(b)  The  oxygen  diffusivity  of  grain  boundaries,  Dg,  is  a  function  of  the  stress 
intensity  factor  range,  AK,  and  the  loading  frequency,  f.  Through  this 
relationship,  the  intergTanular-transgranular  transitional  frequency,  fc,  can 
identify  as  a  AK-dependent  parameter. 

The  quantitative  crack  growth  model  developed  in  this  work  on  the  basis  of  the 
two-stage  crack-tip  oxidation  mechanism  correctly  characterizes  the  environment- 
dependent  fatigue  crack  growth  behavior.  Two  assumptions,  however,  were 


made  in  order  to  simplify  the  mathematical  treatments.  The  first  is  that  the 

crack  tip  zone  is  subjected  only  to  small-scale  yielding  and  the  second  is  that  the 

initial  conditions  of  the  crack  tip  region  at  the  start  of  each  load  cycle  are  # 

identical.  These  assumptions  limit  the  applications  of  this  model  to  conditions 

which  do  not  consider  loadings  with  a  long  hold  time  at  the  maximum  load  level. 

In  these  types  of  applications,  the  spread  of  crack  tip  plasticity  due  to  creep  ® 

associated  with  hold  time  durations  may  not  justify  the  use  of  AK  as  a  crack  tip 
driving  force. 

# 
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APPENDIX  I 


MATERIAL  DESCRIPTION 

The  objective  of  this  appendix  is  to  describe  in  detail  the  aspects  of  physical 
metallurgy  of  Alloy  718.  This  includes  the  phases  of  the  microstructure,  the  properties, 
the  strengthening  mechanisms  and  the  associated  deformation  mechanisms  in  Alloy  718. 

As  a  precipitation-hardenable  nickel-chromium  alloy  developed  by  Eiselstein  [  1 22] 
in  the  late  50’s,  Alloy  718  combines  corrosion  resistance  and  high  strength  with 
outstanding  weldability.  This  alloy  contains  significant  amount  of  iron,  niobium  and 
molybdenum  along  with  lesser  amounts  of  aluminum  and  titanium,  its  nominal 
composition  Is  given  in  Table  Al.l.  The  matrix  of  Alloy  718  is  nickel  plus  Iron  fomiing 
the  austenitic  y  phase.  Elements  which  partition  largely  to  matrix  are  chromium  and 
molybdenum.  The  chromium  provides  oxidation  resistance  while  both  chromium  and 
molybdenum  contribute  to  solid  solution  strengthening.  The  niobium,  dtanium  and 
aluminum  govern  the  phase  precipitation  characteristics  of  Alloy  718  while  the  small 
additions  of  carbon  and  boron  are  intended  to  enhance  grain  boundary  strength  and 
ductility.  Addition  of  Nb  promotes  the  y"  precipitation  [123],  but  the  increase  in 
(Ti+Al)/Nb  ratio  inay  favor  y'  rather  than  y1'  precipitation  [124],  The  effects  of  several 
elements  in  nickel-base  supcralloys  are  summarized  in  Table  A1.2  [125,126]. 

I.I  Phases  In  Alloy  718 

Alloy  718  has  a  large  number  of  phases  which  either  have  characteristic 
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Table  Al.l  The  Nominal  Composition  of  Alloy  718  (wt%) 
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Table  A1.2  Role  of  Elements  in  Nickel-base  Supcralloys 


morphologies  or  will  form  in  a  specific  temperature  range.  Besides  the  FCC  y  phase 
matrix,  the  phases  normally  found  in  Alloy  718  are  disk-shaped  y"  (Ni3Nb,  Al); 
spheroidal  y'  (Ni3Al,  Ti);  a  needle/plate-like  6  phase,  also  called  p  phase  in  some 
literature  (Ni3Nb);  a  round,  island-like  Laves  phase  [(Ni,  Fe,  Cr)2(Mo,  Nb,  Ti)];  discrete 
MC  carbides;  M^C  carbides  which  may  form  as  a  boundary  film  during  service;  M23C6 
carbides  (Cr23C6)  were  also  Identified  at  grain  boundaries  [124].  The  MC  phase  forms  on 
solidification  and  remains  stable  to  about  1204°C.  The  Laves  phase  forms  on 
solidification  in  high-niobium  interdendritic  areas  and  is  stable  up  to  1177°C,  it  tends  to 
deplete  the  matrix  of  the  elements  required  for  precipitation  hardening,  thereby  reducing 
strength  [127].  The  5  phase  forms  on  cooling  during  solidification  in  high-niobium  areas 
and  is  stable  up  to  1010°C  in  wrought  and  1 121°C  in  cast  718,  precipitating  from  843°C 
to  982°C.  The  y "  phase  forms  during  cooling  or  heat  treatment  at  732°C  to  900°C.  The 
Y’  phaso  forms  on  cooling  or  when  heat  treatment  is  in  the  range  of  718°C  to  621°C. 

It  has  been  shown  that  in  Alloy  718  under  the  long  time  (>5000  Hr)  thermal 
exposure  (530°C  -  718°C)  with  stress  two  microstructural  groin  boundary  changes  will 
be  produced  [47,128,129],  The  first  change  is  the  precipitation  of  8  phase  along  the  grain 
boundaries  with  the  formation  of  a  denuded  zone  of  7"  precipitates.  The  second  change 
is  the  precipitation  of  blocky  particles  which  were  identified  us  a  phusc  rich  in  chromium, 
termed  as  aCr  phase.  This  phase  precipitated  along  grain  boundaries  and  began  to  form 
as  low  as  593°C.  The  formation  of  aCr  phase  is  very  dependent  on  the  stress  state  of  the 
material,  it  can  be  formed  at  much  shorted  time  under  creep  conditions  [47].  Radavich 
even  found  o  phusc  existed  during  long  service  conditions  at  high  temperatures.  The 


long-time  structural  transitions  in  Alloy  718  are  as  follows  [130]: 

538°C  y//  +  Y'^"  +  Y'  +  S 

593°C  y"  t  y1  -*y"  +  y*  +  S  +aCr 

621  °C  y"  +  y'-V'  +  y'S  +  aCr 

649°C  y"  +  Y/  -^Y//  +  Y;  +  8  +  aCr  +  Trace  a 

704°C  Y/y  +  y>  -fflf  +  y'  +  8  +  aCr  +  o 

The  time-temperature-transformation  (TIT)  diagram  for  Alloy  718  was  studied  by 
number  of  authors  [47,131-133]  and  the  typical  ITT  diagrams  are  shown  in  Fig.  Al.l, 
From  these  TTT  diagrams  the  following  phase  transforms  can  be  occurred.  Over  the 
temperature  range  of  about  649°C  to  871°C,  precipitation  of  y\  y"  and  6  phases  is  shown 
to  follow  y*-*  Y"  5  uansformation  sequence.  Solutions  g  of  the  Y"  phuse  begins  to 
occur  between  843°C  and  871°C.  This  event  is  associated  with  the  formation  of  large 
amount  of  5  phase.  The  y'  solvus  is  somewhat  lower  and  the  solutioning  of  this  phase 
is  complete  at  843°C.  This  TTT  diagram  also  shows  that  8  phase  forms  most  rapidly  at 
about  941°C  and  has  a  solvus  temperature  near  982°C.  Fig.  A  1.1(b)  ulso  shows  MC 
carbides,  M,C  carbides,  and  Laves  phase  that  precipitate  in  the  alloy  base  metal  over  the 
temperature  range  of  about  700°C  to  1040°C.  All  of  these  phases  go  into  solution  around 
1065°C.  Because  the  grain-boundary  phases  are  going  into  solution  find  are  therefore  no 
longer  pinning  the  boundaries,  grain  growth  occurs  rapidly  at  temperuture  about  beyond 
1040°C.  Based  on  the  knowledge  ofTTT  diagram  of  Alloy  718,  the  common  used  heat 
treatments  of  this  material  are  listed  in  Table  A  1.3. 
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Fig.  Al.l  TimeTemperature-Tram  formation  diagrams  for  Alloy  718. 
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Fig.  Al.l  Time-tomperature-trunsformation  diagrams  for  Alloy  718  (com,) 


Conventional  Heat  Treatment  (CHT) 

Solutioning: 

(Annealing) 

954°C  for  1  hr  then  air  cool  or  faster. 

:  Precipitation  Aging: 
(Duplex  Aging) 

718, 3°C  for  8  hrs,  furnace  cool  to  621. 1°C 
at  a  rate  not  to  exceed  37.8°C/hr,  hold  at 

6?.1°C  for  8  hrs  and  air  cool  to  room 
temperature. 

Modified  Heat  Treatment  (MHT) 

Solutioning: 

1093°C  for  1  hr,  furnace  step  cooled  at  38°C/hr  to 
below  538°C,  He  cooled  to  room  temperature. 

Precipitation  Aging: 

718. 3°C  for  4  hr  step  cooled  at  37.8°C/hr  to  621°C  for 

16  hr,  then  He  cooled  to  room  temperature. 

Table  A1.3  Heat  Treatments  of  Alloy  718 
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1.1.1  Strengthening  Phases 

Analysis  of  the  conventionally  heat-treated  Alloy  718  showed  that  both  y '  and  / 
precipitates  enhance  the  mechanical  properties  of  the  alloy,  the  volume  fraction  of 
Y'phnsc  has  been  observed  to  be  about  15%  while  that  of  y*  phase  was  about  4%  [135]. 
The  major  strengthening  phase  in  Alloy  718  has  been  reported  to  be  the  y"  phase  [134- 
136].  y’  phase  has  a  FCC  (H2)  structure  with  spherical  shape.  Although  a  small  amount 
of  y*  is  present,  it  is  believed  that  any  strength  contributed  by  this  phase  is  incidental 
[134],  The  y"  phase  is  coherent  and  has  an  DO:a-ordered  body-ccntered-tetragonal  (BCT) 
structure.  The  stoichiometry  of  y"  is  given  as  Ni5(Nb,  Ti,  Al)  [137,138].  The  lattico 
parameters  for  y"  phase  are  Bq  -  0.3624  nm  and  c0  *  0.7406  nm,  and  the  unit  cells  for 
y’  and  y"  are  shown  in  Fig.  A  1.2.  For  the  convendonal  heat-treated  Alloy  718,  the 
volume  fraction  of  the  strengthening  phase  (primarily  y")  is  about  20  %  [135],  The  y" 
precipitates  as  dlsc"shaped  particles  and  with  the  following  orientation  relationship: 
(001)/ II  (001)7  and  [1001/11  <100>y 

and  the  c  axis  is  normal  to  the  plane  of  the  disc.  Kiiinan  und  Warrington  [136]  have 
pointed  out  that  the  matrix  /  /  precipitate  misfit  is  concentrated  in  the  direction  parallel 
to  the  c  axis,  and  they  estimated  thut  this  misfit  is  about  2.5%  whereas  that  in  the  plane 
of  the  plate  is  only  about  0.8%.  It  has  been  observed  [124]  that  the  /  panicles  arc 
always  seen  adjacent  to  the  y’  precipitates,  and  they  believed  that  y'  panicles  are 
potential  nuclearion  sites  for  /  phase.  Examination  of  the  two  unit  cells  drawn  on  Fig. 
A  1.2  shows  that  the  yVmatrix  interfaces  can  be  very  favorable  nuclearion  sites  for  / 
phase.  Except  for  the  measured  slight  distortion  along  the  a  axis  and  the  higher  distortion 


along  the  c  axis  of  D022-ordered  structure,  the  two  phases  represent  ordered  variations  of 
the  FCC  lattice.  In  the  y*  phase,  Nb  (Al,  Ti)  atoms  are  located  at  the  cube  comers,  while 
in  the  y"  phase,  Nb  atom  occupies  the  center  and  the  comers  of  a  cell  are  formed  by 
putting  two  FCC  cells  one  upon  the  other.  The  sizes  of  y*  and  y"  phases  for  different 
heat  treatment  conditions  are  listed  in  Table  Al,4  [135],  During  prolonged  thermal 
exposure,  the  growth  of  y*  and  y"  were  observed  by  several  investigators  [9,16,1351. 
Paulonis  et  al  [135],  Pedron  and  Pineau  [139]  reported  that  there  was  no  apparent  change 
in  volume  fraction  and  only  slight  coarsening  of  the  y"  phase  (5%-10%)  during  100  hr 
at  649°C.  However,  considerable  coarsening  of  the  y"  phase  occurred  during  the  100  hr 
age  at  760°C,  and  the  particle  volume  fraction  is  sufficiently  lower  [135]. 

It  worth  mentioning  that  a  peculiar  precipitation  microstructure,  compact 
morphology,  has  been  Identified  by  Cozar  and  Pineau  [124].  It  appears  to  be  cube-shaped 
y'  particles  coated  with  a  y"  shell  over  their  six  faces.  This  morphology  has  proved  to 
have  a  very  slow  rute  of  coarsening.  The  condition  to  achieve  this  morphology  requires 
a  minimum  (Tl+Al)/Nb  ratio  between  0.9  and  1.0.  Moreover  the  ageing  temperature  has 
to  be  chosen  in  such  a  way  that  initially  isolated  y’  precipitates  reach  a  size  of  about 
200A  before  production  of  y"  precipitates. 

Oblak  et  al  [134]  studied  the  ageing  effect  on  the  y"  phase  under  stress  conditions. 
They  found  that  the  y"  precipitate  forms  as  ellipsoidal  discs  lying  on  ( 100) .  Because  the 
coherent  y"  panicles  produce  a  tetragonal  distortion  of  the  matrix,  the  specific  variants 
present  can  be  controlled  by  application  of  stress  during  ageing.  For  example,  only  the 
[001]  y"  variant  is  formed  when  a  tensile  stress  is  applied  along  [001].  Reversing  the 


7’  (diameter  of  sphere) 


CHT 

250  A 

CHT+649°C/100  hr 

300  A 

CHT+760°C/1U0  hr 

600  A 

CHT+843°C/16  hr 

750  A 

CHT+871°C/16  hr 


7"  (diameter  of  disc) 
(thickness:  50-90  A) 

600  A 
725  A 

_ 3000  A 

3500  A 
5500  A 


Table  A1.4  Dimension  of  7’  and  7"  in  Alloy  718 


stress  suppressed  the  [001]  7"  variant  and  promotes  formation  of  the  [100]  and  [010] 
variants.  The  effect  upon  strength  is  significant.  The  critical  resolved  shear  stress  of  the 
former  structure  is  25°’-  ugher  than  that  of  standard  aged  material,  while  that  of  the  latter 
is  6%  less.  This  stress-enhanced  coarsening  of  the  7"  precipitates  also  observed  by 
Radavich  [128]  and  Molins  et  al  [129]. 

1.1.2  8  Phase 

Both  7’  and  7"  phases  in  Alloy  718  are  metastable  [135],  while  the  7"  phase  is 
less  stable  phase  than  the  7’  phase  [139].  Upon  prolonged  thermal  exposure  the  7"  phase 
will  be  replaced  by  the  stable  precipitate  8  phase.  It  has  an  incoherent  orthorhombic 
structure  (DOt)  with  a  stoichiometry  of  Ni3Nb  [16,140].  The  8  phase  precipitates  as  laths 
(in  a  Widmenstatten-type  morphology)  along  { 1 1 1 } 7 planes  with  the  following  orientation 
relationship: 

(010)5||  {111)7  and  [100]6  ||  <110>y 

The  8  phase  may  also  form  by  a  cellular  precipitation  along  grain  boundaries  [136]. 
Cozar  and  Pineau  [124]  have  pointed  out  that  the  metastable  y1'  particles  are  dissolved 
near  the  8-Ni3Nb  lathe  and  only  7’  phase  remains.  Thus,  formation  of  8  phase  is 
detrimental  to  the  strength  of  the  material  [139],  Moreover,  the  7"  ->  8  phase 
transformation  is  largely  accelerated  by  prior  deformation,  since  it  has  been  shown  that 
planar  defects  present  in  the  7"  particles  after  deformation  can  be  act  as  efficient  nuclei 
for  the  precipitation  of  the  stable  8  phase  [15,141].  This  concept  has  been  applied  to 
develop  a  etching  technique  for  measuring  the  monotonic  zone  interface  to  the  fracture 
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surface  [140]. 

Additional  to  the  amount,  Wilkinson  [142]  found  that  the  morphology  of  the  8 
phase  is  important  to  the  mechanical  properties.  His  work  showed  that  during  forging  5 
phase  tends  toward  spheroidization,  which  appears  to  be  a  preferred  morphology.  If  strain 
is  occurring  in  a  temperature  range  where  5  phase  is  precipitating,  then  growth  of  this 
phase  tends  to  be  discontinues  and  occurs  in  short  lath  type  particles  within  and 
emanating  from  the  grain  boundaries,  growing  into  the  grains  along  crystallographic 
planes.  If  material  is  heated  into  a  temperature  range  for  5  phase  precipitation  and  no 
deformation  or  inadequate  deformation  is  carried  out,  then  8  precipitation  tends  to  occur 
in  a  more  continuous  fashion  along  the  grain  boundaries  and  in  long  needles  across  the 
grain  width.  This  morphology  has  been  shown  to  have  an  adverse  effect  on  stress  rupture 
ductility  and  fracture  toughness.  Brooks  and  Bridges  [47]  obtained  similar  conclusion, 
and  they  believed  that  a  microstructure  which  contains  on  even  distribution  of  discrete 
grain  boundary  8  particles  is  required  to  prevent  notch  sensitive  stress  rupture  behavior. 

The  effects  of  8  phase  on  the  mechanical  properties  are  quite  subtle.  As  is  just 
mentioned,  deformation  of  Alloy  718  at  elevated  temperatures  promotes  the  formation  of 
8  platelet.  Presumably  the  deformed  regions  in  the  y  phase  act  as  nuclei  for  8  phase 
transformation.  The  local  re-ordering  of  the  y  lattice  along  with  the  similarity  of  the  5 
and  y  compositions  should  provide  for  a  rapid  nucleation  event  as  well  as  rapid  growth 
since  only  short  range  diffusion  would  be  required  to  convert  y  phase  into  5  phase.  The 
formation  of  deformation-induced  8  phase  will  naturally  reduce  the  strength  since  Nb  is 
drawn  from  the  principal  strengthening  phase,  y".  However,  the  effect  of  8  phase  is  not 
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totally  negative,  since  at  grain  boundaries  it  can  provide  some  oxidation  resistance  or  by 
providing  oxygen  traps  through  an  increase  in  the  boundary  area  [35,143].  Such  a 
mechanism  would  tend  to  reduce  the  effective  diffusion  rate  and  to  require  higher 
effective  concentrations  of  oxygen  in  order  to  form  deleterious  oxides.  The  ready 
formation  of  8  phase  by  deformation  processes  at  650°C  and  above  is  a  major 
consideration  for  Alloy  718  in  limiting  its  use  to  temperature  below  650°C.  Another  use 
of  the  8  phase  to  pin  down  the  boundaries  during  forging  is  well  known  [139],  For 
wrought  Alloy  718  processed  below  the  5  solvus,  grain  boundary  migration  is  inhibited 
by  the  presence  of  this  phase  and  microstructures  retain  a  fine  grain  size.  Subsequent 
thermal  exposure  above  this  temperature  dissolves  the  8  phase  and  allow  more  rapid  groin 
growth. 

1.1.3  Laves  Phase 

Laves  phase  which  is  a  hexagonal  A2B  type  intermetallic  compound  has  been 
identified  in  Alloy  718.  The  usual  fonnulae  of  this  phase  have  been  found  to  be 
Fe2(Nb,Ti,Mo),  Co2Ti  and  (Ni03Fe0i23Cr0,M)2(NbMoSi)  [144].  The  lattice  parameter  was 
determined  to  be  ay  =  0.47  nm  and  c0  =  0.75  nm.  Burke  et  al  [144]  found  that  this  phase 
formed  either  during  the  solidification  process  or  during  the  anneal  at  1038°C. 
Investigations  [37,139]  have  shown  that  the  Laves  islands  are  the  crack  initiation  sites. 
The  presence  of  the  Laves  phase  in  wrought  materials  is  a  result  of  incomplete 
homogenization  of  the  alloy  prior  to  final  heat  treatments.  It  has  been  found  that  the 
formation  of  the  Laves  phase  requires  10-12%  Nb,  the  8  phase  requires  6-8%  Nb,  they 


phase  needs  4%  Nb,  while  the  y’  phase  can  form  with  Nb  level  below  4%  [139],  For 
cast  Alloy  718,  the  7  matrix  rejects  selectively  the  large  atoms  of  elements  such  as  Nb, 
Mo  and  Ti  to  the  interdendritic  regions  during  solidification.  These  elements  can  form 
MC  carbide  (NbC)  and  the  Laves  phase  in  the  high  Nb  areas  as  the  metal  solidifies  from 
about  1260°C  to  about,  1093°C.  In  this  case,  the  alloy  is  heavily  segregated  and  consists 
of  two  type  of  regions:  the  dendritic  regions  which  are  concentrated  in  Fe,  Cr  and  Ni 
while  the  interdendritic  areas  are  very  rich  in  Nb,  Mo  and  Ti.  Thus,  the  Laves  phase 
forms,  in  general,  preferentially  In  the  interdendritic  region. 

It  should  be  noted  that  the  Laves  phase  contains  about  25%  Nb  and  the  formation 
of  this  phase  depletes  the  content  of  Nb  in  matrix  which  is  needed  for  7"  and  7’ 
precipitation,  thus  the  presence  of  the  Laves  phase  is  adverse  to  mechanical  properties. 
Practically,  in  order  to  reduce  the  presence  of  the  Laves  phase  arid  produce  a  unifoim 
material  -for  forging,  the  homogenization  process  is  critical.  Basically,  this 
homogenization  process  consists  of  the  Laves  phase  solutioning  and  uniform  Nb 
distribution.  Solution  of  the  Laves  phase  is  to  carefully  heat  up  to  the  Laves  solvus 
(about  1200°C)  and  hold  there  for  a  period  of  time  (e.g.  2  hrs),  then  rapidly  cooling  from 
solutioning  temperature  to  prevent  the  reprecipitation  of  the  Laves  phase.  To  uniform  the 
Nb  distribution,  a  duplex  aging  treatment  is  usually  applied  [37],  Current  technology  of 
homogenization  process  includes  two  types:  thermal  treatment  and  thcrmal-mechunical 
treatment  (hot  isostatic  pressing  or  HIP).  The  work  of  Smith  and  Michel  [37]  showed 
that  both  thermal  treatment  and  HIP  can  improve  high  temperature  fatigue  crack  growth 
resistance  remarkably. 


A 15 


1.1.4  Castries 


V  ^irbides  are  present  in  Alloy  718  [124,145,146],  The  MC  carbides 
(usually  (Ti,Nb)C)  have  a  cubic  structure  with  a  lattice  parameter  a^  -  0.43  -  0.47  nrn. 
They  are  very  stable  compounds  and  are  formed  during  cooling.  They  exist  in  blocky  or 
script  morphologies  and  are  distributed  heterogeneously  throughout  the  alloy,  both  in 
intergranular  and  trunsgranular  positions,  often  intcrdendritically.  The  MC  carbides  are 
a  major  source  of  carbon  for  subsequent  phase  reactions  during  heat  treatment  and 
service.  It  should  be  mentioned  that  carbonitrides  (  (Ti,Nb)(C,N)  )  sometimes  are  also 
classified  into  MC  type  [145,146], 

The  M<,C  carbides  have  a  FCC  structure  with  a  lattice  parameter  Uo  =  1.085  -  1 . 175 
ntn.  The  usual  formulae  of  M^C  are  Fc3Mo3C,  Fe3Nb3C  and  Nb3Co3C.  They  can  be 
formed  during  heat  treatment  or  service  in  material  containing  relatively  high 
molybdenum  content  and  can  also  be  formed  from  degeneration  of  MC  carbides.  They 
usually  have  a  blocky  intergranular  morphology. 

The  M23CS  carbides,  with  a  FCC  structure  and  a  lattice  parameter  a0  =  1.05  -  1.07 
nrn,  ore  formed  during  lower- temperature  heat  treatment  and  service  (760°C  to  980°C) 
from  degeneration  of  MC  carbides  and  from  soluble  residual  carbon  in  the  alloy  matrix. 
The  common  compositions  of  M23Ca  are  Cr21Mo2C6,  Cr23C6  and  (Cr,Fe,Mo)23C6.  Although 
usually  seen  ut  grain  boundaries  as  a  platelet  morphology,  they  occasionally  occur  along 
twin  bands. 

The  M7Cj  carbides  have  a  hexagonal  structure  with  lattice  parameters  ^  -  1.398 
nm  and  c0  =  0.4523  nm.  They  typically  have  blocky  intergranular  morphology,  and  the 


most  common  formula  of  M7C3  is  Cr7C3. 

In  addition  to  the  y"  phase  strengthening,  investigators  [47,122,124,137,140,146] 
have  suggested  that  niobium  carbide  (NbC)  precipitation  may  also  enhance  the  properties 
of  Alloy  718.  Carbide  type,  distribution  and  morphology  have  been  shown  to  affect 
properties  of  this  alloy,  particularly  ductility  and  stress  rupture  [137]. 

1.2  Mechanical  Properties  of  Alloy  718 

Since  Alloy  718  has  been  used  in  all  kinds  of  applications  which  involves  static 
and  cyclic  loadings  and  where  exposure  at  elevated  temperatures  for  extended  times  is 
required,  thus  the  mechanical  properties  of  this  material,  including  the  behaviors  of 
tensile,  creep,  low  cycle  fatigue,  fracture,  crack  growth  and  time-dependent  notch 
sensitivity,  has  been  studied  extensively.  Also,  some  important  physical  properties  of 
Alloy  718  are  summarized  in  Table  A1.5  [125]  and  Fig.  A1.3  [125]. 

1.2.1  Tensile  Strength  and  Ductility 

Tensile  strength  and  ductility  are  the  most  basic  parameters  for  the  design  concern. 
Thus,  a  large  number  of  research  work  associated  with  these  properties  hus  been  carried 
out  and  the  results  are  listed  in  Table  A1.6  [14,15,147,148]. 

The  Young’s  modulus,  E,  torsional  modulus,  G,  and  Poisson’s  ratio  v  ore  tabulated 
In  Tabic  A  1.7  [120,149]. 

Fournier  et  al  [150]  and  Clavel  et  al  [13]  also  studied  the  cyclic  tensile  behavior 
of  Alloy  718  at  25°C  and  650°C.  They  recorded  the  changes  in  load  during  the  course 
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Fig.  A1.3  Variations  of  coefficient  of  thermal  expansion  with  temperature. 
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of  the  fatigue  tests.  Fig.  A  1.4  reported  the  variation  of  the  cyclic  tensile  stress  with  the 
number  of  cycles.  This  figure  shows  that,  at  room  temperature,  Alloy  718  gives  rise  to 
a  short  period  of  rapid  hardening  followed  by  softening.  The  maximum  load  is  obtained 
at  a  certain  number  of  cycles  which  is  smaller  for  the  higher  the  applied  strain  amplitude. 

At  elevated  temperature,  the  evolution  of  the  load  with  the  number  of  cycles 
was  found  to  be  different  from  that  encountered  at  room  temperature.  The  period  of 
hardening  was  absent  and  the  specimens  soften  from  the  first  cycle.  A  comparison  of  the 
cyclic  behavior  of  specimens  tested  at  various  temperatures  is  plotted  In  Fig.  A1.5. 

The  comparison  of  the  monotonic  and  cyclic  tensile  properties  of  Alloy  718  can 
also  be  made  using  monotonic  and  cyclic  stress-strain  curves,  they  are  shown  in  Fig. 
A1.6.  These  results  indicated  that,  at  room  temperature,  Alloy  718  gives  rise  to  softening 
when  applied  plastic  strain  amplitude  is  lower  than  about  1%.  Above  this  value  this  alloy 
experiences  a  slight  hardening  as  compared  to  the  monotonic  tensile  curve.  At  550°C  and 
650°C,  only  softening  is  observed  in  all  the  plastic  strain  range  investigated.  This 
softening  is  quite  significant  020%)  and  is  larger  than  that  measured  at  room 
temperature.  The  cyclic  softening  behavior  of  Alloy  718  is  related  to  its  peculiar 
strengthening  mechanism,  and  this  will  be  discussed  later. 

Tensile  properties  of  Alloy  718  thermully  aged  up  to  50,000  hours  were 
investigated  by  Korth  and  Trybus  [150],  The  tensile  properties  of  Alloy  718  at  room 
temperature  after  long  term  thermal  exposure  ure  listed  in  Table  A1.8  [150].  Fig.  A1.7 
shows  the  effects  of  long  term  aging  at  593°C  and  649°C  on  yield  and  ultimate  tensile 
strength,  and  Fig.  A1.8  illustrates  the  effects  on  ductility.  The  results  indicated  an 
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Fig.  A1.4  Evolution  of  the  cyclic  tensile  stress  with  the  number  of  cycles  at 
25°C.  Numbers  denote  plastic  strain  amplitude. 
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Fig.  A1.S  Comparison  of  the  cyclic  stress-strain  response  of  Alloy  71 8  at  various 
temperatures  for  a  plastic  strain  amplitude  of  0.42%  and  a  frequency 
of  0.05  Hz. 
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Fig.  A1.6  Comparison  of  the  cyclic  stress-strain  curves  and  monotonic  tensile 
curves  of  Alloy  718  at  25°C,  550°C  and  6.*>0°C. 
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Fig.  A  1.7  Effects  of  long  term  thermal  aging  of  Alloy  71.8  on  the  room 
temperature  (a)  yield  strength  and  (b)  the  ultimate  strength. 
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Fig.  A1.8  Effects  of  long  thermal  aging  of  Alloy  718  on  the  room  temperature  (a) 
total  elongation  and  (b)  reduction  in  area. 
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increase  the  hardening  after  593°C  aging  up  to  50,000  hr,  with  most  of  it  occurring  in  the 
first  few  thousand  hours.  Aging  at  649°C  results  in  an  increase  in  strength  in  the  first 
few  thousand  hours  of  aging  and  then  a  gradual  softening  is  observed. 

1.2.2  Crack  Growth  Rate 

A  detailed  review  on  fatigue  crack  growth  rate  of  Alloy  718  is  already  given  in 
Chapter  1,  so  in  this  section  the  attention  will  be  focussed  on  the  creep  crack  growth  rate 
of  Alloy  718.  The  typical  creep  crack  growth  rates  versus  stress  intensity  factor  K  for 
538°C  and  650°C  are  plotted  [151,152]  in  Figs.  A1.9  and  A1.10.  The  grain  size  effect 
on  creep  crack  growth  in  Alloy  718  at  650°C  is  illustrated  by  Liu  et  al  [152]  and  Pedron 
et  al  [9]  in  Figs.  A  1.11  and  A  1.1 2.  Research  work  has  indicated  that  the  distribution  of 
the  second  phase  in  Alloy  718  could  influence  the  creep  crack  growth  behavior.  Pedron 
et  al  [9]  found  that  the  presence  of  the  5  phase  along  the  groin  boundaries  leads  to  an 
improvement  in  crack  growth  resistance  (shown  in  Fig.  A  1.1 3), while  Liu  et  al  [152] 
observed  that  carbide  distribudon  also  alters  the  creep  crack  growth  behavior  significantly 
(see  Fig.  A1.14).  They  believed  that  the  carbide  distribution  determines  the  local  fracture 
toughness  and  affects  the  creep  crack  growth  behavior  of  Alloy  718  at  650°C  via  the 
tearing  process. 

1.2.3  Fatigue  Strength  and  Life 

A  particular  attractive  property  of  this  alloy  is  its  fatigue  strength  at  elevated 
temperatures.  A  number  of  studies  [14,148,149]  have  been  conducted  on  fatigue  strength 
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Fig.  A1.9  Creep  crack  growth  rate  as  a  function  of  stress  intensity  factor  K  at 
650°C. 
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Fig.  A1.10  Creep  crack  growth  rate  as  a  function  of  stress  intensity  factor  K  at 
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Fig.  Al.ll  Grain  size  effect  on  creep  crack  growth  in  Alloy  718  at  650°C. 
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Fig.  A1.12  Creep  crack  growth  rate  as  a  function  of  microstructure. 
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Fig.  A1.13  Effects  of  presence  of  intergranular  6  phase  on  creep  crack  growth  rate 
at  650°C. 
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Fig.  A1.14  Effects  of  carbide  distribution  on  creep  crack  growth  rate  at  650°C. 


A34 


and  life  of  Alloy  718.  Over  700  fatigue  tests  on  several  heats  of  Alloy  718  [149]  were 
performed  in  Idaho  National  Engineering  Lab  at  various  temperatures  and  the  loading 
frequencies  are  from  1  Hz  to  167  Hz.  The  results  are  shown  in  Fig,  A1.15. 

Clavel  et  al  [14]  and  Sanders  et  al  [148]  have  examined  the  property  of  low  cycle 
fatigue  in  Alloy  718.  The  variation  in  the  fatigue  life  Nf  with  the  plasric  strain  Ae,,  or 
plastic  strain  amplitude  Ae^/2  at  various  temperatures  are  illustrated  in  Figs.  A1.16  and 
A  1.17.  It  should  be  noticed  that  two  distinct  regions  of  life  are  existed. 

Sanders  et  al  [148]  discussed  the  possible  reason  for  the  presence  of  a  break  in  the 
curve.  They  thought  that  the  nature  of  deformation  process  can  be  strongly  affected  by 
plastic  strain  amplitude,  temperature  and  frequency.  They  postulated  that  the  cyclic 
plastic  behavior  below  538°C  is  dominated  by  deformation  twinning,  and  above  538°C 
the  predominate  mode  of  deformation  is  slip.  Thus,  the  specimens  tested  above  538°C 
will  have  different  strain-life  response  than  those  tested  below  that  temperature, 

1.2.5  Fracture  Toughness 

The  fracture  toughness  is  one  of  important  mechanical  properties  for  design 
concern.  The  room  temperature  and  elevated  temperature  Jl0  and  R-curve  of  conventional 
and  modified  Alloy  718  are  illustrated  in  Figs.  A  1.1 8,  A  1.1 9  and  A  1.20  and  tabulated  in 
Table  A1.9  [147], 

These  results  indicated  that  fracture  resistance  for  the  MHT  lots  was  consistently 
superior  to  that  for  their  CHT  counterparts,  and  the  further  work  of  Mills  [147]  found  that 
fracture  toughness  for  the  CHT  material  was  independent  of  temperature,  whereas 
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Fig.  AL15  Fatigue  data  of  Alloy  718  at  538CC  (cone). 


Fig.  A  1.15  Fatigue  data  of  Alloy  718  at  593°C  (conL). 
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Fig.  A1.16  Variation  in  fatigue  life  N,  with  plastic  strain  amplitude  AEp/2  at 
various  temperatures. 
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Fig.  A.  1.18  R-curves  for  conventional  and  modified  Alloy  718  tested  at  room 
temperature. 
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Fig.  A  1.19  R-curves  for  conventional  and  modified  Alloy  718  tested  at  427°C. 
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toughness  values  for  the  MHT  materials  decreased  by  20%  as  temperature  increased  from 

24°C  to  538°C.  Fractographic  examinations  revealed  [147,153]  that  primary  microvoids 

nucleated  at  broken  carbonitrides  during  the  initial  stages  of  plastic  strainning,  regardless  # 

of  heat  treatment.  The  inferior  toughness  for  the  CHT  alloy  was  associated  with  the 

premature  formation  of  a  second  population  of  dimples  that  were  nucleated  by  coarse  8 

particles.  The  MHT  alloy  dissolved  the  coarse  8  particles  which  suppressed  secondary  ^ 

microvoid  formation  and  thereby  enhanced  the  fracture  resistance.  Mills  [147]  also 

studied  the  effect  of  aging  on  the  fracture  toughness.  He  found  that  aging  at  566°C  for 

33,000  hours  caused  a  30%  reduction  in  JIa  at  24°C  and  538°C  for  both  the  CHT  and  ^ 

MHT  materials. 


1.2.5  Stress  Rupture  and  Time-Dependent  Notch  Sensitivity 

The  stress  rupture  and  time-dependent  notch  sensitivity  are  important  properties 
of  this  material  thus  have  been  investigated  by  many  researchers  [16,47,154-159].  The 
normal  stress  versus  time  of  failure  at  different  temperatures  are  shown  in  Fig.  A  1.21. 
It  was  found  [143]  that  exposure  at  575°C  for  10,000  hr  then  creep  tested  at  705°C  with 
a  stress  of  415  MPa  exhibited  a  longer  life  than  the  uncx.posed  material.  However,  the 
material  exposed  at  650°C  showed  a  significant  decrease  in  the  life  when  tested  under  the 
same  conditions.  Similar  results  also  obtained  from  the  work  of  Molins  et  al  [129]. 

The  time  •dependent  notch  sensitivity  is  the  susceptibility  for  premature  crack 
initiation  to  occur  at  elevated  temperatures,  in  the  vicinity  of  a  notch  or  stress  raiser,  as 
a  function  of  time  at  a  constant  external  stress  below  the  yield  stress  as  compared  to  an 
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Fig.  A1.21  Creep-rapture  data  for  Alloy  718  at  various  temperatures. 


identical  unnotched  specimen.  The  occurrence  of  time-dependent  notch  sensitivity  was 
due  to  more  rapid  initiation  of  intergranular  cracks  in  notched  than  in  smooth  specimens. 
This  was  associated  with  the  local  deformation  at  the  base  of  the  notches  accompanying 
the  relaxation  of  the  stress  concentrations. 

Wilson  [16]  extensively  studied  the  time-dependent  notch  sensitivity  of  Alloy  718. 
In  his  work,  the  results  for  the  material  heat  treated  10  hr  at  1066°C  plus  48  hr  at.  732°C 
provide  an  example  of  severe  time-dependent  notch  sensitivity  (shown  in  Fig.  A1.22). 
At  the  shorter  times  at  the  lower  test  temperatures,  the  rupture  strengths  for  notched 
specimens  were  somewhat  lower  than  those  for  smooth  specimens,  At  varying  time 
periods,  the  notched  specimen  rupture  curves  exhibited  drastic  increases  in  steepness  so 
that  the  rupture  strengths  for  notched  specimens  decreased  with  time  to  values 
considerably  below  those  obtained  in  tensile  tests,  i.e.  the  material  exhibited  time- 
dependent  notch  sensitivity.  The  results  for  the  material  heat  treated  lOhr  at  927°C  plus 
48  hr  at  732°C  are  in  contrast  to  the  foregoing  behavior  and  are  typical  of  the  heat  treated 
material  not  susceptible  to  time-dependent  notch  sensitivity  (shown  in  Fig.  A  1.23). 

Wilson  [16]  found  that  variations  in  notch  sensitivity  were  correlated  with  changes 
in  the  dislocation  motion  mechanism,  the  y'/y"  panicles  smaller  than  a  critical  size  were 
sheared  by  dislocations.  This  gave  rise  to  localized  defonnation  and  time-dependent 
notch  sensitive  behavior.  Lager  particles  were  bypassed  by  dislocations,  under  these 
conditions,  the  deformation  was  homogeneous  and  no  time-dependent  notch  sensitivity 
was  observed.  He  further  pointed  out  that  decreasing  the  solution  temperature  or 
increasing  the  time  and/or  temperature  of  the  aging  treatment  decreased  the  susceptibility 
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Fig.  AU2  Larson-Miller  pararaters  curves  showing  rime-dependent  sensitivity. 


to  the  time-dependent  notch  sensitivity. 

Feldstein  and  Mendoza  [158]  obtained  similar  results,  but  they  thought  that  time- 
dependent  notch  sensitivity  is  due  to  the  way  dislocations  break  from  their  precipitate 
barrier  in  reference  to  a  average  critical  distance  between  the  precipitates.  When  the 
average  distance  between  precipitates  is  less  than  or  equal  to  the  critical  distance, 
dislocations  will  shear  through  the  coherent  particles  thus  will  create  localized  plastic 
deformation.  When  this  is  combined  at  identical  locations  with  high  temperature 
creep  deformation  a  crack  will  nucleate  prematurely,  and  the  notch  sensitivity  will  occur. 
Oil  the  other  hand,  when  the  average  distance  between  precipitates  Is  greater  than  the 
critical  distance,  dislocations  will  pass  between  the  precipitates  without  cracking  them  and 
the  notch  sensitivity  does  not  occur. 

1.3  Deformation  Behavior  of  Alloy  718 
1.3.1  Mode  of  deformation  in  Alloy  718 

Wilson  [16]  has  examined  the  deformation  behavior  of  Alloy  718,  and  he  observed 
that  the  precipitate  particle  size  (y’/y")  governs  the  mechanism  of  dislocation  motion  in 
this  material.  When  y’/y"  particles  are  smaller  than  a  critical  size,  they  tire  sheared  by 
dislocation  and  this  gives  rise  to  a  localized  deformation.  For  the  particles  larger  than  a 
critical  size,  they  are  bypassed  by  dislocation.  Under  this  condition,  the  deformation  is 
homogeneous.  He  found  that  the  critical  size  is  associated  to  several  factors,  such  as  the 
fraction  of  y ’/y"  precipitate,  ageing  time  and  temperature.  Decreasing  the  amount  of 
y’/y"  precipitate  lowers  the  critical  size,  while  longer  ugeing  time  or  higher  ageing 


temperature  will  increase  the  critical  size.  Also,  the  critical  size  decreased  with 
decreasing  solution  temperature.  Fournier  and  Pineau  [15]  conducted  similar 
investigation.  In  their  work,  the  microstructure  induced  by  cyclic  deformation  was 
examined.  They  noticed  that  shear  of  y"  precipitates  occurs  during  the  course  of  cyclic 
straining.  The  shearing  process  leads  to  faults  in  the  precipitates.  Plastic  deformation 
proceeds  very  heterogeneously  by  the  propagation  of  intense  deformation  bands  which 
were  identified  as  mechanical  twins.  They  also  explained  that  the  cyclic  softening  in 
Alloy  718  is  attributed  primarily  to  the  shear  of  precipitates  which  locally  alters  their 
crystallographic  structure  from  a  BCT  metastable  structure  (y"  phase)  to  a  orthorhombic 
stable  structure  (5  phase)  and  secondly  to  the  reduction  of  the  mean  particle  size  in 
intense  deformation  bands. 

Using  TEM  investigations  by  Sanders,  Frishmuth  and  Emblcy  [148]  confirmed  the 
presence  of  microtwins  which  formed  during  the  deformation  process.  Principally, 
twinning  was  confined  to  the  lower  temperatures  204,  316  and  427°C.  The  density  of 
microtwins  appeared  to  increase  with  increasing  plastic  strain  range  as  reported  by 
Fournier  and  Pineau  [15].  At  higher  temperatures  (538°C  and  649°C),  the  primary  mode 
of  deformation  was  thought  to  be  by  slip.  When  slip  bands  occurs,  incoherent  particles 
were  often  associated  with  these  bands,  suggesting  that  dislocations  were  generated  at  the 
interface  between  these  particles  and  the  mauix.  Occasionally,  however,  a  favorably 
oriented  grain  may  twin  at  the  higher  temperatures.  If  after  the  twinning  event  occurs, 
successive  deformation  leads  to  slip,  which  is  the  predominate  mode  of  deformation,  the 
twin  will  become  damaged.  The  damage  of  the  twin  does  not  occur  uniformly  since  the 
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growth  of  a  deformation  twin  generates  lattice  defects  which  are  not  uniformly  distributed 
in  the  vicinity  of  the  twin  interface. 

Investigations  of  Alloy  718  by  Fournier  and  Pineau  [15]  have  provided 
considerable  information  relative  to  the  effect  of  frequency  as  well  as  temperature  and 
plastic  strain  amplitude  on  the  deformation  mechanism.  In  their  work,  microtwins  were 
observed  at  room  temperature  and  550°C  as  well.  They  noticed  that  more  intense 
deformation  bands  formed  by  twinning  were  accumulated  at  high  temperatures  (550°C) 
than  at  low  temperatures.  However,  deformation  by  twinning  was  promoted  by  a  decrease 
in  frequency  at  high  temperatures.  The  results  from  Sanders  et  al  [148]  supported  this 
conclusion  since  very  little  twinning  was  observed  at  a  frequency  of  0.33  Hz.  This  range 
of  behavior  from  no  twinning  to  intense  deformation  bands  was  exactly  what  has  been 
observed  in  the  work  of  Clavel  and  Pineau  [14].  At  high  stress  intensity  factor  range  AK, 
no  twinning  was  observed  for  frequencies  of  0.67  Hz  and  above,  whereas  the  intense 
deformation  bands  were  noted  for  frequencies  of  0.05  Hz  and  below.  They  concluded 
that  the  deformation  mode  at  538  and  649°C  was  characterized  as  slip  bands,  however, 
the  deformation  mechanism  shifts  form  slip  band  deformation  at  high  frequency  (0.33  Hz) 
to  increasing  twin  band  formation  at  low  frequency  (0.05  Hz). 

1.3.2  Deformation  Mechanisms  of  7"  Phase 

The  deformation  mechanisms  as  well  as  die  deformation  microstructures  pertaining 
to  7"  strengthened  alloys  have  been  investigated  extensively  [15,134,136,148,160,161]. 
For  7"  D022-ordered  structure,  the  close-packed  (112)y  planes  are  governed  by  the 


stacking  rule  that  all  Nb  atoms  shall  be  in  contact  with  all  Ni  triangular  interstices  in 
neighboring  layer,  so  that  each  Nb  atom  is  coordinated  with  Ni  atoms  only.  In  order  to 
conform  with  this  rule,  the  superlattice  must  span  two  three-layer  stacking  cycles,  i.e. 
these  planes  are  stacked  in  a  six-layer  sequence,  ABCDEFABCDEF,  to  make  the  D022- 
ordered  structure.  Electron  microscopy  observations  for  deformed  Alloy  718  specimens 
revealed  that  the  Y'  precipitates  were  sheared  by  dislocations  and  that  the  cutting 
mechanisms  lead  faults  within  7"  particles,  These  faults  are  frequently  observed  as  bands 
lying  along  close-packed  (112)7  planes.  Several  researchers  [15,16,136,160,162]  have 
suggested  that  these  faults  are  formed  by  a  dislocation  reaction  at  the  y"  precipitates 
interface  as,  for  example 

i[110]-»i[2U]  +  i[I2T] 

2  0  0 

and  the  passage  of  the  second  partial  dislocation  through  the  D022-ordered  structure.  This 
shearing  mechanism  changes  that  six-layer  stacking  sequence  of  they"  particles  and  gives 
the  lower  energy  stacking  sequence  ABCDEFABABCDEF  which  is  associated  with  the 
stable  5  phase.  It  is  worth  noting  that  the  observation  of  faults  inside  the  precipitates  by 
dark  field  technique  using  the  extra  reflection  and  the  diffusion  streak  associated  with 
twins  in  the  matrix  is  in  agreement  with  die  above  analysis  for  7"  faulting.  Kirman  and 
Warrington  [136]  believed  that  this  reaction  may  be  an  important  factor  in  the 
deformation  of  y"  dispersions. 

This  shearing  mechanism  by  partial  dislocations  can  promote  deformation  by 
twinning  [15].  The  above  dislocation  reaction  leaves  partial  dislocation  loops  around 
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precipitates  which  can  eventually  coalesce  under  the  influence  of  internal  stresses  to 
propagate  micro-twins  in  the  matrix.  Furthermore,  a  noteworthy  feature  of  the 
microstructure  of  the  solution  treated  as  well  as  deformed  alloy  was  that  dislocations  were 
invariably  arranged  in  planar  arrays.  Apparently  cross  slip  of  dislocations,  leading  to 
tangling  and  cell  formation,  was  almost  absent  [162].  Such  situadon  would  arise  if  the 
stacking  fault  energy  (SFE)  of  the  alloy  was  sufficiently  low  to  inhibit  cross  slip  [163]. 
The  low  SFE  structure  also  favors  the  formation  of  mechanical  twinning  in  FCC  metals 
[164,165].  Though  pure  nickel  has  a  high  SFE  (about  240  mJ/m2  [166]),  it  is  known  that 
alloying  additions  of  elements  like  Cr,  Mo,  A1  and  Ti  can  bring  down  the  value  of  the 
SFE  substantially  [165,167,168].  From  the  knowledge  of  the  various  elements  present  in 
this  alloy,  the  SFE  was  roughly  estimated  by  Sundararaman  et  al  [162]  to  be  75  mJ/mz 
and  by  Fournier  and  Pineau  to  be  SOmJ/m2,  which  is  lower  than  that  of  pure  nickel. 
Although  these  values  should  be  regarded  as  approximate  ones,  they  point  to  the  fact  that 
Alloy  718  is  one  of  the  intermediate  SFE  nickel  base  alloys.  A  common  feature  of  these 
alloys  is  that  dislocations  in  these  materials  are  arranged  in  planar  arrays  [165],  and 
del  irmation  by  twinning  has  effectively  been  observed  for  such  values  of  the  SFE  in  the 
case  where  high  tensile  stresses  were  applied  [164].  Moreover,  Stacking  faults  were  quite 
frequently  observed  within  y"  precipitates  in  deformed  as  well  as  undeformed  specimens. 
This  suggested  that  the  SFE  of  the  y"  phase  in  Alloy  718  is  quite  low.  However,  no 
further  work  has  been  done  for  determining  the  value  of  the  SFE  in  the  y"  phase. 

It  is  interesting  to  note  that  a  great  majority  of  deformation  twins  remained 
confined  within  the  y"  particles  and  did  not  propagate  into  the  matrix  [162].  This 


observation  was  suggestive  of  the  fact  that  the  value  of  the  critical  resolved  shear  stress 
for  twinning  in  the  y  matrix  was  higher  than  that  in  the  y"  precipitates.  As  pointed  out 
earlier,  the  SFE  of  the  y  phase  was  higher  than  that  of  the  7"  phase  and  possibly  this  was 
the  reason  why  the  propagation  of  deformation  twins  through  the  7  phase  was  relatively 
difficult. 

Sundararaman  and  associates  [162]  have  discussed  the  mode  of  dislocation- 
precipitate  interaction  with  reference  to  the  room  temperature  deformation  behavior  of 
Alloy  718.  They  observed  that  matrix  dislocations  often  appeared  in  pairs  and  the  pair 
of  dislocations  were  mainly  associated  with  7"  particles,  which  suggested  that  these 
dislocations  are  superdislocafions  that  glided  through  the  ordered  precipitates.  When  a 
pair  of  l/2[fl0]  dislocations  pass  through  7"  particles  and  emerge  in  the  matrix,  no  order 
faults  are  produced.  Thus  the  paired  dislocations  observed,  they  believed,  was 
corresponding  to  this  situation.  They  found  that  the  7"  particle  size  is  the  parameter 
which  controls  a  slip-to-twin  transition,  i.e.  the  deformation  through  the  7"  precipitates 
occurs  by  the  passage  of  a  group  of  dislocation  which  restores  order  in  them  for  particle 
mean  radius  smaller  than  IOOA,  while  for  precipitates  radii  larger  than  100A,  precipitate 
shearing  occurs  by  deformation  twinning  of  precipitates.  In  their  study,  the  bypassing 
mechanism  was  not  observed  for  the  range  of  7"  particle  size  (up  to  600A).  Contrast  to 
the  work  of  Sundararaman  et  al,  Wilson  [16]  has  observed  that  dislocations  were 
entangled  with  7"  particles  when  specimen  tested  at  760°C,  and  this  indicated  that 
dislocations  bypassed  the  precipitates.  The  discrepancy  of  their  results  might  related  to 
the  thermal  activation  by  allowing  dislocations  to  cross  slip  and  climb  out  of  their  original 
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slip  planes  [169]. 
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1.4  Strengthening  Mechanisms  in  Alloy  718 

Oblak  et  al  [134]  and  Chaturvedi  et  al  [170]  have  studied  the  strengthening 
mechanisms  in  Alloy  718.  Oblak  et  al  [134]  have  calculated  the  critical  resolved  shear 
stresses  based  on  the  order  strengthening  model  and  coherency  strengthening  model, 
respectively.  By  comparing  these  results,  they  found  three  major  reasons  why  the 
coherency  strengthening  appears  to  be  predominant  rather  than  the  order  strengthening  as 
follows: 

a)  the  increase  in  critical  resolved  shear  stress  predicted  by  the  coherency 
strengthening  model  exceeds  that  expected  for  order  strengthening  and 
closely  matches  the  experimental  results: 

b)  the  temperature  dependence  of  the  critical  resolved  shear  stress  follows 
that  of  the  shear  modulus  as  predicted  by  the  coherency  strengthening 
model; 

c)  if  the  order  strengthening  mechanism  were  predominant,  it  would  suggest 
that  the  alloy  should  deform  at  relatively  low  stressed  by  the  motion  of 
dislocation  quadruplets.  However,  this  is  inconsistent  with  the  observation 
of  dislocation  pairs  rather  than  quadruplets  in  the  material. 

Thus,  they  suggested  that  the  strength  of  Alloy  718  depends  primarily  upon  the  coherency 
strengthening  arising  from  the  tetragonally-distorted  D022  y"  precipitate. 
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It  worth  mentioning  that,  after  calculating  the  antiphase  boundary  energy  (APBE) 
of  the  y*  phase,  they  found  the  value  of  this  APBE  was  surprisingly  low,  which  strongly 
suggested  that  the  contribution  of  the  7'  phase  to  the  strength  of  Alloy  718  is  minimal 
and  can  be  ignored, 

The  work  of  Chaturvedi  and  Han  [170]  confu,»’  i  that  the  major  strengthening  is 
provided  by  the  coherency  hardening  mechanism  by  the  7"  precipitate,  and  they 
considered  the  strengthening  contribution  by  the  7’  precipitate  particles,  although 
relatively  small  amount,  either  by  the  coherency  strengthening  mechanism  or  by  the  order 
hardening  mechanism. 


APPENDIX  II 


METHODS  OF  EXPERIMENT 


11.1  Specimens  Specifications 

The  compact  tension  (CT)  specimens  have  been  used  in  this  research  program. 
The  specimen  specifications  are  following  ASTM  standard  (E647-88a)  and  are  shown  in 
Fig.  A2.1. 

Fig.  A2.2  shows  the  microstructural  features  of  testing  materials  and  the  average 
grain  size  is  about  130  |im. 

11.2  Apparatus 

An  automated  test  system  based  on  a  microcomputer  was  used  in  conducting  the 
experiments.  All  tests  in  air  was  carried  out  on  a  MTS  mode  810  servo-hydraulic  test 
station  with  a  250  KN  load  cell  and  a  MTS  model  880  servo-hydraulic  test  station  with 
a  100  KN  load  cell  was  employed  when  tested  in  vacuum  chamber.  Inside  the  vacuum 
chamber,  the  elevated  temperature  is  reached  by  using  a  10  KHz  induction  heating 
system.  During  crack  growth  tests,  the  degree  of  vacuum  is  monotored  with  a  residual 
gas  analyzer  and  the  testing  vacuum  is  better  than  5  x  10'®  torr  total  pressure  which 
consists  primarily  of  nitrogen  and  oxygen  with  partial  pressure  less  than  10‘9  torr.  The 
vacuum  testing  system  is  shown  in  Fig.  A2.3. 

The  digital  PDP  11/34  microcomputer  through  the  MTS  433  processor  interface 
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Fig.  A2.2  The  microstructure  of  the  testing  material. 
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was  used  in  real  time  for  feedback  control  of  the  test  parameters  and  to  acquire  data.  An 
IBM  PC-AT  computer  was  utilized  to  reduce  data.  The  specimens  were  heated  with  a 
resistance  two-zone  furnace  when  tested  in  air  and  with  a  10  KHz  induction  heating  unit 
when  tested  in  vacuum.  Chromel/Alumel  thermocouple  was  used  to  measure  the 
temperature,  and  the  temperature  at  a  point  on  the  specimen  was  controlled  at  650°C  ± 
2°C  by  a  CN1204  temperature  controller. 

II.3  Techniques  of  Measurement  of  Crack  Length 
II.3.1  The  Elastic  Compliance  Method 

For  this  method,  crack  length  is  determined  indirectly  from  elastic  compliance 
measurements  using  crack  mouth  opening  displacement  (CMOD)  data  obtained  by  a 
extensometer  having  mechanical  extension  (shown  in  Fig.  A2.4).  During  a  fatigue  crack 
growth  test,  CMOD  and  load  data  are  acquired  simultaneously  for  the  unloading  part  of 
the  histerics  loop  by  the  computer,  then  a  straight  line  is  fitted  according  to  the  least- 
square  error  minimization  procedure  to  this  data  set,  the  load  data  of  unloading  part  is 
chosen  from  45%Pm„  to  95%Pmtx  in  order  to  obtain  good  linear  relationship  (shown  in 
Fig.  A2.5).  Thus,  the  elastic  compliance,  C,  can  be  obtain  as  a  slop  of  this  straight  line: 

C  .  ™.OD  (2.1) 

p 

By  using  the  initial  crack  length  a0  and  elastic  compliance  C,  the  apparent  Young’s 
modulus  E  can  be  estimated  as  [171]: 
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Fig.  A2.4  Extensometer  with  mechanical  extension  used  in  the  experiments. 
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where  B  is  the  thickness  of  the  specimen,  b  is  the  ratio  of  a^/w.  Now,  the  current  crack 
length,  a,  can  be  calculated  by  the  following  expression  [171]: 


a=  W(  1.001  -4.6695U  + 18.46U2-236.82U3 

+  1214.9U*-2143.6U5)  (2-3> 

where 

•U  =  (v^BEC"  +l)  * 

Thus,  by  continually  measuring  the  current  elastic  compliance,  C,  the  current  crack  length, 
a,  can  be  determined  according  to  eq.  (A2.3). 


II.2.2  The  D.  C.  Potential-Drop  Method 

The  D.  C.  potential-drop  method  for  measuring  crack  growth  is  relied  on  the 
increase  in  electrical  resistance  of  the  specimen  when  its  cross-sectional  area  is  reduced 
[172].  The  basis  of  this  method  is  to  pass  a  known  direct  current  generated  by  a  stable 
power  supply  through  a  test  specimen  and  to  monitor  the  change  in  potential  drop 
between  one  active  probe  pair  which  span  a  crack  as  the  crack  growth  occurs.  These 
probes  will  react  to  changes  in  crack  length,  temperature,  current  and  deformations.  A 
reference  probe  pair  is  located  at  a  position  remote  to  the  crack.  This  probe  pair  will 
react  in  the  same  manner  to  all  variables  except  crack  length.  As  the  crack  grows,  the 
increase  in  electrical  resistance  of  the  specimen  will  result  in  the  potential  drop  to  rise, 
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thus  the  measurements  of  the  potential  drop  across  the  crack  region  can  be  related  by 
experimentally  determined  calibration  curve  to  crack  length.  The  D.  C.  current  applied 
for  CT12.5  is  5  A  while  for  CT31.5  is  12  A.  Fig.  A2.6  shows  a  schematic  of  D.  C. 
potential  drop  system. 

An  thermal  electrometric  force  (EMF)  due  to  thermocouple  action  may  be 
generated  at  the  junction  of  the  potential  probe  and  the  specimen.  If  the  junction 
characteristics  of  the  both  potential  probes  are  identical,  the  net  EMF  across  the  probes 
will  be  zero.  However,  in  practice  there  is  always  a  contribution  of  EMF  to  the  output 
voltage  which  should  be  accounted  for.  In  order  to  eliminate  the  effect  of  EMF,  the 
corrected  voltage  V#  was  used  which  is  as  follows 


V«V-V0  (2.4) 

C  0 

where  V  is  the  voltage  measurement  when  current  is  on,  while  V0  is  the  voltage 
measurement  when  current  is  of,  both  are  measured  when  current  is  stabilized.  The 
current  configuration  during  cycles  is  shown  in  Fig.  A2.7.  Furthermore,  to  minimize  the 
effects  of  the  variations  of  alloy  resistivity,  temperature  and  applied  current,  the 
normalization  of  V0l  and  Vo2  was  utilized  here: 


V  V  -V 

vcl  _  v  l  V01 

v  v  -v 

v  c2  2  Y  02 


(2,5) 


where  Vol  is  the  corrected  voltage  for  active  probe  pair,  and  Ve2  is  the  corrected  voltage 
for  reference  probe  pair.  Thus,  the  ratio  V0l/V02  versus  the  number  of  cycle,  N,  is 
recorded  on  the  computer  disk,  and  this  set  of  data  can  be  converted  to  crack  length 
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Fig.  A2.6  Schematic  of  the  D.C.  potential  drop  system. 
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versus  N  by  experimentally  determined  calibrated  relationship  between  the  ratio  V0,/Vc2 
and  crack  length  a. 


II.2.3  The  Method  of  Obtaining  da/dN  and  AK 

From  D,  C.  potential  drop  method,  a  data  set  of  N  versus  Vcl/Vo2  is  obtained,  so 
is  has  to  be  converted  to  a  data  set  of  N  versus  a,  this  can  be  done  by  using  calibration 
curve  experimentally  determined  from  high  temperature  test  (650°C),  for  CT12.5  this 
calibrated  relationship  is  as  follows: 


-1  =  1.0986-0.03799 

-0.02092 

rv  T 

vol 

W 

lv-J 

vr 

’  v2 

1  J 

-0.001976 


f  N» 

V 

vct 

V  1 J 


+  0.0004657 


f  V 

vcl 

V.. 

V  61 J 


(2.6) 


It  should  be  mendoncd  that  since  initial  crack  length,  temperature,  crack  front 
configuration  are  not  exactly  the  same  for  calibration  test  and  real  tests,  thus  for 
relationship  of  (a/W)  versus  (V0,/Vcj),  there  exits  a  deviation  from  calibration  test  curve 
and  real  test  curve,  schematically  shown  in  Fig.  A2.8.  In  general,  they  are  related  as 
follows: 


g(v)  ■  f(v)  +  h(v)  (2*7) 

where  v  »  Vel/Va2,  g(v)  represents  the  real  test  curve,  f(v)  represents  the  calibration  test 
curve,  and  h(v)  is  a  certain  function  of  v,  for  the  sake  of  convenience  h(v)  usually  is 
taken  as  polynomial  form,  i.e. 
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g(v)  =  f(v)  +  Aj  +  A2v  +A3v2+  •••  (2.8) 

For  any  tests,  it  is  certain  that  the  initial  reading  (a0,  v0)  and  the  final  reading  (af, 
vf)  can  be  measured,  thus  for  simplicity,  h(v)  is  often  chosen  as  a  linear  form: 

h(v)  =  Aj+A2v  (2.9) 

so  we  have 


g(v0)=f(v0)+A}+A2v0  = 


g(vf)  »f(vf)+A1  +  A2vf  = 


(  \ 

ao 

\^ } 

f  \ 

h. 

w 

V  J 


(2.10) 


by  solving  eq.  (A2.ll)  simultaneously,  A,  and  A2  can  be  obtained,  and  the  relationship 
g(v),  i.e. 


g(v)  *f(v)  +  Aj  +  AjV  (2.11) 

can  be  used  to  convert  Vj/V2  to  crack  length  a.  More  readings,  except  initial  and  final 
readings,  obtained  during  test  will  certainly  benefit  to  increase  the  accuracy  of  this 
converting  process.  The  typical  curve  of  crack  length  versus  number  of  cycles  is  shown 
in  Fig.  A2.9. 

For  data  set  (N,a),  the  total  number  of  data  points  usually  is  around  serval 
hundreds.  However,  based  on  experience  of  analyzing  potential  drop  data  sets  [173],  it 
is  desirable  to  reduce  the  total  number  of  data  points  to  about  100.  The  so-called  block 
averaging  approach  was  used  in  the  study,  and  the  computer  program  is  listed  in 
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Fig.  A2.9  The  curve  of  crack  length  versus  number  of  cycle  obtained  from  D.C.  poten^al  drop  method. 


Appendix  IV. 


From  constant  AK  tests,  the  data  set  is  already  in  the  form  of  N  versus  a.  Thus 
in  this  case  the  data  set  just  need  to  be  smoothed,  and  then  the  linear  regression  analysis 
for  this  data  set  can  be  applied  to  obtain  the  slope  which  is  the  crack  growth  rate  da/dN. 

The  seven  point  incremental  polynomial  technique  [174]  has  been  used  to  obtain 
fatigue  crack  growth  rate  da/dN,  and  the  stress  intensity  factor  range  AK  is  calculated 
according  to  the  following  expression  [174]: 

AK  »  -J— (0.886  +  4.64  a 

Bv/W"  (l-a)u  (2.12) 

-13.32aM4.72a3- 5.5a4),  for  as0.2 

where  a  «  a/W. 

Finally,  in  order  to  maitain  crack  growth  in  the  plane  strain  condition,  for  CT 
specimen  the  following  requirement  should  be  fulfilled  [174]: 

W-a&iiS^L  (2.13) 

* 

or 

AK£  Jit  (W-a)  (2-14) 

2 

where  aYc  is  0.2%  offset  cyclic  yield  strength,  R  is  the  load  ratio. 
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DERIVATION  OF  THE  EXPRESSION  FOR  CALCULATING 
OXIDATOIN  CHARACTERISTIC  PARAMETER  (3 


During  the  grain  boundary  oxidation  process,  the  chromium,  iron  and  nickel  oxides 
could  form  as  follows: 


2Cr  +  ^.02  ~  Cr203 


Fe  =  i.CL  =  FeO 
2  i 


(A3.1) 


(A3.2) 


Denoting  m02°r,  m02F*  and  m^’1  as  the  amount  of  oxygen  required  to  form  oxides 
Cr203,  FeO  and  NiO,  respectively.  We  assume  that 


m 


Cr 


Ft 


0, 


mn  +  m0  +  m( 


.  Nl 

0. 


(A3.4) 


From  eqs,(A3.2)  and  (A3. 3),  the  linkage  between  m^0,  m02Fil  and  m02NI  can  be 
established  as  follows.  According  to  eq.(A3.2)  we  know  that  1  mole  of  Fe  comsumes  1/2 
mole  of  02  to  form  FeO,  i.e.  xF,  moles  of  Fe  need  x02F*  moles  of  02  to  form  FeO: 
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(A3 .5) 


the  value  of  moles  for  a  substance  is  the  ratio  of  the  amount  of  the  subatanee  to  the 


atomic  or  molecular  weight  of  the  substance.  Thus  according  to  this  definition,  we  have 


where  mFi  is  the  amount  of  iron  required  to  form  the  iron  oxide  FeO,  and  MFe  is  the 
atomic  weight  of  iron.  Therefore,  we  have 


mFe  0  m0, 

MPe-  M0i 


or 


mF#  = 


2M,  p. 

m0, 


3.491  mj; 


(A3.6) 


Similarly,  from  eq.(A3.3)  we  can  get 


^MN,  Nl 

mM,  =  . . .  Illg 

N1  Mn 


N1 


3.6694  m^ 


Now,  if  eqs.(A3.6)  and  (A3.7)  are  combined  together,  we  can  obtain 


(A3.7) 
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(A3.8) 


rn0 

— —  =  0.9513 — I 
m 


m  ni 
m0, 


For  the  sake  of  convenience,  we  denote 


then  eq.(A3.8)  could  be  written  as 

0.9513  m£-Ymo!a0 

By  solving  eqs.(A3,4)  and  (A3.9)  simultaneously,  one  can  obtain: 


nu  -  m 


0.9513  +  y 


and 


m 


NI 

°a 


a  0.9513 


mo, 

0.9513  +  V 


From  eq.(A3.1),  the  relationship  between  m^0,  and  could  be  determined 


=  0.4615  mCr 

Substituting  eq.(A3.12)  into  eqs.(A3.10)  and  (A3. 11),  we  have 

Fe  mO,  “  °'4615 

m  q  -  y  — - - 

0>  '  0.9513 +Y 


and 


(A3.9) 


(A3.10) 


(A3.ll) 

as 

(A3.12) 


(A3.13) 
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m^1  =  0.9513 


mn  -  0.4615  m 

ua 

0.9513+7 


Cr 


(A3.14) 


For  the  same  reason  as  in  eq.(A3.6),  we  have 


m, 


FeO 


M 


FeO 


or 


•  _  «  ^FeO  Fe 

ITlfj  a  ®  '  ■  ■  in 

Fe0  M 
m0j 


O. 


-4.49  m  % 


(A3.15) 


Now,  combining  eqs.(A3.14)  and  (A3. 15),  one  obtains 


mFe0  =  4.49  7  —4 


mQ  -0.461 5  mCr 


0.9513+7 


In  the  same  way,  we  can  get 


mN!o  =  4.442  7 


m0j  -  0.461 5  mCr 
0.9513+7 


and 


m 


CfsO, 


1 .462  mCr 


(A3.16) 


(A3.17) 


(A3.18) 


m 


Since  parameter  (3  is  defined  as 
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(A3.19) 


(5  = 


mCr,0, 


m 


ox 


m, 


Cr,0, 


mCrjO,  +  mPeO  +  mNiO 


by  substituting  eqs.(A3.16),  (A3. 17)  and  (A3. 18)  into  (A3. 19),  the  expression  of  parameter 
P  in  terms  of  m^,  and  y  can  be  obtained  as  follows: 


P- 


CtmCr 


C2  mo,  -  C3  mCf 


(A3.20) 


where 


C,  -  1.462(0.951  +  Y) ,  C2  =  4.442  +  4.498 Y  ,  C3  -  0.659  +  0.610y 
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APPENDIX  IV 


COMPUTER  CODE  FOR  PREDICTING  da/dN  BASED  ON 
THE  PROPOSED  MODEL 


Q  **+*  id**  id***  #1(1*  *<  it***  ****>(i***>ti>i»tM  id**  ***>!<>♦« 

C  MODEL.FOR  IS  A  PORGRAM  THAT  SIMULATES  THE  FATIGUE  CRACK 
C  GROWTH  RATE  BASED  ON  THE  PROPOSED  MODEL 

^  |||  ^  j||  ^  |||  ^  ^  |||  |||  ||^  |||  ||^  ^  j||  ||^  j|^  |||  |||  j||  ^  |||  iff  ||f  j||  ||| j|g  ||^  |||  |||  ||^  |||  j||  ||^  ^  )||  ||f  j||  |||  ^  ^ ||| 


DIMENSION  DK(50),RES (50), GA  1(50), GB 1  (50), GC 1  (50), GD 1  (50) 
COMMON/AA1/  ETAJ3TA1.BTA 
COMMON/AA2/  DTK,R0,D1,R,TH 
COMMON/AA3/  DLPDS 
COMMON/AA4/GA,GB,GC,GD 
C  INPUT  r*,d, hold-time, frequencyDELTA-K,Ratio-R 

OPEN  (20,FILE=  ’DA  1  D  AT’, STATUS = ’  OLD  ’ ) 
READ(2O,*)R0,DD,TO,FQ,R,STEP 

C  INPUT  COEFFICIENTS  OF  GRAIN  BOUNDARY  DUCTILITY 
C  INPUT  A  SET  OF  STRESS  INTENSITY  FACTOR  RANGES  DELTA-K 
OPEN(22,FILE=’DA2.DAT\STATUS=’OLD’) 
OPEN^.FnJE^’DASDAT’.STATUSs’OLD’) 

DO  9  1=1,50 

READ(22,*,END=1 1 )  GA  1  (I),GB  1  (I),GC  1  (I),GD  1  (I) 
READ(24,*,END=11)  DK(I) 

9  CONTINUE 
11  NPTS=I-1 
RR0=R0 

DO  12  JY=1,NPTS 
WRITE(6,300)  NPTS.JY 
300  F0RMAT(1X,’NPTS=’,I2,2X,’JY=\I2) 

DTK=DK(JY) 

GA=GA1(JY) 

GB=GB1(JY) 

GC=GC1(JY) 

GD=GD1(JY) 

R0=RR0+(JY-1)*STEP 

RI0=R0 
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c 

3 


501 

105 

34 


Il»l 

JJ=0 

WRITE(6,2) 

2  FORMAT (4X, ’  r*  (uM)  d(uM)  Th(scc.)  fq(Hz)  DTK(MPaml/2)  R’) 

WRITE(6,2 1)R0*DD,THJFQ,DTK,R 
21  FORMAT(6F9.3) 

DLP=8.3 14*923,0 
RR=l./FQ-5.0 
1F(RR.GT.0.0)  THEN 
GOTO  3 
ELSE 
TR-TH 
GO  TO  105 
END  IF 

CALCULATION  OF  THE  CYCLE  EFFECTIVE  OXIDATION  TIME  t* 
TR=TH+(1/FQ*IFIX(RR/ABS(RR))) 

ETA  1«DD/SQRT(TR*  1 .2E-6) 

WRITE(6,501)ETA1 
FORMAT(2X,F12.6) 

ETA«RO/SQRT(TR*  1 .2F.-6) 

WR1TE(6,34)  II 

FORMAT(5X,  ’The  Number  of  Iteration  II  :\I3) 

CALL  ENGY 

CALL  SIMP2(2,1 .0 JDS,  1  .E»6,CR) 

CALL  SIMP2(  1 , 1 .0J9 1 , 1  .E-6,02) 

C02I«ETA*02*0.28209479 
CCRI=ETA1  *CR*0.28209479 
C  THE  AMOUNT  OF  OXYGEN 

C02«C02I*97.181 

C  THE  AMOUNT  OF  CHROMIUM 

CCR=CCRI*  1 .4988E6 
WRITE(6,10)  C02.CCR 
UV=3.1273*CO2-0.4433*CCR 
IF(UV.EQ.0.)  THEN 
WRITE(6,71) 

71  FORMAT(3X,’The  model  does  not  apply  for  this  case’) 

STOP 

ELSE 

CONTINUE 
END  IF 
BTA=CCR/UV 
WRITE(6,99)BTA 
99  FORMAT(5X,’BTA»’,F10.6) 

CALL  CMO(DR,CP) 
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EF(ABS(CP).LT.3.5E-2.AND.ABS(CP).GT.l.E-2)  THEN 
DRR=ABS(DR) 

END  EF 

WRITE(6,51)  DRR 
51  F0RMAT(5X,  ’DRR=’  ,E20. 10) 

ADR=ABS(DR)/FL0AT(3*I1) 

DF(CP.GT.O.O)  THEN 
R01=R01+ADR 
ELSE 

R01=R01-ADR 
END  IF 

IF(ABS(CP).LT.5.E“3)  THEN 
GOTO  202 
END  IF 

WRrTEf6,42)R01 

42  FORMAT(3X,’R0='E20.10) 

IF(I1.GT.15)  THEN 
GOTO  202 
END  IF 
R0»R01*1.E6 
Il-Il+l 
GOTO  105 
202  CONTINUE 

EF(JJ.EQ.l)  THEN 
GOTO  25 
END  IF 
RI^ROl 

WRITE(6,37)  RL 

37  FORMAT  (5X,  ’RL=  ’  ,E20. 1 0) 

WRITE(6,28) 

28  FORMAT(10X, ’Iteration  from  the  other  side’) 

DP=RL*1.E6~RI0 
WRITE(6,63)  DP 

63  FORMAT  (5X,  ’D  P=  ’  ,E20. 1 0) 

DF(DP.GT.0.0)  THEN 
R0-RL+DRR 
ELSE 

R0=RL-DRR 
END  IF 

WRITE(6,64)  R0 

64  FORMAT(5X, ’inial  value  iO  for  other  sidel:’E20.10) 

R0=R0*1.E6 

11=1 

JJ=1 
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n  n 


GOTO  105 
25  RKaROl 

C  THE  FINAL  RESULT  r* 

R02»(RL+RK)  *0.5  *  1 E-6 
102  WRrTE(6,35)  R02 

35  FORM  AT(3X,  ’  da/dN»  ’  ,E20. 1 0) 

RES(JY)=R02 

10  FORMAT(2X,’CO2=\E20.10, 1X,’CCR=’ JE20. 10) 

12  CONTINUE 

OPEN^S.FTLEs'RESULT.DAT’.STATUS^’NEW’) 

WRITE(25,31) 

31  FORMAT(4X, ’Fatigue  Crack  Growth  Rate«da/dN  (mm/cyc)’) 
WRHE(25,32)  (RES(J),J»1,NPTS) 

32  FORMAT(8X,E20. 10) 

STOP 

END 


SUBROUTINE  SIMP2(J,A,B,EPS,SUM) 

DIMENSION  F(2,30),FM(2,30),E(2,30),KRTN(30) 

SUM-0.0 

T-1.0 

ABSA-1.0 

EST-1.0 

DA=B-A 

FA=FCT(A,J) 

FB=FCT(B,J) 

RA1»0.5*(A+B) 

FP-4.0*FCT(RA1,J) 

X»A 
L=0 
1  K*i 
L-L+l 
T-T*1.7 
DX=D  A/3.0 
SX-DX/6.0 
RA2=X+0.5*DX 
FM1=4.0*FCT(RA2,J) 

RA3-X+DX 

F1=FCT(RA3,J) 

FM(1,L)=FP 

RA4=X+2.0*DX 

RA5=X+2.5*DX 

F(1L)“FCT(RA4,J) 


on  s  no 


FM(2,L)«4.0*FCT(RA5,J) 

F(2,L)=FB 

E1-SX*(FA+FM1+F1) 

E(1,L)®SX*(F1+FP+F(1,L)) 

E(2,L)«SX*(F(1,L)+FM(2,L)+FB) 

S=E1+E(1,L)+E(2,L) 

ABSA=ABSA-ABS(EST)+ABS(E1)+ABS(E(1,L)) 
*  +ABS(E(2,L)) 

IF(EST.EQ.  1.0)  GOTO  5 
IF(T*AB.S(EST-S).LE.EPS*ABSA)  GO  TO  2 
IF(L.LT.30)  GO  TO  5 

2  SUM-SUM+S 

3  L-L-l 
T-T/1.7 
K-KRTN(L) 

DX«DX*3.0 
IF(K.EQ.3)  DF(L>1)  6,6,3 
EST=E(K,L) 

EP=FM(K,L) 

FA-FB 

FB«F(K,L) 

K-K+l 

X«X+DA 

4  DA-DX 
KRTN(L)»K 
GOTO  1 

5  EST=E1 
FP=FM1 
FB=F1 
GO  TO  4 

6  RETURN 
END 


FUNCTION  FCT(X,J) 

COMMON/AA1/  ETA3TA1.BTA 
IF(J.EQ.2)  GO  TO  10 

FCT»EXP(-0.25*ETA*ETA/X)/(X*SQRT(X)) 

RETURN 

I  FCT«EXP(-0.25*ETA1*ETA1/X)/(X*SQRT(X)) 

RETURN 

END 

CALCULATION  OF  INELASTIC  STRAIN  ENERGY  Wp, GRAIN 
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BOUNDARY  DIFFUSIVITY  Dg  AND  DISLOCATION  DIFFUSIV1TY  Dd 

SUBROUTINE  ENGY 
COMMON/AA2/  DTK,R0,D1,R,TH 
COMMON/AA3/  DLP,DS 
R01»R0*lE-6 
DWP-0.0 

IF(TH.EQ.0.0)  GOTO  4 
CONVERTING  DELTA-K  TO  Kmin 
EK=R*DTK/(  1  .-R) 

DWC»EK*EK/R01 

CREEP  PART  OF  Wp 

DWP»DWC*2.881542E-7 

D  WF-DTK*DTK/RO  1 

FATIGUE  PART  OF  Wp 

DWPl«DWF*2.9762E-9 

TOTAL  INELASTIC  STRAIN  ENERGY 

DWP2»D  WP+DWP I 

WRITE(6,7)  DWP2 

FORMAT(3X, ’inelastic  strain  energy  density  Wp  \E20.10) 

G1»(5.61144*DWP2+2.17105)**2 

G2--14.98601349-25.8863*EXP(0.520202/G1) 

GRAIN  BOUNDARY  DIFFUSIVITY  Dg 

DG-EXP(G2)*1. 47736 1802E6 

THE  RATIO  OF  Dg/Dm 

D1-DG/1.2E-18 

WRrTE(6,3)Dl,DG 

FORMAT(3X, ’The  Rado  of  Dg/Dm:’ (2E20. 10) 

DIFFUSIVITY  Dd  FOR  Cr 
WEl«(35.7822*DWP2-37.572)**2 
WE2»32.007*EXP(0.37 1 1 87/WE  1) 

DDL»EXP(-  16.69870073-WE2)*  1.43000365779E7 

DS=DDU1.2E-18 

WRrTE(6.6)DS.DDL 

FORMAT(3X, 'The  Ratio  of  Dd/Dm:’,2E20.10) 

RETURN 

END 


CALCULATION  OF  INELASTIC  STRAIN  AND  GRAIN  BOUNDARY 
DUCTILITY 

SUBROUTINE  CMO(DR,CP) 

COMMON/AA2/DTK,RO,D  1  ,R,TH 
COMMON/AA 1/ETA3TA 1  ,BTA 


t 


C0MM0N/AA4/GA.GB  ,GC,GD 
C  CONVERTING  SECONDS  TO  HOURS 
TH1 -TH/3600. 

RPO-RO*  IE-6 

C  GRAIN  BOUNDARY  DUCTILITY 

EC=GA*ATAN(-GB*BTA*BTA+GC)+GD 

WRITfi(6,89)EC 

89  FORMAT  (5X,  ’EC- '  ,E20. 10) 

C  TOTAL  INELASTIC  STRAIN 
EMC*0,0 

IFCra.EQ.0.0)  GOTO  57 
AK«DTK*R/(1.-R) 

EMI  -TH 1 *  *0.0909090909 1 
EM2-AK**  1.81 8 1818182 
EM3  -RP0*  *  (-0.9090909 1 ) 

EMC«(EM1  *EM2*EM3)*8.34E-9 

57  EMF»((DTK*DTK/RPO)*,*0.909091)*6,202E-11 
EM-EMC+EMF 

WRITE(6,10)EM 

10  FORMAT  (3X,  ’EM- '  ,E  1 8 . 1 0) 

CP-(EM-EC)/EC 

AM-0.0 

EF(TH.EQ.0.0)  GOTO  58 
AM-EMl*EM2*8.34E-9 

58  BM1-DTK**  1.8 18182 
BM=»BM1*6.202E-1 1 

R 1  -(( AM+BM)/EC)  *  *  1 . 1 

C  THE  DIFFERENCE  BETWEEN  r*  AND  r-trying 
DR=RP0-R1 
RETURN 
END 
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APPENDIX  V 


COMPUTER  CODES  FOR  DATA  REDUCTION 


Q  **************************************************1**************** 

C  CT122.BAS  IS  A  PROGRAM  THAT  CONVERTS  A  SET  OF  DATA  N  VS 
C  V1/V2  TO  A  SET  OF  DATA  N  VS  A  FOR  CT12.5  SPECIMEN  BASED  ON 
C  TWO  MEASURED  PINTS 

0  ****  *41  *  *  *******  *  **  *  *  *  *  *  *  «  *****  4n«i  *  *  *****  ************  ***  *  *  if****  i(i  ** 


10  PRINT  "ENTER  THE  NAME  OF  THE  INPUT  FILE":INPUT  S$ 

20  PRINT  "ENTER  THE  NAME  OF  THE  OUTPUT  FmE":INPUT  T$ 

30  PRINT  "ENTER  V1/V2  (BEGINNING)  AND  A/W":INPUT  V1.AW1 
60  PRINT  "ENTER  V1/V2  (FINAL)  AND  A/W":INPUT  V2.AW2 
70  PRINT"ENTER  THE  WIDTH  OF  SPECIMEN"  :INPUT  W 
80  OPEN  S$+"DAT"  FOR  INPUT  AS  #1 
90  OPEN  T$+"DAT"  FOR  INPUT  AS  #2 
100  C 1 » 1 .007 :C2=-0.0357:C3*-0.0 1 3898 

1 10  AWC1-(C1+C2*V1+C3*V1A2):AWC2-(C1+C2*V2+C3*V2A2) 

1 20  A«((A  W2 1  -AW  1 )-( AWC2-A  W2))/(V2»  V 1) 

1 30  B=(AW1-AWC1)-A*V  1  :G  1»C1  +B:G2=C2+A:G3=C3 

140  INPUT  #1,  N,V 

150  A«W+(G1+G2*V+G3*VA2) 

160  PRINT  #2,  USING  "#######";N, 

170  PRINT  #2,  USING  "####.#####"; A 
180  EFEOF(1)=0  GOTO  140 
190  CLOSE  #l:CLOSE#2:END 


0  ******************************************************************** 
C  CT124.BAS  IS  A  PROGRAM  THAT  CONVERTS  A  SET  OF  DATA  N  VS 
C  V1/V2  TO  A  SET  OF  DATA  N  VS  A  FOR  CT12.5  SPECIMEN  BASED  ON 
C  FOUR  MEASURED  PINTS 

0  ******************************************************************** 


10  PRINT  "ENTER 

20  PRINT  "ENTER 

30  PRINT  "ENTER 


THE  NAME  OF  THE  INPUT  FILE": INPUT  SS 
THE  NAME  OF  THE  OUTPUT  FILE": INPUT  T$ 
V1/V2  (BEGINNING)  AND  A/W’’:INPUT  V1AW1 
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40  PRINT  "ENTER  V1/V2  (SECOND)  AND  A/WiINPUT  V2.AW2 
50  PRINT  "ENTER  Vim  (THIRD)  AND  A/W";INPUT  V3AW3 

60  PRINT  "ENTER  Vim  (FINAL)  AND  A/W":INPUT  V4.AW4 

70  PRINT'HNTER  THE  WIDTH  OF  SPECIMEN" :INPUT  W 
80  OPEN  S$+"DAT"  FOR  INPUT  AS  #1 
90  OPEN  T$+  ’DAT"  FOR  INPUT  AS  #2 
100  0=1 .Q0'7  :C2*-0.0357 :C3=-0.0 15898 

1 10  AWC1«(C1+C2*V1+C3*V1^2):AWC2-(C1+C2*V2+C3-  V2A2) 

120  AWC3=(C1+C2*V3+C3*V3A2):AWC4=(C1+C2*V4+C3*V4A2) 

130  U1*AW1-AWC1:U2»AW2-AWC2:U3=AW3-AWC3:U4=AW4-AWC4 

140  S1**V1-V2:Q1=V1A2-V2A2:F1=U1-U2 

150  S2=V3'V4:Q2»V3A2-V4A2:F2»U3-U4 

160  S3=V1"V3:Q3=V1A2-V3A2:F3“U1-U3 

170  E1«S2*Q1-S1*Q2:Y1=S2*F1-S1*F2 

180  E2=S1*Q3-S3*Q1:Y2=S1*F3-S3*F1 

190  C-V1/F.1 

240  B«(S2*F1-S2*Q1*C)/(S1-+S2) 

250  A«U1-B*V1-C  V  iA2 

260  GUC1+A:G°-C2+B:G3=C3+C 

270  INPUT  #1,  N,V 

280  A=W’\G1+G2*V+G3*VA2) 

290  PRINT  #2,  USING  "#######";N, 

300  PRINT  #2,  USING  "####.#####"; A 
310  IF  EOF(  1  )=0  GOTO  270 
320  CLOSE  #l:CLOSE  #2:END 


£  ******************************************************** ********** 

C  SMOOTH. BAS  IS  A  PROGRAM  THAT  SMOOTHS  THE  DATA  POINTS 

C  BASED  ON  SEVEN  POINT  MOVING  AVERAGE  METHOD 
£  ****************************************************************** 

10  DIM  X(3000),Y(3000),Y1(3000) 

20  1*0 

30  PRINT  "ENTER  THE,  NAME  OF  THE  INPUT  DATA  FILELJNPUT  S$ 

40  PRINT  "ENTER  THE  NAME  OF  THE  OUTPUT  DATA  FILE": INPUT  Z$ 

50  OPEN  Z$+".DAT"  FOR  OUTPUT  AS  #2 

60  OPEN  S$+".DAT"  FOR  INPUT  AS  #1 

70  INPUT  #1,  NA 

80  X(I)»N:Y(I)=A:I=I+1 

90  IF  FOF(IM)  GOTO  70 

100  CLOSE  #1:M=1-1:M1=M-3:PRINT  M 


110  FOR  J«3  TO  Ml 

120  Y 1  (J)=(Y  (J-3)+Y  (J-2)+Y  ( J- 1  )+Y  (J)+Y  (J+ 1  )+Y  (J+2)+ Y  (J +3))/7 

130  NEXT  J 

140  Yl(0)»(13*Y(0)+10*Y(l)+7*Y(2)+4*Y(3)+Y(4)-2*Y(5)-5*Y(6))/28 

150  Y1(1)=>(5*Y(0)+4*Y(1)+3*Y(2)+2*Y(3)+Y(4)-Y(6))/14 

160  Yl(2)»(7*Y(0)+6*Y(l)+5*Y(2)+4*Y(3)+3,(,Y(4)+2*Y(5)+Y(6))/28 

170  B  l=(Y(M-6)+2*Y(M-5)+3*Y(M-4)+4*Y(M-3)+5*Y(M-2)) 

1 80  Y1  (M-2)-(B  1  +6*Y(M-l)+7*Y(M))/28 

190  B2=(-Y(M-6)-2*Y(M-5)+Y(M-4)+4*Y(M-3)+7*Y(M-2)) 

200  Y1(M-1)=(B2+5*Y(M))/14 

210  B3=(-5*Y  (M-6)-2*Y  (M-5)+Y(M-4)+4*  Y(M-3)+7*  Y  (M-2)) 

220  Yl(M)«(B3+10*Y(M-l)+13*Y(M))/28 

230  FOR  J1=0  TO  M 

240  PRINT  #2,  USING  "#######";X(J1), 

250  PRINT  #2,  USING  "###.#####";' Y1(J1) 

260  NEXT  Jl:CLOSE  #2:END 


Q  ****#>M******ifc*****************>fi*********>M>H***>t<>li*>M*>H4<***>k**4<  +  **4< 

C  BLOCK.BAS  IS  A  PROGRAM  THAT  REDUCES  THE  NUMBER  OF  DATA 
C  POINTS  DOWN  TO  A  APPROPRIATE  NUMBER 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 


PRINT  ‘'ENTER  3 HE  NAME  OF  THE  INPUT  FILE":INPUT  Z$ 

PRINT  "ENTER  THE  NAME  OF  THE  OUTPUT  FILE":INPUT  S$ 

OPEN  S$  FOR  OUTPUT  AS  #1 

I=0:Vl=0:Il=O 

OPEN  Z$  FOR  INPUT  AS  #3 

INPUT  #3,  N,V:I1=11+1 

IF  EOF(3)=0  GOTO  60 

CLOSE  #3:PRINT  "TOTAL  DATA  POINTS' Ml 
PRINT  "THE  SMALLEST  BLOCK  SIZE  IS":INPUT  B2 
PRINT  "THE  INCREMENT  OF  EACH  BLOCK  IS":INPUT  D1 
Sl=»(l-2*B2/Dl):S2sSlA2:S3“8*Il/Dl:Pl=»(51+SQR(S2+S3))/2 
OPEN  Z$  FOR  INPUT  AS  #3 
I3«FIX(P1):BB=0 
IF  13  <0  GOTO  160 
B 1  =FIX(B2+(I3- 1 )  *D  1 ) 

INPUT  #3,  N,V 
I=I+1:V1=V1+V 
IF  1=1  ’l’HEN  M=N 
IF  EOF(3)=-l  GOTO  220 
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200  IF  I=B  1  GOTO  220 

210  GOTO  1 60:  V 1  =V  1/B 1  :GOTO  230 

220  Vl-Vl/I 

230  M 1  **N:  N 1  *FIX(M+M  1  )/2) 

240  PRINT  #1,  USING  "######";Nl, 

250  PRINT  #1 ,  USING  "####.######"; V 1 
260  I=0:V1=0:I3=I3*1:IF  I3<0  THEN  B  1=*1 

270  BB=BB+B1:IF  BB>I1  GOTO  290 
280  IF  EOF(3)=0  GOTO  140 
290  CLOSE  #l:CLOSE  #3:END 


0 

C  DADNDK.FOR  IS  A  PORGRAM  THAT  TAKES  N  VS  A  DATA  AND 
C  COMPUTES  DA/DN  USING  THE  SEVEN  POINT  INCREMETNAL 
C  POLYNOMIAL  TECHNIQUE  THAT  IS  RECOMMENDED  BY  THE  ASTM 
C  STANDARD  E647-88. 

C  **  ******iM*#*M  ***********  ill*********** 'fr*********************  ****** 

DIMENSION  A(900),  N(900),  BB(3),  DADN(900),  DELK(900),  ID(7) 
DIMENSION  AA(10),  NN(10) 

REAL  N 
REAL  NN 
INTEGER  QQ 
CHARACTER  NAME*  12 
1  FORMAT  (A) 

WRTTE(V)  ’ENTER  THE  NAME  OF  THE  DATA  FILE  CONTAINING 
*N/A  DATA’ 

READ(M)  NAME 
OPKN(l,FILE»NAME) 

READ  (1,*,END=10)  (N(I),A(I), 1*1, 10000) 

10  CONTINUE 

WRITE(*,*)  'ENTER  THE  MINIMUM  LOAD  AND  THE  MAXIMUM  LOAD' 
READ(*,*)  PMIN.PMAX 

WRITE(*,*)  ’ENTER  THE  THICKNESS  AND  THE  WIDTH  OF  THE  SPECIMEN’ 

READ(*,*)  B,W 

K=*0 

R=>PMIN/PMAX 
PP»PMAX-PMIN 
NPTS»I-1 
DO  11  I*1,NPTS 

11  CONTINUE 


WRITE(V)  ’ENTER  THE  NAME  OF  THE  OUTPUT  DATA  FILE’ 

READ(M)  NAME1 

OPEN(2,FrLE»NAMEl,STATUS=,NEW’) 

WRITE(V)  ’OBS.NO.  CYCLES  A(MEAS)  A(REG)  M.C.C.  DK  DA/DN’ 
DO  110  1=1,3 

D12=(A(I+1)-A(I))/(N(I+1)-N(I)) 

A1«(A(I+1)+A(I))*0.5 

T1=A1/W 

TPa(2.0+TD/(1.0-Tl)**(1.5) 

FTl=TP*(0.886-f4.64*Tl-13.32*Tl**2+14.72*Tl**3-5.6*Tl**4) 

DELK1-(31.624*(FT1*PP))/(D*SQRT(W)) 

C  WRITE(2,78)  I,N(I),A(I),DELK1,D12 
WRITE(2,79)  DELK1.D12 
WRITE(*,78)  I,N(I),A(I),DELK  1  ,D  12 
110  CONTINUE 
NPTS=NPTS-6 
DO  100  M.NPTS 
L=0 
K=K+1 
Kl=K+6 
DO  60  J=K,K1 
L=L+1 
AA(L)«A(J) 

NN(L)=N(J) 

60  CONTINUE 
C1=0,5*(NN(1)+NN(7)) 

C2»0.5*(NN(7)-NN(1)) 

SX=0 

SX2=0 

SX3=0 

SX4=0 

SY=0 

SYX=0 

SYX2=0 

DO  70  J-1,7 

X=(NN(J)-C1)/C2 

YY»AA(J) 

SX=SX+X 

SX2=SX2+X*X 

SX3=SX3+X*X*X 

SX4=SX4+X*X*X*X 

SY=SY+YY 

SYX=SYX+X*YY 

SYX2=>SYX2+YY*X*X 


70  CONTINUE 

DEN=7.0*(SX2*SX4-SX3**2)-SX*(SX*SX4-SX2*SX3)+SX2*(SX*SX3-SX2**2) 

T2=SY*(SX2*SX4-SX3**2)-SYX*(SX*SX4-SX2*SX3)+SYX2*(SX*SX3-SX2**2) 

BB(1)=T2/DEN  - 

T3=7.0*(SYX*SX4-SYX2*SX3)-SX*(SY*SX4-SYX2*SX2)+SX2*(SY*SX3  9 

*-SYX*SX2) 

BB(2)=T3/DEN 

T4=7.0*(SX2*SYX2-SX3*SYX)-SX*(SX*SYX2-SX3*SY)+SX2*(SX*SYX 

*-SX2*SY) 

BB(3)=T4/DF,N  # 

YB=SY/7.0 

RSS-0 

TSS=0 

DO  75  J=l,7 

X=(NN(J)-C1)/C2 

YHAT=BB(1)+*BB(2)*X+BB(3)*X**2  ft 

RSS=RSS+(AA(J)-YHAT)**2 
TSS=TSS+(AA(J)-YB)**2 
75  CONTINUE 

76FORMAT(I4,2X,F9.0,2X^8.2,2X,F8.2,2X,F8.4,2X,F8.2,2X,E15.5) 

78  FORMAT(I4,2X,F9.0,2X,F8.2,2X,8X,2X>8X,2X,F8.2,2X,E15.5) 

79  FORMAT(F8.2,2X,E15.5)  • 

R2=1.0-RSS/TSS 

DADN(I)=BB(2)/C2+2.0*BB(3)*(NN(4)-C1)/C2**2 

X=(NN(4)-C1)/C2 

AR=BB(1)+BB(2)*X+BB(3)*X**2 

QQ=I+3 

T=AR/W  • 

TP=(2.0+T)/(1.0-T)**(1.5) 

FT=TP*(0.886+4.64*T13.32*T**2+!4.72*T**3-5.6*T*'*4) 

DELK(I)=(3 1 .624*(FT*PP))/(B*SQRT(W)) 

C  WRITE(2,76)  OQ,N(QQ),A(QQ),AR,R2,DELK(I),DADN(I) 

WRITE(2,79)  DELK(I)JDADN(I) 

WRITE(*,76)  QQ,N(QQ),A(QQ),AR,R2DELK(I),DADN(I)  • 

100  CONTINUE 
J=NPTS+4 
K--NPTS+6 
DO  120  I=J,K 

D12=(A(I+1)-A(D)/(N(I+1)-N(U)  m 

A2=(A(I+1)+A(I))*0,5 

T2=A2/W 

TP=(2.0+T2)/(  1 .0-T2)**(  1 .5) 

FT=TP*(0.886-f4.64‘*,T2-13.32*T2**2+14.72*T2**3-5.6*T2**4) 

DELK2=(3 1 .624*(FT2*PP))/(B*SQRT(W)) 
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C  WRITE(2,78)  I,N(I)tA(I)fDELK2tD12 
WRITE(2,79)  DELK2.D12 
WRITE(*,78)  I,N(I)A(I),DELK2,D12 
120  CONTINUE 
STOP 
END 
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